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ABSTRACT 

Boyer,   Robert  Elston,   Ph.D.,   Purdue  University,   August   1972. 
Predicting  Pavement  Performance  Using  Time-Dependent  Transfer  Functions. 
Major  Professor:      Dr.   Milton  E.   Harr. 

An   approach  using  a  real  world  modeling  technique  was   developed  to 

predict  the  performance  of  flexible  pavement  systems  subjected  to 

dynamic   load  inputs.      "A  posteriori"  modeling  was  presented  in   contrast 

to  the   normally  accepted  "a  priori"   modeling  technique.      With  this 

approach,  the  initial  hypothesis  was: 

"that  time-dependent  functions   could  be  calculated  from  a 
pavement's   dynamic   input  and  output   deflection  responses;    that  these 
functions   could  be  obtained,    in  a  mathematical  sense,  without  the 
need  to  simulate   respective  material  performance  to  determine  values 
for  preselected  descriptors;    and  that  these  time-dependent    functions 
could  be  employed  to  predict   deflection  output   responses  when  the 
pavement  was   subjected  to  any   imposed  load". 

The   investigation  was   carried  out  by  extending  transfer  function 
theory  in   connection  with  a  decretized  convolution  procedure  to  define 
the  pavement's  time-dependent  transfer  functions.      The  time-dependent 
transfer  functions  were   calculated  from  data  measured  during   full  scale 
tests   in   a  service  environment   (three  pavement   cross-sections   and  three 
aircraft   of  different  weights   and  gear  configurations  were  used). 
Subsequently,    a  method  to  predict   deflection  output   responses  was 
developed  for  use   in   conjunction  with  the   calculated  time-dependent 
transfer  functions. 

It  was   concluded  that  the  characteristics  of  a  flexible  pavement 
could  be  represented  by  a  series  of  time-dependent  transfer  functions 
obtained  from  a  series  of  dynamic   tests   using  a  common  prime 
mover.      It  was  also  concluded  that  the  time-dependent  transfer 
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functions  contained  certain  isolated  descriptors  which  were  related 
to  the  aircraft  gear  configurations,   and  other  descriptors  which  were 
only  indicative  of  the  pavement's  behavior.      In  addition,   agreement 
was  achieved  between  the  calculated  and  measured  time-dependent 
deflection  functions  resultant  from  employing  the  procedures  in  the 
proposed  performance  evaluation  method. 


Pavement  evaluation  techniques  commonly  employed  in  practice 
follow  two  approaches.   The  more  common  method  (empirical  approach) 
is  to  base  the  evaluation  on  noted  changes  of  certain  material 
descriptors  of  the  pavement  system;  e.g.,  changes  in  strength  as 
measured  by  CBR  tests.   An  evaluation  is  then  made  by  correlating  the 
measured  descriptor  or  descriptors  with  similar  data  obtained  from 
previous  experience  on  a  similar  structure  in  a  similar  geographic 
locale.   These  procedures  consider  only  a  few  material  descriptors, 
and  offer  great  difficulty  in  extrapolating  to  greater  magnitudes,  or 
different  types  of  loadings  in  new  environments. 

Another  approach  used  to  evaluate  pavement  system  requires  the 
theoretical  postulation  of  an  analytical  model,  or  combinations  of 
models ,  followed  by  field  and/or  laboratory  tests  to  determine  the 
system  parameters.   This  approach  does  provide  a  mechanism  for 
extrapolating  to  greater  magnitudes  of  loadings,  and,  in  general,  is 
capable  of  incorporating  a  greater  number  of  material  descriptors; 
however,  the  degree  of  complexity  associated  with  a  theoretical 
approach  is  proportional  to  the  number  of  material  descriptors 
contained  in  the  selected  model.   A  shortcoming  of  this  technique  is 
measured  by  the  degree  to  which  the  preselected  model  actually 
represents  the  real  system. 

A  new  approach  to  the  pavement  evaluation  problem  is  advanced  by 
this  investigation.  The  objective  of  this  approach  is  to  overcome 
some  of  the  shortcomings  of  the  two  discussed  above.   In  this  method 
it  is  hypothesized  that  time-dependent  functions  can  be  calculated 
from  a  pavement's  input  and  output  deflection  responses;  that  these 


functions  (transfer  functions)  can  be  obtained,  in  a  mathematical  sense, 
without  the  need  to  simulate  respective  material  performance  to 
determine  values  for  preselected  descriptors;  and  that  these  time- 
dependent  functions  can  be  employed  to  predict  deflection  output 
responses  when  the  pavement  is  subjected  to  most  imposed  loads. 


REVIEW  OF  LITERATURE 

The  more  common  empirical  method  of  design  and/or  evaluation  of  a 
pavement  system  emerged  nearly  a  century  ago  in  the  form  of  rule -of - 
thumb  procedures.   Around  the  turn  of  the  century  prototype  test  roads 
were  used  to  determine  the  pavement  needs  to  accomodate  given  traffic; 
then,  between  World  Wars  I  and  II  highway  engineers  introduced  the 
classification  of  soils  to  correlate  performance  with  subgrade  types. 
During  World  War  II,  however,  a  new  problem  surfaced.   Engineers  were 
faced  with  the  need  to  predict  the  performance  of  pavement  systems  for 
greater  wheel  loads  and  frequencies  than  they  had  ever  before 
experienced  (references  2  and  3).   The  urgency  for  better  design 
techniques  led  to  consideration  of  available  theoretical  developments; 
thus,  a  rational  design  procedure  was  introduced  (reference  k) .      Since 
that  time  many  investigators  have  attempted  to  model  various  pavement 
components  and  systems  in  their  search  to  be  better  able  to  design 
and/or  to  evaluate  a  pavement. 

The  discussion  of  literature  which  follows  places  primary 
emphasis  on  solutions  derived  from  the  classical  theories  of  elasticity 
and  viscoelasticity.   The  nature  of  applied  loads  (i.e.,  static  and 
dynamic)  is  treated  separately.   Lastly,  some  developed  uses  of  the 
so-called  mechanistic  models  are  explored. 

Elastic  Solutions 


The  application  of  analytical  models  to  pavement  systems  has 
evolved  over  the  years  principally  through  the  use  of  the  theory  of 


elasticity.   Common  to  the  approaches  to  stress  analysis  problems 
using  the  theory  of  elasticity  are  the  assumptions  that  the  medium  is 
homogeneous  and  isotropic,  and  that  its  material  properties  can  be 
defined  by  linear  elastic  relationships. 

Static  Loading 

In  the  latter  part  of  the  19th  century,  J.  Boussinesq  developed  the 
stresses  and  deflections  within  a  semi-infinite  mass  under  a 
concentrated  vertical  force  by  applying  the  theory  of  elasticity. 
Beginning  in  1926,  Westergaard  (references  5,  6,  7,  and  8)  modeled  a 
pavement  as  a  flexible,  elastic  beam  supported  by  a  subgrade  consisting 
of  isolated  elastic  springs  (Winkler  foundation).  A  later  refinement 
was  introduced  by  Hogg  (reference  9)  who  used  a  flexible,  elastic 
plate  on  a  semi-infinite  foundation  to  model  a  pavement  system.  Harr 
and  Leonards  (reference  10)  first  employed  the  model  of  a  partially 
supported  slab.   Their  analysis  included  provisions  for  predicting 
stresses  and  deflections  produced  by  moisture  and  temperature  gradients 
within  concrete  slabs  in  addition  to  the  stresses  and  deflections 
induced  by  some  superimposed  vehicular  loadings.  Numerous  investigators 
(references  11,  12,  13,  1*+,  and  15)  have  presented  variations  to  refine 
Westergaard' s  model. 

In  19U3,  Burmister  (reference  16)  extended  the  Boussinesq  solution 
to  layered  systems.  His  models  (reference  IT)  proved  to  predict 
deflections  better  than  the  earlier  Boussinesq  model.  Subsequently, 
the  theory  of  elasticity  has  been  invoked  to  provide  a  multitude  of 
models  which  consider  additional  layers  and  variable  thicknesses  with 


consequent  multiplication  of  more  material  parameters  (references 
18,  19,  20,  21,  and  22). 

Dynamic  Loading 

Elastic  solutions  for  steady-state  dynamic  loadings  have  received 
considerable  attention  during  the  past  twenty  years.  In  1951,  Van  der 
Poel  (reference  23)  measured  velocities  of  elastic  waves  on  the  surface 
of  soils  and  roads.   Later,  in  195**,  Pickett  (reference  2k)   provided  a 
theoretical  explanation  for  the  propagation  of  waves  on  the  surface  of 
a  concrete  slab.   Only  in  the  past  decade,  however,  has  emphasis  been 
placed  on  the  so-called  surface  wave  technique.   The  procedure  involves 
measuring  the  wavelength  and  velocity  of  waves  at  the  surface  of  the 
pavement  under  test.   The  more  recent  theoretical  developments,  advanced 
by  Jones  (reference  25)  and  Jones  and  Thrower  (reference  26),  explain 
the  propagation  of  waves  in  a  multilayered  pavement  structure  in  terms 
of  the  elastic  properties  and  thicknesses  of  the  composite  layers.   In 
1967,  practical  applications  of  the  surface  wave  method  were  provided 
by  Jones,  Thrower,  and  Gatfield  (reference  27).  Good  correlation 
between  estimates  of  thickness  and  strength  of  in-situ  pavement 
sections  were  obtained  using  this  method.   Current  attention  is  being 
focused  toward  refinement  of  this  technique  (reference  27). 

Elastic-Viscoelastic  Solutions 

Recent  studies  (references  28,  29,  30,  31,  32,  3*+,  35,  36,  and  37) 
have  incorporated  time-dependent  attributes  to  a  portion  of  a  pavement 
system  by  using  either  a  viscoelastic  model  to  represent  a  pavement 


layer  over  a  Winkler  foundation,  or  an  elastic*  pavement  layer  over  a 
viscoelastic  foundation.   These  investigations  (for  both  static  and 
dynamic  loadings),  like  those  using  purely  elastic  methods,  assume 
component  materials  to  be  homogeneous  and  isotropic. 

Static  Loading 

Freudenthal  and  Lorsch  (reference  28)  analyzed  an  infinite, 
elastic  beam  on  three  different  types  of  linear  viscoelastic 
foundations  under  various  loadings.   Reissner  (reference  29)  provided 
the  solution  for  a  model  consisting  of  a  plate  resting  on  a 
viscoelastic  foundation  and  incorporated  shear  interactions  between 
the  plate  and  foundation  interface.   Hoskin  and  Lee  (reference  30) 
showed  that  for  some  special  but  important  cases  the  time-dependent 
problem  could  be  reduced  to  an  elastic  model  with  the  same  geometry  as 
the  parent  problem  ("correspondence  principle",  after  Lee  (reference 
31)).   Pister  and  Williams  (reference  32)  modified  Reissner's  (reference 
29)  concept  to  obtain  the  continuity  of  horizontal  displacements  at  the 
interface  between  an  elastic  plate  and  a  viscoelastic  foundation;  i.e., 
in  addition  to  permitting  horizontal  shearing  stresses  at  the  plate/ 
foundation  interface,  the  horizontal  displacements  were  also  preserved 
(not  taken  equal  to  zero).   Here,  also  the  elastic  plate  equation  was 
coupled  with  the  viscoelastic  modeling  of  the  foundation  through  the 
use  of  the  correspondence  principle.   A  Voigt  (spring  and  dashpot  in 
parallel)  model  was  employed  to  represent  the  foundation  behavior. 


*Here  and  in  subsequent  uses  the  term  "elastic"  is  understood  to  plead 
as  assuming  the  assumptions  of  the  classical  "theory  of  elasticity". 


Pister  and  Monismith  (reference  33)  used  a  3-element  viscoelastic 
model  to  represent  an  asphalt  concrete  beam  supported  on  a  Winkler 
foundation.   The  viscoelastic  solution  was  obtained  from  the  classical 
solution  of  an  elastic  beam  on  a  Winkler  foundation  through  the  use  of 
the  correspondence  principle.   Monismith  and  Secor  (reference  3k)    later 
employed  a  ^-element  viscoelastic  model  to  represent  an  asphalt  concrete 
layer  supported  on  a  Winkler  foundation.   In  contrast  to  the  earlier 
work  of  Pister  and  Monismith  (reference  3b),   an  additional  free-dashpot 
was  added  to  the  model  to  account  for  irrecoverable  viscous  flow.   The 
correspondence  principle  was  again  invoked  to  solve  the  problem,  and 
laboratory  tests  were  conducted  which  demonstrated  that  the  model  could 
provide  a  reasonable  approximation  of  the  properties  of  the  asphaltic 
concrete  material. 

Barksdale  and  Leonards  (reference  35)  investigated  long  term 
effects  due  to  repeated  loads  by  applying  the  correspondence  principle 
to  the  solution  of  a  multilayer  elastic  model.  A  series  of  Kelvin 
(spring  and  dashpot  in  parallel)  elements  with  both  short  term  and 
long  term  creep  compliance  characteristics  were  used  to  obtain  an 
operational  modulus  for  the  elastic  equations.   Predictions  of 
immediate  deflection  and  rutting  of  asphalt  pavements  at  the  AASHO 
Road  Test  compared  favorably  with  the  measured  data. 

Eynamic  Loading 

In  1962,  Pister  and  Westman  (reference  36)  introduced  a  quasi- 
static  approximation  of  a  moving  load  for  the  static  load,  elastic- 
viscoelastic  model  presented  earlier  by  Pister  and  Monismith  (reference 


33) .   Many  aspects  of  their  results  were  found  to  be  in  good 
qualitative  agreement  with  field  pavement  performance;  e.g.,  maximum 
deflections  varied  with  the  velocity  of  load,  viscoelastic  profiles 
were  not  symmetrical  as  in  the  elastic  case,  and  deflections  preceding 
the  load  were  found  to  be  less  than  those  that  followed. 

Harr  and  Lewis  (reference  37)  developed  a  theory  whereby  stresses 
and  deflections  could  be  calculated  for  a  series  of  rectangular  slabs 
lying  on  a  viscoelastic  foundation  and  subjected  to  a  moving  load. 
Inertia  forces  and  partial  support  consequent  to  linear  temperature 
variations  within  the  slab  were  also  taken  into  account.   This 
investigation  revealed  that,  contrary  to  prevailing  beliefs,  an 
increase  in  velocity  of  load  could  produce  an  increase  in  deflection 
and  stress. 

Viscoelastic  Solutions 

A  number  of  investigators  have  modeled  pavement  systems  as  being 
entirely  viscoelastic.   These  approaches  normally  consider  inertia 
forces  in  contrast  to  the  aforementioned  elastic  and  elastic- 
viscoelastic  solutions  (with  the  exception  of  Harr  and  Lewis  (reference 
37)).   The  foundation  for  these  studies  was  provided  by  Biot  (reference 
38)  in  1955  when  he  introduced  new  principles  for  the  formulation  of 
dynamic  problems  in  viscoelasticity  for  the  most  general  case  of 
anisotropy.   Biot's  approaches  were  based  on  a  variational  principle 
derived  from  irreversible  thermodynamics  and  the  correspondence 
principle. 
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Static  Loading 

Pister  (reference  39)  presented  the  solution  of  a  viscoelastic 
plate  on  a  viscoelastic  foundation,  wherein  the  effect  of  transverse 
shear  deformations  were  accounted  for  within  the  plate.   Viscoelastic 
material  properties  were  included  in  both  the  plate  and  the  foundation; 
however,  inertia  forces  were  neglected. 

Dynamic  Loading 

Harr  (reference  1+0 )  used  a  simple  Voigt  element  to  show  the 
influence  of  vehicle  speed  on  pavement  deflection.   The  analytical 
results  compared  favorably  with  AASHO  Road  Test  data  secured  at 
Ottawa,  Illinois  (reference  kl) .      Perloff  and  Moavenzadeh  (reference 
1+2)  presented  an  analysis  for  the  deflection  due  to  a  point  load 
moving  across  the  surface  of  a  semi-infinite,  homogeneous,  linear 
viscoelastic  medium.   A  Kelvin  model  was  used  to  characterize  the 
material;  inertia  forces  were  neglected.   Recently,  Szendrie  and 
Freeme  (reference  1+3),  in  conjunction  with  vibration  tests,  modeled  a 
pavement  system  using  seven  system  parameters  derived  from  a  series  of 
springs,  dashpots ,  and  layered  mass  representations. 

Viscoelastic  models  have  been  employed  recently  in  conjunction 
with  the  load- deflect  ion  stiffness  method  (references  1+1+ ,  1+5,  U6 ,  and 
1+7).   In  this  method  the  dynamic  input  is  in  the  form  of  a  steady-state 
vibratory  load,  and  the  ratio  of  the  applied  load  magnitude  to  the 
measured  deflection  is  termed  the  dynamic  stiffness  of  a  pavement 
structure. 
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Mechanistic  Solutions 

Each  successive  investigation  appears  to  add  insight  into  the 
action  of  pavement  systems;  however,  none  of  the  models  thus  far 
presented  are  able  to  represent  field  situations  with  sufficient 
predicative  capability  to  satisfy  either  existing  or  contemplated  needs. 
For  example,  in  the  classical  theory  of  elasticity,  the  requirements  of 
equilibrium  and  the  definitions  of  strain  (and  the  compatibility 
relations),  even  at  a  point,  provide  less  equations  than  the  unknowns 
contained  within  them.   Hence,  additional  assumptions  are  required  to 
provide  a  tractable  system.   Additional  independent  relations  are 
generally  introduced  at  a  point  by  inferring  certain  mechanistic 
properties  to  the  materials  —  commonly  referred  to  as  "constitutive 
relations".   In  elastic  theory  the  "constitutive  relations"  are 
assumed  to  relate  stresses  and  strains  in  a  linear,  homogeneous 
manner;  i.e.,  by  generalizing  Hooke's  law.   The  action  of  real  materials 
under  field  conditions  indicate  that  the  true  relationship  between 
stresses  and  strains  are  neither  linear  nor  homogeneous  and  depend 
greatly  on  time,  temperature,  and  previous  history. 

The  present  investigation  offers  an  alternate  approach,  which  is 
based  on  the  observed  input  and  output  responses  of  real  materials 
under  actual  field  conditions.   This  method  incorporates  a  far  greater 
degree  of  versatility  to  the  description  of  component  materials,  while 
at  the  same  time,  it  introduces  a  global  concept  to  the  response  of 
these  materials  to  loading  (in  lieu  of  the  point  concept  and  subsequent 
assumption  of  homogeneity  and  isotropy  required  in  the  case  of  the 
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classical  theories). 

The  approach  has  its  basis  in  three  earlier  investigations: 
Busching,  Goetz,  and  Harr  (reference  U8),  concerning  the  mechanistic 
description  of  a  deformable,  transversely  anisotropic  bituminous 
mixture;  Lai,  Goetz,  and  Harr  (reference  1+9),  obtaining  properties  of 
sheet-asphalt  mixtures  under  conditions  of  plain  strain;  and  Swami , 
Goetz,  and  Harr  (reference  50),  predicting  behavior  of  bituminous  mix 
under  a  variety  of  conditions  including  temperature  differentials. 
As  a  result  of  the  latter  study,  it  was  concluded  that  the  viscoelastic, 
or  time-dependent,  characteristics  of  an  asphalt  concrete  could  be 
represented  by  a  time  dependent  expression  called  a  "transfer 
function";  that  the  transfer  function  was  unique  for  the  sample  of  the 
material  at  a  given  temperature;  and  that  the  parameters  inherent  in 
the  transfer  function  were  better  indicators  of  the  material 
characteristics  than  previously  employed  models  (such  as,  Young's 
Modulus,  Poisson's  ratio,  etc.). 

Transfer  Function  Concept 

By  definition,  a  transfer  function  is  the  ratio  of  an  operational 
output  (expressed  as  the  Laplace  transform  of  the  output)  to  an 
operational  input  (Laplace  transform  of  the  input)  of  a  time-dependent 
system  (references  51,  52,  and  53).   If  a  system  is  subjected  to  a 
dynamic  input,  l(t),  and  evidences  a  resulting  output,  0(t),  then 
mathematically,  the  transfer  function  is  depicted  as 

G(s)=^M  (1) 

i(s) 
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where, 

G(s)  =  transfer  function, 

T(s)  =  Laplace  transform  of  input  function,  l(t),  and 

0(s)  =  Laplace  transform  of  output  function,  0(t). 

Note  that  the  transfer  function  is  the  ratio  of  two  Laplace  transforms, 

and  hence  is  a  function  of  the  complex  variable,  "s"  (references  5k 

and  55). 

"A  Priori"  Modeling 

Use  of  the  classical  (elasticity  and  viscoelasticity)  theories 
normally  use  "a  priori"  modeling;  i.e.  a  model  is  employed  to 
simulate  the  response  of  the  system  to  loading,  as  depicted  in  Figure  1. 
Inputs  of  interest  are  induced  in  the  model  after  material  and  system 
parameters  are  identified  and  predictions  are  made  of  the  behavior  of 
the  system.   Since  the  model  is  chosen  "a  priori",  the  number  of 
system  parameters  are  predetermined.   Inherent  in  this  approach  is  the 
existance  of  specifically  designed  laboratory  and/or  field  tests  to 
evaluate  the  necessary  parameters.   Comparisons  with  "real  world" 
behavior  indicates  the  degree  to  which  success  is  achieved. 

"A  Posteriori"  Modeling 

When  "a  posteriori"  modeling  inputs  are  set  forth,  consequent 
outputs  are  measured  and  the  two  are  related  in  a  ourely  mathematical 
way  (transfer  functions)  without  regard  to  simulating  anticipated 
effects.   Hence,  there  is  neither  the  preselection  of  system  parameters 
nor  a  modeling  in  the  "simulation  world"  (see  Figure  2).   In  contrast 
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to  "a  priori"  modeling,  the  technique  involves  the  desired  use  of  the 
real  world  system  itself. 

Determining  Transfer  Functions 

Transfer  functions  for  predicting  system  behavior  can  be  obtained 
by  either  "a  priori"  or  "a  posteriori"  modeling  techniques.   In  the 
case  of  "a  priori"  modeling,  generally,  a  mechanical  system  is 
hypothesized  and  a  mathematical  expression  is  obtained  relating  the 
transfer  of  input  to  output.   An  example  of  how  the  transfer  function 
is  obtained  using  "a  priori"  modeling  is  given  in  Appendix  A  (reference 
56). 

In  "a  posteriori"  modeling,  several  methods  -an  be  used  to 
determine  a  transfer  function.   One  method  is  to  approximate  the 
transfer  function  from  input  and  output  data  obtained  from  field  tests 
(viz.,  by  analyzing  the  frequency  spectrum  of  the  measured  inputs  and 
outputs)  as  proposed  by  Tsein  (reference  53)  and  Crafton  (reference  51 ) 
for  certain  types  of  input  and  output  data  from  laboratory  tests . 
Swami  (reference  57)  obtained  the  transfer  function  of  asphaltic 
concrete  specimens  in  the  absence  of  any  prior  knowledge  of  their 
dynamic  properties  by  using  input  and  output  data  in  the  manner 
indicated  above.   The  transfer  functions,  thus  obtained,  define  an 
"a  posteriori"  type  experimental  model.   A  second  way  to  determine  the 
transfer  function  for  a  system  is  to  find  the  output  of  the  system  due 
to  an  unit  impulse,  or  Dirac  Delta  function  input.   In  this  case  the 
measured  output  provides  directly  the  transfer  function.   A  third 
method  for  finding  a  system's  transfer  function  is  to  implicitly  invoke 
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the  process   of  real   convolution   (to   input   and  output   data)    in   the 
time   domain,    and  subsequently,    if  desired,   to  invert  the   function  to 
a  complex  domain.      Appendix  A  presents   some   details   of  each   of  the 
above  methods   for   finding   a  system's  transfer   function. 
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FORMULATION  OF  A  TESTING  PROGRAM 

Figure  3  represents  an  aircraft  moving  on  a  pavement  system. 
Figure  Ua  depicts  a  detail  of  this  figure  which  provides  a  series  of 
three  dimensional  plots  of  time-dependent  deflection  responses  (along 
a  line  normal  to  the  direction  of  aircraft  movement).   A  time-dependent 
deflection  response  function  represents  the  vertical  movement  at  a 
Doint  on  the  surface  of  pavement  system  as  a  "prime  mover"  passes  in 
the  vicinity  of  the  point.   The  deflection  resronse  may  be  a  measure 
of  the  surface  movement  in  terms  of  an  acceleration,  a  velocity,  or  a 
displacement.   In  the  present  work  only  displacement  responses  shall  be 
considered  in  further  developments . 

Formulation  of  Pavement  Transfer  Functions 

The  "signature"  of  a  prime  mover  is  the  vertical  response  at  some 
convenient  point  on  a  pavement's  surface  (defined  at  x  =  x  =  0) 
corresponding  to  the  input  of  a  prime  mover's  tires.   In  this  regard 
the  prime  mover  can  be  thought  of  as  providing  a  "quasi  input  force" 
equivalent  to  the  lumped  prime  mover  input  (e.g.,  truck  duals, 
aircraft  gear,  single  tire,  etc.).   In  subsequent  developments,  the 
signature  is  taken  to  be  the  time-dependent  input  fui  ction  to  the 
pavement  system. 

The  time-dependent  deflection  response  at  x  =  x  represents  a 
time-dependent  output  function  at  that  point  corresponding  to  the 
signature  input  function  at  x  =  x  .   Similarly  the  response  at  x  =  x, 
represents  the  time-dependent  output  function  at  the  point  x  =  xg  due 
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to  the  signature  at  x  =  x  ,  etc.   In  essence,  each  successive  time- 
dependent  deflection  function  is  defined  to  be  a  time-dependent  output 
function  of  the  signature  at  the  corresponding  points,  x  =  x  , 
x  =  x,...,  x  =  x.,  respectively. 

In  general,  following  equation  (l),  the  transfer  function  at  a 

point  on  the  pavement  at  x  =  x.  due  to  the  signature  at  x  =  x  can 

i  o 

be  written  as 

0x  (s) 

5   (s)  =  — i (2) 

ix  (s) 

o 


where, 


I   (s)  =  Laplace  transform  of  the  signature  at  x  =  x  (an  input 

J\.  o 

o 


response) , 


0   (s)  =  Laplace  transform  of  an  output  response  function  at 
i 

x  =  x.  ,  and 

G   (s)  =  the  transfer  function  between  the  operational  signature 
i 

at  x  =  x  and  an  operational  output  response  function 

at  x  =  x.  . 

i 

The  transfer  function  is  seen  here  to  be  expressed  as  deflection 

response  functions  at  two  different  points  due  to  a  quasi  input  force 

at  x  =  x  .   An  examination  of  equation  (2)  also  shows  that  to  be 
o 

meaningful  the  transfer  function  must  reflect  pertinent  material 

parameters  which  characterize  the  pavement  system  between  points 

x  =  x  and  x  =  x..   Hence,  in  philosophy,  the  transfer  function  provides 
o  1 

a  "body  concept"  to  the  pavement  in  contrast  to  the  "point  concept" 
inherent  in  "a  priori"  modeling  by  classical  (elasticity  and 
viscoelasticity)  theories. 
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In   summary,   the   signature   reflects  both  the   input  of  the  prime 
mover  and  the   consequent   feedback   of  the  pavement    system.      The  time- 
dependent  transfer   functions   are   seen   to  act   as    filters  which 
attentuate  the   available  energy  of  the   signature.      As   such,   they  must 
contain   component   parameters   of  the   pavement  which,    in   turn,    indicate 
the   performance  of  the   system. 

An   "A  Priori"   Model   Technique 

The   simple  theoretical   spring-dashpot  model   depicted  in   Figure    Ub 
will  be  used  to   further   illustrate  the   concept   of  a   transfer   function 
for  an   "a  priori"    system.      The   differential   equation  cf  motion   for  this 
particular  model  may  be  written   as    (references   53   and   57). 

0      (t)+£5      (t)+^0      (t)   =  M      (t)   +  £l      (t)  (3) 

x.  m     x.  m     x.  m     x  in      x 

1  1  l  o  o 

Taking  the  Laplace  transform  of  both  sides  of  equation  (3)  and 

rearranging, 

0      (s)  c  k 

x  —  s   ■+■  — 

s      i    \  i  m  "        m  ,  I  V 

G      (s)    =  =  -5 -  (M 

x.  =      /    v  £        c  k 

1  I      (s)  s     +  —  s   +  — 

x  mm 

o 


or, 


2 
-      /    \   _  2coas   +  co /    x 

Gx    (s)        ~2 2  (5} 

i  s "   +   2  was   +  u 


where, 


-{ 


and 


CO 

"  m 


c 

a  = 


cimu) 


23 


The   inverse  Laplace   transform  produces   the  time-dependent   transfer 
function   for  the   theoretical,    "a  priori"  model   as 

G      (t)    =        ;  jjj        •    e_aU)t    [2a  sin    (co\/l-a2t  +   0)   +   sin    Ul/l-a2t]      (6) 


and 


tan   0  =Y"^- 


2 

(7) 


A  characteristic  plot  of  G   (t)  versus  t  (for  a  <  l)  is  shown  in 

i 
Figure  5. 

"A  Posteriori"  Modeling  Techniques 

In  the  early  stages  of  this  research,  it  was  thought  that  a  low 
frequency  vibrator  could  excite  a  pavement  system  sufficiently  to 
provide  the  transfer  function  from  a  frequency  spectrum  analysis. 
To  investigate  this  possibility,  a  testing  program  was  undertaken 
at  the  Civil  Engineering  Research  Facility,  Kirtland  Air  Force  Base, 
New  Mexico,  using  a  low  frequency  vibrator.   The  vibrator  was  capable 
of  inducing  a  constant  force  between  50  and  1500  pounds  over  a 
frequency  range  of  approximately  15  to  50  cycles  per  second.   The 
output  responses  were  measured  with  Columbia  Model  1107-^ 
accelerometers.   A  typical  set  of  data  with  amplitude  of  surface 
motion  plotted  versus  frequency  of  input  force  is  depicted  in  Figure  6. 
Figure  7  represents  calculated  responses  for  a  theoretical  second  order 
viscoelastic  model.   A  comparison  of  the  two  plots  led  to  the  conclusion 
that  the  vibrator  used  in  the  tests  at  Kirtland  Air  Force  Base 
contained  frequencies  less  than  that  of  the  natural  frequency  of 
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Maximum 
Range  of  Data 


Data   Required  at    Least  to  This 
Point  ( i.e.,  Vibrator    Capable  of 
Greater   than    100  cps.) 


a   =   Driving    Frequency 
u)  =    System    Undamped    Nat. 
Frequency 


0.0 


.0  2.0 

Frequency   Ratio  ■  a/a> 


3.0 


40 


FIGURE    6.      FREQUENCY    SPECTRUM   PLOT  FOR   VIBRATOR   DATA 

F0  sin(at)    _ 


M 


Driving    Frequency 
Model    Undamped   Nat. 
Frequency 


1.0  2.0  3.0 

Frequency    Ratio  =  a/w 


4.0 


FIGURE    7.     FREQUENCY   SPECTRUM    PLOT  FOR   SIMPLE    MODEL 
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the  pavement  system  tested.   To  determine  the  transfer  functions  for 
a  system  from  its  frequency  spectrum,  it  is  necessary  to  excite  the 
system  to  frequencies  three  to  four  times  that  of  its  natural  frequency; 
that  is,  it  is  necessary  to  reach  magnitudes  of  frequency  which  display 
a  constant  slope  at  the  upper  bound  of  the  frequency  spectrum  curve 
for  the  system.   This  observation  coupled  with  the  lead  time  and  cost 
required  to  design  and  fabricate  an  adequate  vibrator  negated  the 
frequency  spectrum  technique. 

Attention  was  then  directed  to  the  "a  posteriori"  impulse  input 
method  (see  Appendix  A).   Unfortunately,  this  method  also  necessitated 
new  and  expensive  equipment;  hence,  it  was  also  deemed  unfeasible. 

Finally,  the  "a  posteriori"  real  convolution  technique  (see 
Appendix  A)  was  explored.   This  technique  provides  a  time  domain 
solution  which  permits  the  deflection  responses  to  be  used  directly 
in  calculating  a  time  domain  transfer  function.   The  method  also 
eliminates  the  need  to  obtain  Laplace  transforms  and  inverse 
transforms. 

Beginning  with  the  definition  of  a  transfer  function  for  the 
physical  situation  as  given  in  equation  (2),  it  follows  that 


0x  (s)  =  Gx  (s)  •  Ix  (s)  (8) 

i        i        o 

As  equation  (8)  is  seen  to  represent  the  product  of  two  Laplace 

transforms,  0   (t)  is  said  to  be  the  "real  convolution  of  G   (t)  and 
Xi  Xi 

I   (t)".   Invoking  the  real  convolution  theorem  (references  5U,  58, 
x 
o 


and  59)  equation  (8)  can  be  expressed  as 
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(9) 


An  implicit  solution  of  equation  (9)  was  devised  to  yield  G   (t), 

x. 

l 

termed  the  inverse  transform  of  the  transfer  function,  or  the  time 
domain  transfer  (TDT)  function,  hereinafter  referred  to  as  a  TDT 
function.   Equation  (9)  was  discretized  to  facilitat  numerical 
solutions  by  a  computer.   The  integral  was  solved  implicitly  for 
G   (t),  provided  0   (t)  and  I   (t)  were  available,  say,  as  arrays  of 

X .  X .  X 

1  1  o 

numerical  values.   Conversely,  an  explicit  solution  was  obtained  for 

0   (t),  when  G  (t)  and  I   (t)  were  known.   Recursion  formulas  used  to 

i  i  o 

develop  computer  subroutines  for  the  implicit  and  explicit  solutions 

of  the  real  convolution  integral  are  provided  as  equations  (10)  and 

(11). 

Implicit  convolution — 

k-1 


Gx.^=-~ J=1     I1    (1)    •    At 


°x  <k>  -5_GX  <j)  •  ix  u-j+D  •  At 

i j=l       i    o_ 

"1      (1)    •    At 
i  xx. 

o  1 

Explicit  convolution — 

k 
°v  (k)  =51  Gv  (J}  *  Zv  (k-J+1)  *  AtY  (11) 

X.        *■  -    X.  X  X 

1        j=l    1  O  1 

Accelerations  and/or  velocities  could  be  used  to  determine  the  TDT 
functions,  although  displacements  have  been  indicated  here.   For 
example,  if  the  output  responses  are  given  as  the  accelerations  of  a 
pavement  systems  surface, 
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s20   (s)   0   (s) 

X  X 

G  (s)=_ 1 .-J (12) 

Xi      s%  (s)   I  (s) 
X         x 
o        O 


Signature  and  TDT  Functions 

Equation  (10)  establishes  the  TDT  function  once  a  pavement's 
response  functions  are  known  at  two  points.   To  obtain  the  requisite 
response  functions,  the  time  history  of  the  vertical  motion  of  a 
pavement's  surface  is  required  at  specific  locations  adjacent  to  the 
vehicle's  path.   The  actual  location  of  a  prime  mover's  signature  is 
arbitary;  however,  once  selected,  its  position  must  be  fixed  relative 
to  the  gages  used  to  measure  the  deflection  responses.  This  is 
necessary  because  a  TDT  function  serves  as  a  filtering  mechanism  to 
transform  a  signature  at  one  point  into  a  deflection  response  at  an 
adjacent  point;  therefore,  the  distance  over  which  a  TDT  function  is 
calculated  must  remain  constant. 

It  was  not  considered  feasible  to  confine  the  path  of  a  prime 
mover  to  a  precise  location  relative  to  the  installed  gages;  so,  a 
method  was  developed  to  accomodate  the  shifting  of  the  response  function 
measurements.   The  technique  employed  to  accomplish  the  minor  spatial 
adjustments  for  the  deflection  response  functions  is  predicated  on  the 
assumption  that  the  peak  dynamic  deflection  basin  attentuates  according 
to  a  normal  distribution  curve.  The  spatial  adjustments  were  considered 
to  be  small  compared  to  the  extent  of  the  dynamic  deflection  basin, 
so  that  any  error  in  the  modeling  process  would  thereby  be  minimized. 
This  was  later  found  to  be  the  case.   The  correction  necessary  to  shift 
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a  response   function  must,    in   essorce,    scale   the   response   function 
spatially;   yet,    it  must   also  retain   its   shape.      This  was    accomplished 
by   forming  a  scaling  factor  defined   as 


x.+A,x.    =   V+A(Deal0/0x>eak)  (13) 


c 
Ji    ■■-!  ~i    "  -1 


where, 

A   =   Offset   distance   of   innut    (i.e.,    distance  between   the    actual 

ooint  where   a  prime   mover's    signature    is    applied    and  a 

defined  signature   input   point   on  the   pavement    relative   to 

the   gages ) , 

0    .(peak)  =  The  peak  value  of  the  deflection  response  function 
1 

at  a  point  located  at  x  =  x.+A  from  the  defined 

signature  input  point, 

i  =  0,  1,  2,...,  i  (i  =  0  represents  the  signature), 

and, 

n      ,        ,  x    .     -B(.0Ul7  x. )2  ,.,  > 

0   (peak)  =  A  •  e         l  (In) 

i 

or, 

n    (        n  8  a  .   -3(.ol4l7(x. +  >))''  (15) 

0   ,  '  peak  -  A   e  l 

x.+A 

l 

The  A   and  B   parameters    associated  with   the  ncrma_L   distribution  model 
are    obtained   by    a  regression    analysis    or.    neak   defle^ti'>n   values.      Figure 
8  shows   a  2-dimensional  representation  of  how  the   method   is   employed. 
For  example,    consider  that   a  prime  mover  passes   the  gag;es   at   a  distance 
x  =   A  to  the   left   of  the  first   gage.      The  peak  dynamic   deflections   are 
fitted  to  a  normal   distribution  model;    the   factors, 


30 


"«/> 

■o 

01 

*■ 

<*- 

!£ 

in 

i_ 

o 

.  *-            «fl 

«/>«/>            £ 

c  c            .9 

• 

.2  .2          ■  S 

</> 

lect 
lect 

ions 
efle 
oad 

z 
o 

0)  a>         n  Q 

1- 

u 

c 
o 

4— 

a 

/namic  D 
ynamic  D 
ocations. 
re  Appiio 
ignature 
of  Input 
pacing  . 

o 

z 

D 
U. 

UJ 

IT) 

1                  QQJjtf)         ifl 

</> 

eak 

eak 

age 

igna 

eak 

ffse 

age 

z 
o 

Q_ 

1 

V   " 

Q_  Q.  O  (/>  Q-  O  O 

UJ 

X 

\l     * 

i 

1      •  =,*  <  "<    T3 

QC 

-  \ 

1         = 

O 

f  ' 

1 

\\ 

U. 

1 

„\\ 

K 

X 

\\ 
\\ 

Z 
Ui 

f ' 

1 

■o  «~»  \  \ 

D 

in 

CM 

J£       \    \ 

-> 

O 

X 

o      \ 

■*  \\ 

Q 

< 

r 

13 

1 

T 

1  \\ 

3       »T"          \      \ 

00 

ro 

ox       \\ 

O 

X 

1       \  ^ 

</> 

UI 

a: 

CM 

•*i 

.X   1, 

1      \\ 

o                          \    \ 

CO                                    \    > 

x                                         \ 

o                         \ 

H- 

in 

± 
1  ° 

1  * 

c 
o 

o 
a> 

Q 

u. 

r 
J, 

">    h < H 

31 


c  c 

x  +A,x  .....   x.+X,x. ,  are  computed;  and  the  correction  process  is 
o    o        1    1 

employed. 

The  testing  program  had  as  an  objective  the  identification  of 
the  individual  tests  necessary  to  validate  the  hypothesis  for  this 
investigation.   The  hypothesis,  as  initially  set  forth,  stated  that 
TDT  functions  could  be  calculated  from  a  pavement's  input  and 
output  deflection  responses;  and  that  these  TDT  functions  could  be 
obtained,  in  a  mathematical  sense,  without  the  need  to  simulate 
respective  material  performance  to  determine  values  for  preselected 
descriptors. 

In  summary,  the  deflection  response  functions  at  x  =  x. 

(i  =  0,  1,  2,...,i)  relative  to  a  signature  applied  at  x  =  x  ,  as 

transferred  from  x  =  x  +X,  are  given  when  the  measured  deflection 

o 

responses  are  multiplied  by  the  factors  defined  from  equation  (13). 
Then,  equation  (10)  is  invoked  to  determine  the  pavement's  TDT 
functions.   Once  deflection  response  functions  are  measured,  the 
pavement's  respective  prime  mover's  signature  and  TDT  functions  may 
be  determined  as  follows : 

(1)  Fit  peak  dynamic  deflections  to  a  normal  distribution  model  - 
least  squares  fit  of  equation  (lU). 

(2)  Determine  signature  and  deflection  responses  with  respect  to 
a  common  origin  -  apply  correction  factor  defined  by  equation  (13). 

(3)  Calculate  TDT  functions  -  apply  technique  defined  by 
equation  (10). 

A  listing  of  computer  code  NMLFIT  for  accomplishing  item  (l)  is 
provided  in  Appendix  B.   Appendix  C  provides  a  listing  of  computer 
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code  CONVLI  which  contains  maneuvers  for  items  (2)  and  (3). 

Predicting  Performance 

As  a  secondary  feature,  and  dependent  on  being  able  to  calculate 
a  pavement's  TDT  functions,  it  was  further  hypothesized  that  these 
functions  could  be  employed  to  predict  output  deflection  response 
functions  when  the  pavement  was  subject  to  any  imposed  load.   In  the 
previous  section,  the  procedures  for  ascertaining  a  prime  mover's 
signature  and  a  pavement's  TDT  functions  were  outlined.   Once 
sufficient  catalogues  of  the  signatures  and  TDT  functions  of  different 
prime  movers  and  pavement  systems  are  obtained,  estimates  can  be  had 
of  combinations  of  vehicles  and  pavements  for  which  tests  were  not 
conducted.   For  example,  with  a  cataloged  signature  and,  say,  TDT 
functions  obtained  from  a  particular  vehicle,  a  third  computer  code 
CONVLE  (see  Appendix  D  for  a  listing)  may  be  invoked  to  calculate 
deflection  response  functions  at  the  x.  points.   Figure  9  summarizes 
a  complete  procedure  for  predicting  a  pavement's  performance  when 
subjected  to  a  particular  prime  mover.   The  procedure  outlines  the 
combinations  of  catalogued  signatures  and  TDT  functions  necessary  to 
perform  these  predictions. 

Referring  to  Figure  9,  two  arbitary  pavements  are  shown: 
pavement  A  on  the  left  is  the  pavement  for  which  signatures  of  both 
prime  movers  #1  and  #2  have  been  previously  cataloged  (the  pavement's 
TDT  functions  are  also  known) ;  and  pavement  B  on  the  right  is  the 
pavement  whose  deflection  response  functions  for  prime  mover  #2  are 
desired.   Since  the  deflection  response  functions  are  desired  for 
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Pavement   A 


—  f 


Deflection  Responses  Deflection  Responses 

( Prime  Mew  •  I  )  (  Prime  Mover  »  2  ) 


CONVLI    Cod* 


Pavement    8 


Deflection  Responses  Deflection  Responses 

(Prime   Mover  »  I  )  (  Prime  Mover  «  2  ) 


CONVLI   Code 


Signature  ,  IB 

L 


Account  for  Pavement  Feedback  to  Signature 
Gl. 


CONVLE  Code 


Predicted    Deflection 
Responses  (Prime  Mover"  2) 


FIGURE    9-    PERFORMANCE    EVALUATION    METHOD 
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prime  mover  #2,   prime  mover  #1  shall  be  termed  a  "standard"   vehicle 

to  illustrate  the   prediction  method.      As   a  first   step,    consider  the 

standard  vehicle    (prime  mover  #l)    on  pavement  B.      Deflection  response 

functions   are  measured,    and  the   vehicle's   signature   and  the  pavements 

TDT   functions    are   calculated  using   computer   code   CONVLI.      The  next 

step   accounts    for  the  pavement's   feedback   to  the  signature.      In   the 

early  part   of  this   investigation,    it  was   assumed  that  the   deflection 

magnitudes   of  a  particular  prime  mover's   signatures   on  different 

pavements   could  be  treated  linearly.      This  was   later  found  to  be 

the   case.      Therefore,   by  using  the  three  known  signatures,   it   follows 

that  the   fourth  signature   for  prime  mover  #2  on  pavement  B  may  be 

expressed  as 

I    (peak) 

Jx   (t)  =  /^(peak)  •  Tx   (t)  (l6) 

Xo2B      1x     /Peak;    X02A 
olA 

Lastly,  the  fourth  signature,  equation  (l6),  and  the  pavement  B  TDT 

functions,  from  the  first  step,  are  used  in  computer  code  CONVLE  to 

calculate  the  corresponding  deflection  response  functions  for  prime 

mover  #2  on  pavement  B.   Thus,  the  desired  performance  predictions 

are  achieved. 

The  deflection  response  functions  provide  the  true  deflection 

basin  for  any  instant  of  time  desired.   In  addition,  the  deflection 

basin  is  not  predicated  on  use  of  predetermined  parameters  as  is 

done  in  the  "a  priori"  modeling  approaches.   Inherent  in  the  TDT 

functions  are  the  numerous  pavement  "properties";  yet,  none  have  to 

be  evaluated  directly  in  order  to  assess  the  pavement's  performance. 
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FIELD  INVESTIGATION 

The  field  investigation  phase  of  this  study  had  as  its  objective 
the  development  and  use  of  techniques  for  obtaining  all  the  data 
required  to  determine  a  pavement's  TDT  functions  by  the  convolution 
method.   The  main  feature  of  the  technique  was  to  devise  a  system 
which  could  accurately  measure  the  time-dependent  deflection  response 
functions  of  a  pavement  when  subjected  to  a  prime  mover's  dynamic  load 
input . 

Exploratory  field  tests  were  conducted  to  ascertain  the  criteria 
needed  to  design  the  measuring  system.   These  tests  were  necessary 
as  the  deflection  response  functions  were  transient  in  nature,  rather 
than  steady  state,  and  neither  the  magnitude  of  deflections  nor  their 
time  durations  were  known.   Once  the  information  from  these  tests  had 
been  analysed,  it  was  possible  to  design  an  adequate  measuring  system. 
After  fabrication  and  installation  of  the  measuring  components,  a  full 
scale  field  test  program  was  carried  out  to  obtain  the  data  required 
to  validate  the  performance  evaluation  method  summarized  by  Figure  9. 

Exploratory  Tests 

Field  exploratory  tests  were  conducted  on  an  operational  taxiway 
at  Kirtland  AFB,  New  Mexico  (Albuquerque  International  Airport).   They 
had  as  their  primary  objective  acquisition  of  data  necessary  to  select 
a  suitable  energy  sensing  device,  and  to  determine  locations  on  the 
pavement  adjacent  to  a  prime  mover's  path  which  could  adequately 
define  the  dynamic  deflection  basin.   Both  vertical  deflections  and 
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accelerations  of  the  pavement's  surface  were  considered,  and  a  series 
of  tests  were  conducted  using  a  specially  built  load  cart  and  a 
Boeing  727  aircraft  (see  Appendix  E).   Appendix  E  also  provides  details 
concerning  installation  and  calibration  of  the  sensing  devices 
employed  as  well  as  typical  plots  of  deflection  response  functions 
obtained.   Table  1  gives  a  comparison  of  the  test  results  for  each 
prime  mover  with  respect  to  the  particular  sensing  device  used. 
Conclusions  from  the  data  acquired  in  these  tests  indicated  that  a 
deflection  gage  employing  a  linear  variable  differential  transformer 
(LVDT)  was  the  device  best  suited  for  measuring  the  time-dependent 
functions. 

The  peak  magnitudes  of  the  deflection  response  functions  and  their 
relative  attentuation  with  distance  from  the  prime  mover's  path  guided 
the  design  of  the  deflection  gages;  e.g.,  they  established  the 
required  LVDT  sensitivity  and  the  gage  spacing,  respectively. 

Field  Tests 


The  prototype  field  tests  were  also  conducted  at  Kirtland,  AFB, 
N.M.  (Albuquerque  International  Airport).   In  contrast  to  the 
exploratory  tests,  however,  this  testing  program  was  designed  to 
control  certain  variables  which  were  thought  to  influence  a  pavement's 
TDT  functions;  e.g.,  data  was  only  secured  when  pavement  temperatures 
were  between  55  F  and  65  F,  and  the  prime  mover  speeds  were  controlled 
to  approximately  five  miles  per  hour  for  each  data  run. 
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Instrumentation 

Time-dependent  deflection  response  functions  were  measured  with 
specially  designed  deflection  gages  using  LVDT  displacement  transducers. 
The  overall  system  was  designed  to  monitor  a  pavement's  vertical 
movement  at  six  different  locations  simultaneously;  six  locations  on  a 
line  normal  to  a  prime  mover's  path.   Also,  since  it  was  necessary 
to  direct  the  prime  movers  close  to  some  of  the  sensing  devices, 
special  consideration  had  to  be  given  to  protect  these  gages  from  a 
prime  mover's  tire;  these  special  gages  were  placed  within  the 
pavement's  wearing  course,  and  were  labelled  type  I  installations. 
When  provisions  were  not  required  to  protect  the  gages,  a  beam  was 
used  to  serve  as  the  LVDT  mounts.   The  regular  deflection  gages  were 
labelled  type  II  installations. 

Following  from  the  combination  of  theory  and  experience  associated 
with  the  exploratory  tests,  deflection  gage  locations  relative  to  the 
path  of  a  prime  mover  were  selected  as  shown  in  Figure  10  (both  type  I 
and  type  II  installations  are  depicted).   Appendix  F  describes  each 
type  gage,  installation  details,  and  the  method  used  to  calibrate 
the  deflection  gages. 

The  following  provisions  were  incorporated  into  the  overall 
data  acquisition  technique: 

(l)   The  gage  outputs  were  recorded  on  magnetic  tape  with  an 
Ampex  FR  1300  tape  recorder  and,  in  order  to  ascertain  the  validity 
of  the  data,  selected  gages  were  monitored  on  an  oscilloscope  at  the 
time  of  a  test  run. 
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(2)  An  IRIG  time  generator  ( Culton  Industries)  was  used  to  record 
a  time  scale  on  the  magnetic  tape.   This  record  was  later  used  to 
select  time  intervals  for  digitizing  and  analyzing  the  data. 

(3)  The  temperature  of  the  asphaltic  concrete  wearing  course 
was  monitored  with  a  thermometer  embedded  approximately  one  inch  into 
the  pavement's  surface. 

(h)      A  DC  triggering  circuit,  actuated  by  pressure  switches 
located  2k   feet  on  each  side  of  the  deflection  gages,  was  installed 
for  checking  the  velocity  of  the  prime  mover  for  each  test  run.   The 
DC  voltage  was  also  recorded  on  the  magnetic  tape. 

(5)   Adhesive  tape  was  applied  to  the  pavement's  surface 
parallel  to  the  line  of  deflection  gages  to  facilitate  a  physical 
measurement  of  the  prime  mover  tire  relative  to  the  first  deflection 
gage.   The  distance  between  the  location  of  the  first  deflection  gage 
and  the  centroid  of  the  prime  mover  tire's  imprint  on  the  tape 
provided  the  offset  distance. 

Figures  11  through  lU  represent  the  type  I  and  type  II 
installations,  respectively.   Figure  15  shows  the  overall  instrumenta- 
tion setup,  and  Figure  l6  is  a  view  of  the  recording  equipment 
components. 

Scope  of  the  Field  Tests 

The  prototype  investigations  were  conducted  at  three  sites  (see 
Figure  IT)  with  pavement  cross-sections  as  shown  in  Figure  18.  Four 
different  prime  movers  were  used  at  each  test  site;  C135»  C131  and  C130 
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FIGURE    II.       TYPE  I  INSTALLATIONS,  DISTANT 


FIGURE    12.     TYPE  I  INSTALLATIONS,  CLOSE-UP 
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FIGURE    13.      TYPE  H  INSTALLATIONS,  DISTANT 


FIGURE    14.     TYPE  H  INSTALLATIONS,  CLOSE-UP 
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FIGURE     15.      INSTRUMENTATION   SETUP 


FIGURE    16.     RECORDING    EQUIPMENT  COMPONENTS 
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Test 
Site    No. 

Taxi  way 
No. 

®  Asphalt 
Concrete 

®  Base 
Course 

(c)     Subgrade 

1 

6 

3" 

9" 

Silt/ Sand 

2 

2 

6.5" 

8" 

Silt /Sand 

3 

8 

10" 

6" 

Highly     Compacted 
Silt/Sand 

FIGURE     18.     PAVEMENT    CROSS -SECTIONS 
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type  aircraft,  and  a  load  cart.   Figures  19  through  22  depict  each 
prime  mover,  respectively.   A  minimum  of  four  forward  passes  of  each 
prime  mover  was  made  at  each  site;  the  number  of  passes  of  the  prime 
movers  was  limited  by  two  factors  —  cost  and  availability  of  the 
prime  movers,  and  the  number  of  passes  considered  necessary  to  provide 
a  good  definition  of  the  dynamic  deflection  basin.   Table  2  summarizes 
the  physical  characteristics  of  each  prime  mover.   The  main  wheel 
and/or  gear  configurations  are  illustrated  in  Figure  23. 
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FIGURE    19.      PRIME   MOVER  -  C  135   AIRCRAFT 
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FIGURE    20.     PRIME    MOVER  -  C    131    AIRCRAFT 
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FIGURE    21.      PRIME  MOVER  -  C  130  AIRCRAFT 


FIGURE    22.     PRIME   MOVER  -  LOAD    CART 
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Single  Wheel 
(Load    Cart) 


h—   26"—) 
Twin 
(  CI3I  Aircraft  ) 


h—  36" — h 
Twin -Tandem 
(CI35  Aircraft) 


60* 


60" 


Single -Twin 
(CI30  Aircraft) 


FIGURE  23.   PRIME  MOVER  MAIN   GEAR  CHARACTERISTICS. 
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RESULTS  OF  THE  FIELD  TESTS 

Peak  Dynamic  Deflection  Basin 

Peak  deflection  response  values  of  the  dynamic  tests  for  all  the 

prime  movers  at  the  three  test  sites  are  summarized  in  Tables  3 

through  6.   These  values  were  then  fitted  to  a  normal  distribution 

2 
model  of  the  form  0   (peak)  =  A  •  e   *     i   by  conducting  a  least 

Xi 
squares  regression  analysis.   The  resulting  parameter  values  and 

correlation  coefficients  are  provided  in  Table  7.   The  A  factors  were 

adjusted  linearly  to  an  arbitrary  standard  prime  mover  weight  because 

the  weights  of  the  aircraft  varied  from  site  to  site  depending  on 

the  quantity  of  fuel  on-board  at  the  time  the  tests  were  conducted. 

Finally,  a  ratio  factor  defined  as  the  ratio  of  the  A  factor  divided 

by  the  largest  A  factor  at  each  test  site  was  calculated  for  each  prime 

mover,  respectively.  The  adjusted  A  factors  and  the  ratio  factors  as 

defined  above  are  also  provided  in  Table  7. 

The  relationship  between  A  factors  and  the  ratio  factors  are 

shown  in  Figure  2U;  and  the  range  of  the  B  parameters  are  provided 

in  Figure  25. 

Deflection  Response  Functions 

In  each  of  the  dynamic  tests,  deflections  were  recorded  over  a 
one  to  two  minute  period  starting  from  the  time  a  prime  mover  began  to 
accelerate  until  it  had  passed  beyond  the  vicinity  of  the  deflection 
gages.   Oscillographs  of  the  recorded  data  reflected  measurable 
deflections  over  only  an  approximate  two  second  period,  or  whenever 
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the  prime  mover  passed  within  ik   to  18  feet  of  the  gage  locations. 
The  IRIG  time  signal  was  used  to  select  a  9.5  second  interval  over 
which  to  convert  the  gage  output  voltages  to  digital  units.   Prior  to 
ditigitizing,  however,  the  data  were  filtered  electronically  to 
remove  any  interference  having  frequencies  greater  than  100  Hz.   A 
sampling  rate  of  Uo  units /second  was  used  in  the  ditigitizing  process, 
thus,  providing  380  discrete  values  for  each  deflection  response 
function. 

To  achieve  a  common  null  point  for  all  the  data,  the  first  ten 
discrete  values  of  each  deflection  response  function  was  averaged  and, 
in  turn,  this  average  value  was  substracted  from  each  discrete  value 
in  the  array.   This  procedure  was  necessary  since  it  was  observed 
that  the  null  for  the  respective  gages  were  different.   Some  of  the 
null  values  also  changed  after  each  test  run  because  of  the  pavement's 
viscoelastic  behavior.   The  null  for  the  gages  close  to  the  point  of 
signature  application  changed  most,  whereas,  it  changed  very  little 
for  gages  at  adjacent  locations.   In  the  extreme  case,  it  was  noted 
that  the  null  could  change  by  as  much  as  ten  percent  of  the  peak 
dynamic  deflection  for  a  gage  located  at  a  point  where  the  signature 
was  applied. 

The  data  contained  in  the  deflection  response  function  arrays  were 
truncated  to  reduce  the  computer  time  required  to  calculate  TDT 
functions  and  to  predict  the  desired  deflection  functions.   This 
procedure  follows  from  using  the  9-5  second  interval  over  which  to 
ditigitize  the  recorded  data.   An  approximate  two  second  interval 
would  have  been  large  enough  had  sufficient  experience  been  obtained 


6o 


vith  the  ditigitizing  process.  Therefore,  it  was  decided  to  truncate 
the  deflection  response  function  arrays  containing  380  discrete 
values  to  retain  only  the  significant  data.   The  truncation  was 
accomplished  by  deleting  all  initial  values  in  a  signature  array  which 
were  smaller  than  .002  inch.   Subsequently,  the  number  of  discrete 
values  between  the  first  array  value  that  was  retained  and  the  peak 
array  value  was  counted.   The  length  of  the  array  was  then 
established  as  three  times  this  counted  number.  After  the  length  of 
the  array  was  established  in  this  manner,  the  latter  array  values 
(which  were  of  no  particular  significance)  were  deleted.   The  interval 
of  the  signature  array  containing  only  the  significant  values  was 
taken  as  the  common  interval  over  which  to  retain  the  discrete  values 
for  successive  deflection  response  functions  that  were  adjacent  to  the 
point  where  the  signature  was  defined. 

Linear  interpolation  between  successive  discrete  values  of  the 
response  functions  was  then  employed  so  that  a  common  time  base  of 
1.8U  seconds  containing  UT  discrete  values  could  be  achieved.  The 
reason  for  providing  a  common  time  base  was  twofold:  (l)  to 
standardize  the  number  of  discrete  values  for  the  calculated  TDT 
functions;  and  (2)  to  account  for  fluctuations  in  prime  mover 
velocities. 

Lastly,  the  deflection  response  function  arrays  were  shifted 
with  respect  to  their  respective  offset  distance  as  defined  in 
Figure  8.   The  scaling  factors  given  by  equation  (13)  were  multiplied 
by  the  deflection  array  values  for  each  set  of  test  data.  This 
shift  would  not  have  been  necessary  had  the  deflection  functions  been 
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measured  relative  to  the  prime  mover  instead  of  at  fixed  locations 
on  the  pavement. 

Computer  code  CONVLI  (Appendix  C)  has  been  provided  with  all 
the  instructions  required  to  accomplish  the  procedures  outlined  above. 
Figures  26  through  3^  show  plots  of  the  deflection  response  functions 
after  these  procedures  were  applied.   It  is  to  be  noted  that  the 
ordinate  scales  are  variable  in  these  plots,  and  because  a  number  of 
data  runs  were  available  for  each  prime  mover  at  each  test  site,  the 
plots  represent  an  average  of  the  measured  data.  Appendix  G  contains 
computer  sheets  listing  the  discretized  data,  printed  out  to  six 
decimal  places. 

TDT  Functions 


The  denominator  in  equation  (10)  is  always  equal  to  the  product 

of  the  first  discrete  signature  array  value  times  the  discretized 

time  interval;  i.e.,  I  (l)  •  At.   It  was  found  necessary  to 

o 
substitute  a  value  larger  in  magnitude  for  I  (l)  in  order  to  obtain 

o 
convergence  for  the  implicit  convolution  process.  This  substitution 

was  accomplished  with  the  understanding  that  a  reciprocal  value  would 

have  to  be  used  to  obtain  the  suitable  inverse  when  using  the  explicit 

convolution  technique  defined  by  equation  (ll).  A  value  of  ten  times 

the  peak  signature  deflection  was  chosen  to  replace  the  initial 

signature  array  value.   In  essence,  this  substitution  provided  a  value 

for  I   (l)  which  was  near  unity,  or  practically  500  times  it  true  value, 
o 
The  procedure  outlined  above  was  initiated  because  in  the  early 

stages  of  this  analysis,  the  discretized  TDT  functional  values  were 
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FIGURE    26.    DEFLECTION  RESPONSE   FUNCTIONS 
( DATA   FOR    C  135    SITE    I  ) 
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FIGURE     27.        DEFLECTION  RESPONSE  FUNCTIONS 
(DATA  FOR   C  131  SITE   I  ) 


6h 


.10 

[x 

r 

x  =  0   (Signature) 

05 

00 

r         i 

_i...                1                    .                    1                                         1- - 

i 

.4 


°x 

r       I 


.02 


°x 

r     2 


.0008 


.0004 


.0000 


.8 


x=  18* 


1.2 


x=36' 


j  t 
2.0 


_i  t 

2.0 


J  t 


2.0 


=b 


.8 


J  t 


1.2 

x  =  72" 

1 


2.0 


.8 


1.2 


1.6 


2.0 


FIGURE    28.   DEFLECTION    RESPONSE   FUNCTIONS 
(  DATA   FOR  C  130 -SITE   I  ) 


65 


x  =  0  (  Signature  ) 


_i  t 
2.0 


1.6 


_i  t 
2.0 


.008 


.004 


,000 


x*36' 


2.0 


.8 


1.2 


J  t 


2.0 


J  t 


2.0 


FIGURE    29.       DEFLECTION    RESPONE   FUNCTIONS 
(DATA  FOR  C  135  SITE  2) 


66 


.002 


x  »  0  (  Signature  ) 


_i  t 
2.0 


_it 
2.0 


_J   t 
2.0 


.8 


2.0 


L*4 


x=72" 


.8 


2.0 


FIGURE    30.       DEFLECTION  RESPONSE  FUNCTIONS 
( OATA  FOR  C  131  SITE  2  ) 


67 


x  =0  (  Signature  ) 


j  f 


2.0 


02 

rX| 

^-^          x  =  18" 

.01 

00 

r      .       i 

^ —  i.i. 

.002 


.001  - 


.000 


.0003  - 


.8 


'x4 


.8 


x-36' 


_L 


1.2 
x  =  54" 


j  t 


2.0 


1.2 
x  =  72" 


1.6 


±: 


J    t 


2.0 


_J  t 
2.0 


.8 


1.6 


2.0 


FIGURE    31.     DEFLECTION  RESPONSE  FUNCTIONS 
(  DATA  FOR  C  130  SITE  2  ) 


68 


.025 


x-0  (Signature  ) 


x=l8' 


1.6 


j  r 


2.0 


J  t 


2.0 


.04 


x  =  36' 


.8 


ju. . i  t 


2.0 


°x 


x  =  72 


I         i        =fc 


_i  t 

2.0 


.4 


.8 


1.6 


FIGURE   32.    DEFLECTION  RESPONSE  FUNCTIONS 
( DATA  FOR  C  135  SITE  3  ) 


69 


.02  - 


.005  - 


.000 


x  =  0    (Signature) 


_j  t 
2.0 


i         l 


x-36' 


_i  t 
2.0 


J  t 


2.0 


j  r 


2.0 


r 


x  =  72' 


.8 


1.2 


1.6 


j  t 


2.0 


FIGURE   33.      0EFLECTI0N   RESPONSE  FUNCTIONS 
(DATA   FOR  C  131   SITE  3) 


70 


04 

[x 

r 

x  =  0  (Sign 

02 

00 

— ^\       i 

i.i. 

.4 


.8 


J  f 


2.0 


J i L 


.8 


2.0 


.02 

.01 
.00 


.000 


0 


'x4 


i 

i 

m i  f 

.4 

.8 

1.2 

1.6 

2.0 

Ox 
r 

x  =  54" 

■ 

•         " 

•        i 

>         1 

i 

— -        ' 

.8 


1.2 
x  =  72" 


3: 


.4 


.8 


J  t 


2.0 


j  t 


1.6 


2.0 


FIGURE    34.    DEFLECTION  RESPONSE  FUNCTIONS 
(DATA  FOR  C  130  SITE  3  ) 


71 


found  to  increase  in  magnitude,  almost  without  bound,  as  time  increased 
(or  as  consecutive  array  values  were  computed).   Close  examination  of 
the  implicit  convolution  technique  revealed  that  the  constant 
denominator  in  equation  (10)  was  of  the  order  of  magnitude  of  10   ; 
and  since  the  summation  process  was  repeated  U6  times,  there  was 
reason  to  believe  that  repeated  division  by  such  a  small  number  could 
cause  this  difficulty.  It  was  also  noted  that  each  consecutive  step 
carried  over  inflation  of  its  result  to  the  next  step  in  the  summation 
process.  To  examine  the  effect  of  this  division  in  the  procedure, 
discretized  analytical  input  and  output  functions  were  used  to 
simulate  deflection  response  functions.  The  implicit  convolution 
procedure  was  used  to  calculate  a  discretized  TDT  function.  Then, 
explicit  convolution  between  this  calculated  TDT  function  and  the 
original  analtyical  input  function  was  used  in  an  effort  to  reproduce 
the  original  analytical  output  function.  Good  agreement  was  achieved 
when  16  digits  beyond  the  decimal  point  were  retained,  but  when  only 
six  digits  beyond  the  decimal  point  were  retained,  considerable  error 
was  observed  in  attempting  to  reproduce  the  discretized  analytical 
output  function.  The  conclusion  from  this  exercise  was  that  an 
apparent  rounding  error  in  the  consecutive  computer  operations,  due 
to  the  small  magnitude  of  the  divisor  in  equation  (10),  led  to 
increasingly  large  magnitudes  for  the  discretized  TDT  array  values 
(i.e.,  as  the  summation  procedure  was  carried  out). 

Once  the  deflection  response  functions  had  been  obtained  as 
shown  in  Figures  26  through  3U,  the  TDT  functions  were  calculated  by 
implicit  convolution  (equation  (10))  of  the  signatures  with  their 
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respective  output  response  functions.  Figures  35  through  ^3  contain 
plots  of  the  TDT  functions  calculated  for  each  aircraft  at  each  test 
site.  Computer  sheets  containing  the  discretized  functions,  printed 
out  to  six  decimal  places,  are  provided  in  Appendix  G. 
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DISCUSSION  OF  RESULTS 

It  was  initially  hypothesized  that  a  pavement's  TDT  functions 
could  be  calculated  from  a  prime  mover's  time-dependent  deflection 
input  and  its  respective  time-dependent  deflection  outputs  at 
adjacent  locations  on  the  pavement  surface;  that  these  functions 
could  be  obtained  in  a  mathematical  sense,  without  the  need  to 
simulate  material  performance  to  determine  values  for  preselected 
descriptors;  and  secondarily,  that  these  TDT  functions  could  be 
employed  to  predict  deflection  outputs  when  the  pavement  was  subjected 
to  an  arbitrarily  imposed  loading. 

The  secondary  hypothesis  also  formed  the  basis  for  the  non- 
destructive performance  evaluation  method  presented  in  Figure  9-   The 
procedures  for  obtaining  TDT  functions  and  those  set  forth  in  the 
performance  evaluation  method  were  employed  in  conjunction  with  the 
dynamic  deflection  data  from  the  field  tests.   Analysis  of  the  data 
included: 

(1)  modeling  the  peak  dynamic  deflection  basin; 

(2)  presenting  the  measured  deflection  response  functions; 

(3)  calculating  the  TDT  functions;  and 

{h)     predicting  deflection  response  functions  for  arbitarily 
chosen  prime  movers  and  test  sites  in  accordance  with 
the  procedures  of  the  developed  performance  evaluation 
method. 
Items  (l)  and  (2)  had  to  be  accomplished  before  the  primary  hypothesis 
could  be  employed.   Item  (3)  (implicit  convolution)  provided  the 
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mechanism  to  investigate  the  validity  of  the  primary  hypothesis. 
Item  {k)    (explicit  convolution)  permitted  the  secondary  hypothesis, 
concerning  predicting  a  pavement's  performance,  to  be  explored. 

The  discussion  of  these  results  are  contained  in  this  section. 

Peak  dynamic  Deflection  Baain 

The  peak  dynamic  deflections  were  found  to  form  a  normal 

distribution  about  the  point  where  a  signature  was  applied.  The 

validity  of  this  model  was  established  by  the  correlation  coefficients 

calculated  as  part  of  the  regression  analysis  (See  Table  7)-  The 

2 

parameters  in  the  model,  0  (peak)  =*  A*e        i   ,  relate  significant 

Xi 
information  about  the  shape  of  the  deflection  basin  adjacent  to  the 

signature  point,  or  spatial  origin.  The  A  factor  is  merely  the  peak 

dynamic  deflection  of  a  signature  input,  whereas,  the  B  factor  is  a 

measure  of  the  rate  at  which  the  pavement  system  is  capable  of 

dissipating  an  energy  input,  or  a  measure  of  the  rate  at  which  the 

system  deflections  attentuate  spatially. 

From  Figure  2U,  it  is  noted  that  the  relationships  between  the 

A- factors  and  ratio  factors  were  found  to  be  linear  for  a  particular 

test  site.  A  slight  deviation  was  observed  with  respect  to  aircraft 

between  sites,  apparently  due  to  difference  in  prime  mover  weights. 

In  general,  however,  the  ratio  factors  for  each  aircraft  could  be 

taken  to  be  constant.  The  ratio  factors  used  for  the  C135,  C131, 

and  C130  aircraft  were  1.00,  0.31*,  and  O.lh,   respectively.  Once 

these  factors  were  known,  equation  (l6)  was  employed  to  calculate 

the  required  signature  array  (see  Figure  9).  This  technique  provided 
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for  the  influence  of  a  pavement's  feedback  to  an  input  signature. 

The  range  of  B  factors  for  a  prime  mover  at  a  particular  test 
site  was  observed  to  be  small  and  was  taken  to  be  constant ,  with 
exception  of  the  load  cart.  For  the  load  cart,  it  was  determined 
that  a  significant  amount  of  the  single  wheel  load  was  carried  by- 
stabilizing  wheels  located  approximately  30"  adjacent  and  to  each  side 
of  the  load  wheel.   This  was  ascertained  by  the  correspondingly  small 
B  factors  —  which  indicated  that  the  deflections  decreased  with  lateral 
distance  at  a  lesser  rate  than  for  the  other  prime  movers.  In  view 
of  this,  it  was  decided  that  the  load  cart  data  should  be  deleted 
from  further  consideration.  Figure  25  shows,  that  for  the  aircraft 
prime  movers,  a  single  B  factor  could  be  used  at  each  test  site. 

Form  of  the  Deflection  Response  Functions 

The  general  characteristics  of  the  deflection  response  functions 
for  a  particular  aircraft  were  found  to  be  similar  at  each  test  site. 
Some  observations  from  these  results  are: 

(1)  The  C135  aircraft  (twin-tandem  gear)  had  its  peak  deflection 
resultant  from  the  second  set  of  tires  in  the  gear  truck.  This 
reflected  the  pavement's  viscoelastic  behavior. 

(2)  The  C131  aircraft  displayed  only  one  peak  deflection.  This 
was  consistant  with  its  twin  gear  truck.  Some  viscoelastic  action  of 
the  pavement  was  also  evident. 

(3)  The  C130  aircraft  (single-twin  gear)  had  two  peak  deflections 
like  that  of  the  C135,  but  the  largest  peak  resulted  from  the  first 
wheel  in  the  gear  in  contrast  to  the  second  wheel  for  the  C135.  This 
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appears  to  be  due  to  the  load  distribution  on  the  gear.  Again, 
viscoelastic  behavior  of  the  pavement  was  observed. 

Form  of  the  TDT  Functions 


Figures  35  through  1*3  represent  the  TDT  functions  that  were 
calculated  by  computer  code  CONVLI  without  the  need  to  simulate 
material  performance  or  to  determine  values  for  preselected 
descriptors.  Various  characteristics  of  these  functions  are  summarized 
in  Table  8;  viz.  the  maxima,  minima,  and  zero  values  which  may  be 
thought  of  as  the  more  basic  descriptors  in  the  functions.   The 
magnitudes  of  the  maxima  and  minima  values  are  seen  to  be  independent 
of  the  aircraft.  For  example,  these  values  can  be  taken  to  be  constant 
for  all  the  aircraft  at  a  particular  site.   In  view  of  this  independence, 
the  values  can  be  thought  of  as  representing  parameters  which  indicate 
the  performance  of  the  pavement  system  itself.  It  is  also  of  interest 
that,  in  comparison  with  the  pavement  cross-sections  (see  Figure  18) , 
the  absolute  maxima  and  minima  values  increase  with  increasing  thickness 
of  pavement.  The  time  required  to  reach  sequential  zeroes  of  the  TDT 
functions  appears  to  be  independent  of  site;  e.g.  the  time  for  each 
zero  can  be  taken  to  be  constant  for  all  the  sites  for  a  particular 
aircraft.  The  zero  time  values  seem  to  represent  parameters  which 
indicate  how  each  of  the  prime  movers  influence  the  performance  of 
the  pavement  system.   In  considering  the  prime  mover  main  gear 
characteristics  (see  Figure  23),  it  is  reasonable  to  believe  that 
these  sequential  zero  time  values  represent  indicators  of  the 
influence  of  the  different  aircraft  gear  configurations.  This  finding 
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was  not  anticipated  in  the  earlier  stages  of  this  investigation. 
It  is,  however,  considered  to  be  very  significant. 

The  significant  features  of  the  TDT  functions  are  seen  to  be 

(1)  that  indicators  can  be  isolated  which  are  independent  of  the 
prime  mover  (indicators  which  reflect  the  way  a  pavement  system 
attentuates  the  energy  imparted  by  a  prime  mover's  signature);  and 

(2)  that  indicators  can  be  identified  which  represent  a  pavement's 
response  to  the  configuration  of  a  prime  mover's  input. 

Prediction  of  Deflection  Response  Functions 

The  discussion  so  far  has  indicated  the  general  validity  of  the 
concepts  associated  with  the  dynamic  deflection  basin;  time-dependent 
deflection  response  functions  (including  signatures);  and,  time- 
dependent  transfer  functions.  It  was  indicated  previously  ("Formula- 
tion of  Testing  Program"  in  the  section  on  "Predicting  Performance") 
that  once  catalogues  of  the  signatures  and  TDT  functions  of  different 
prime  movers  and  pavement  systems  were  obtained,  estimates  could  be 
had  for  combinations  of  prime  movers  and  pavements  for  which  tests 
were  not  conducted.   In  keeping  with  the  performance  evaluation 
method  as  initially  advanced  by  Figure  9»  it  was  seen  that  one 
aircraft  and  one  test  site  could  be  chosen  arbitrarily  as  standards. 
Subsequently,  predictions  of  the  deflection  response  functions  for 
the  other  two  aircraft  at  the  other  two  sites  could  be  predicted 
from  results  of  only  the  standard  vehicle  at  those  locations  (using 
computer  code  CONVLE).  Figure  hk   represents  a  sample  scheme  where 
the  C131  aircraft  and  site  2  are  taken  as  the  standards . 
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A  typical  set  of  calculated  and  measured  deflection  response 

functions  for  the  C135  and  C130  aircraft,  consistent  with  the  scheme 

in  Figure  UU,  are  shown  in  Figures  1+5  and  k6.      In  all  cases  the 

predicted  peak  deflection  values  are  seen  to  be  in  close  agreement 

with  the  measured  peak  deflection  values;  however,  the  time  history 

of  the  predicted  deflection  response  functions  are  slightly  out  of 

phase  with  the  measured  ones.  This  small  difference  was  thought  to 

be  due  to  not  taking  into  account  the  difference  in  the  sequential 

zero  time  values  for  the  various  gear  configurations.  Therefore, 

in  an  effort  to  improve  upon  the  predictions  of  the  deflection 

response  functions,  computer  code  CONVLE  was  further  refined  so  as 

to  reduce  even  these  small  differences.  This  was  done  by  using 

linear  interpolation  to  predict  a  fourth  set  of  TDT  functions  from 

the  three  known  sets  as  had  been  done  previously  for  the  signatures 

(equation  16).  Hence,  it  was  found  that  the  times  to  the  sequential 

zero  values  of  the  TDT  functions  for  the  respective  aircraft  would 

be  in  better  correspondence.  From  this,  it  follows  that  the  fourth 

set  of  TDT  functions  for,  say,  prime  mover  2  on  pavement  B  may  be 

expressed  as 

G    (first  peak) 

6   (t)  =  pXllB/^  ; rr  •  G   <*>     <17) 

xi2B      G    (first  peak)    x±2k 
ilA 

A  typical  set  of  calculated  and  measured  deflection  response  functions 

for  the  C135  and  C130  aircraft,  consistent  with  the  scheme  in  Figure 

UU,  are  shown  in  Figures  hj   and  U8,  respectively.  In  all  cases  the 

calculated  deflection  response  functions  correspond  well  in  time  to 

the  measured  functions,  and  in  all  cases  the  calculated  peak  deflections 
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FIGURE    45.    TYPICAL    CALCULATED   AND   MEASURED 
DEFLECTION  RESPONSE   FUNCTIONS   FOR 
CI35    USING   PREDICTED   TDT  FUNCTIONS 
WITH   CI3I  AS  STANDARD   PRIME  MOVER  . 
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FIGURE    46. 


TYPICAL  CALCULATED  AND  MEASURED 
DEFLECTION  RESPONSE  FUNCTIONS  FOR 
C 135, USING  TDT  FUNCTIONS  FROM  CI3I 
WITH  CI3I  AS  STANDARD  PRIME  MOVER. 
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were  within  10?  of  the  measured  ones.  A  schematic  of  the  revised 
performance  evaluation  method  is  provided  as  Figure  U9. 

Just  as  various  forms  of  the  constitutive  relations  may  change 
outputs  for  a  given  input,  it  could  he  expected  that  the  TDT  functions 
could  also  reflect  these  characteristics.   Radical  differences  between 
a  standard  vehicle  and  a  prime  mover  in  question  could  lead  to  large 
decrepancies  between  predicted  and  observed  deflections.   The  two 
phenomena  responsible  for  this  are  the  relative  contact  area  of  the 
load  inputs  and  the  comparative  stress  levels  of  the  vehicles;  these 
should  be  matched  to  some  degree  to  insure  good  predictions.   It  is 
not  believed,  however,  that  these  limitations  are  restrictive  in  the 
world  of  real  aircraft  and  real  pavements. 

Some  additional  plots  of  measured  and  calculated  deflection 
response  functions  for  other  combinations  of  the  scheme  depicted  in 
Figure  hk   are  given  in  Appendix  H.   Case  1  in  Appendix  H  represents 
the  scheme  shown  in  Figures  kh   through  1*8.   Cases  2  and  3  use  site  2 
as  standard  with  C135  and  C130,  respectively,  as  standard  prime 
movers.   Case  h   uses  the  C130  and  site  1  as  standards. 
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SUMMARY  AND  CONCLUSIONS 

Transfer  function  theory  in  connection  with  discretized  real 
convolution  can  be  extended  to  define  a  pavement's  time-dependent 
transfer  functions.   The  time-dependent  transfer  functions  can  be 
calculated  from  an  implicit  solution  of  the  convolution  integral 
using  data  measured  in  a  field  environment.      Dynamic  deflection 
response  functions   can  be  calculated  to  predict  pavement  performance 
by  employing  the  convolution  procedure  in  an  explicit  manner. 

A  procedure  was  advanced  to  predict  a  pavement's  performance 
for  a  particular  load  or  situation  which  had  not  been  previously 
experienced.     As  an  example,  two  sites   (labled  site  1  and  site  2) 
and  two  prime  movers   (A  and  B)   are  involved.     One  site  and  one  prime 
mover  are  considered  to  be  standard,   say  site  1  and  prime  mover  A. 
It  is  desired  to  predict  the  pavement's  performance  at  site  2 
incident  to  prime  mover  B.     Signatures  and  time-dependent  transfer 
functions  are  catalogued  for  both  prime  mover*  at  the  standard  site. 
The  standard  prime  mover  is  then  run  at  site  2  so  that  its  signature 
and  time-dependent  transfer  functions  may  be  calculated.     The 
relationships  between  signatures  and  time- dependent  transfer 
functions  at  the  standard  site  and  the  signature  and  time-dependent 
transfer  functions   from  the  standard  prime  mover  at  site  2   (site 
under  question)  permit  prediction  of  the  deflection  response 
functions   for  prime  mover  B  at  site  2  to  be  carried  out. 

To  validate  the  transfer  function  approach  and  the  prediction 
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method,  full  scale  field  tests  were  conducted.  The  time-dependent 
transfer  functions  were  calculated  for  three  different  pavement 
cross-sections  using  the  data  collected  for  three  operational 
military  aircraft. 

The  results  of  these  tests  indicate  the  following  conclusions: 

(1)  There  is  a  unique  set  of  time-dependent  transfer  functions 
at  a  particular  site  for  a  particular  prime  mover.  These  functions 
do  appear  to  provide  the  basis  for  a  proposed  performance  evaluation 
method. 

(2)  The  peak  deflections  form  a  normal  distribution  about 

the  point  where  a  signature  is  applied.  Excellent  correlation 

between  the  measured  data  and  the  normal  distribution  model 

0  (peak)aA  e~  *     i   was  achieved.  The  A  parameters  were 

xi 
found  to  be  dependent  on  the  pavement  and  prime  mover;  whereas, 

the  B  parameters  were  taken  to  be  constant  for  a  given  pavement, 

independent  of  prime  mover. 

(3)  The  time-dependent  transfer  functions  are  independent  of 
magnitude  of  load  input  but  reflect  to  a  minor  extent  some  characteristics 
of  the  prime  mover  gear  configurations . 

(1*)  In  philosophy,  the  transfer  function  provides  a  measure 
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of  "global  properties"   of  a  pavement   in   contrast   to  the   "point-wise 
aspects"   inherent   in   classical   (elasticity  and  viscoelasticity) 
theories. 

(5)      The  time-dependent   characteristics  of  a  flexible 
pavement   can  be  represented  by  a  series   of  time-dependent   transfer 
dynamic  tests  using  a  common  prime  mover. 
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SUGGESTIONS  FOR  FUTURE  RESEARCH 

This  investigation  has  been  a  first  step  in  formulating  and 
testing  the  concept  of  using  transfer  function  theory  to  evaluate  a 
pavement  system.   Conclusions  from  the  full  scale  field  tests,  although 
limited  in  scope ,  clearly  establish  the  validity  of  the  hypotheses  as 
initially  conceived.  In  view  of  the  favorable  results  and  the 
experience  gained  by  conducting  this  investigation,  the  following 
recommendations  for  further  research  are  provided: 

1.  The  transfer  function  approach  advances  a  global  concept 
wherein  the  TDT  functions  are  sensitive  to  the  distance 
over  which  the  input  and  output  functions  are  measured. 
A  more  reliable  deflection  measuring  system  should  be 
developed  to  measure  the  time-dependent  input  and  output 
functions  relative  to  the  prime  mover  instead  of  at 
fixed  locations  on  the  pavement  surface.   Such  a  system 
would  eliminate  the  need  to  shift  the  deflection  response 
functions,  and  thereby  reduce  any  errors  associated  with 
the  modeling  used  in  the  shifting  procedure. 

2.  Develop  a  vibration  testing  device  capable  of  determining 
a  pavements  transfer  functions  by  the  frequency  spectrum 
technique.  The  device  must  be  capable  of  delivering  a 
constant  input  force  between  5  and  150  Hz.  Successful 
employment  of  this  technique  would  negate  use  of  a 
standard  vehicle. 

3.  Perform  a  study  to  isolate  and  identify  system 
parameter  contained  in  the  TDT  functions.  These 
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parameters  could  then  be  correlated  with  current 
theoretical  models  and  existing  evaluation  and 
design  methods. 

k.     Conduct  a  parameter  study  to  incorporate  a  greater 
number  of  variables,  including  temperature,  speed, 
environmental,  and  geographic  conditions,  etc. 
Application  and  correlation  of  the  concept  to  layered 
systems  and  different  materials,  such  as  soil  and 
concrete,  should  be  considered. 

5.  Initiate  a  program  to  monitor  TDT  functions  over  a 
long-time  period  to  permit  correlation  of  parameters 
with  pavement  deterioration.  Emphasis  could  be 
placed  on  establishing  a  distress  criteria  and 
standards  for  maintenance  and  repair  of  pavements . 
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APPENDIX  A 

METHODS  FOR  DETERMINING  TRANSFER  FUNCTIONS 

"A  Priori"  Modeling 

When  "a  priori"  modeling  is  employed,  a  transfer  function  is 
normally  obtained  through  a  three-stage  process.  First,  a  free- 
body  diagram  is  constructed  to  model  the  real  system;  secondly,  an 
appropriate  lav,  or  principle,  is  invoked  (e.g.,  Newton's  Second 
Law  of  Motion,  Conservation  of  Energy,  etc.)  to  obtain  the  governing 
(differential)  equations;  and  lastly,  the  (differential)  equations  are 
solved  so  that  the  ratio  of  input  to  output  can  be  formed.  As  an 
example  of  the  process,  consider  the  schematic  representation  of  a 
simple  model  of  a  pavement  system  as  shown  in  Figure  50.  Assuming 
motion  to  be  in  the  vertical  direction  only,  applying  Newton's  Second 
Law  of  Motion,  the  governing  differential  equation  of  motion  becomes 

mx(t)  +  cx(t)  +  kx(t)  =  f(t)  (18) 

Solving  equation  (18)  by  Laplace  transform  method  (with  initial 
conditions  equal  zero)  results  in  (reference  55) 

x(s)  =  — — •  ?(s)         (19) 

ms  +  cs  +  k 

or, 

Si*!  s   G(s)  =  —g—i (20) 

f(s)         ms  +  cs  +  k 


F(t) 
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FIGURE   50.    SIMPLE   SECOND  ORDER   MODEL 
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FIGURE  51.      ASYMPTOTIC    APPROXIMATION    OF  A 
SYSTEM'S  FREQUENCY  SPECTRUM 
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where , 

G(s)  =  transfer  function,  with  displacement  of  the  mass  as  the 
operational  output  and  the  forcing  function  as  the 
operational  input. 

"A  Posteriori"  Modeling 

In  "a  posteriori"  modeling,  methods  for  determining  a  transfer 
function  are  basically  experimental.  Three  commonly  used  methods  are 
presented  below: 

Frequency  Spectrum 

Suppose  a  sinusoidal  input  of  known  magnitude  and  frequency  is 
applied  to  a  system  whose  transfer  function  is  sought.  The  magnitude 
of  the  steady-state  peak  input  and  the  steady-state  peak  output  (and 
generally  the  phase  angle  between  the  two)  are  measured.  The  ratio 
of  the  magnitude  of  the  peak  output  to  that  of  the  peak  input  is 
plotted  for  a  range  of  frequencies.   The  resulting  plot  is  called 
the  "frequency  spectrum"  for  the  system.  The  frequency  range  need 
cover  only  those  frequencies  wherein  the  slope  of  the  frequency 
spectrum  curve  possesses  a  changing  slope.  The  curve  is  then 
approximated  by  a  series  of  straight  lines  as  shown  in  Figure  51.  The 
straight  lines  are  selected  so  as  to  provide  asymptotes  to  the  curve. 
The  closer  the  approximation  to  the  actual  frequency  spectrum  curve, 
the  better  will  be  the  approximation  of  the  transfer  function  (reference 
51).   The  intersections  of  the  selected  approximating  asymptotes 
determine  the  "break"  frequencies,  or  "corner"  frequencies.  In  Figure 
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51,   for  example,   there  are   five;   shown  as    w    =  0,  u>  ,   w  ,   w,    and  w  . 
The  slopes  of  the  asymptotes  are  given  as  K   ,  K?,  K   ,  K,    and  K   , 
respectively.      Once  the  break  frequencies  and  the  slopes  of  the 
approximating  asymptotes  are  ascertained,   the  modified  transfer 
function  which  describes   the   frequency  spectrum  can  be  expressed 
(reference   51)    as 

G(ju)   =  AUid+oo^I   •   (Ju>+U)2)K2~K1 (ju*a)n)Kn"Kn-l        (21) 

Substituting  for  the  complex  variable  s  =  ju)  in  equation   (21) 
the  experimental  transfer  function  for  the  system  is  obtained  as 

G(s)   =  A(s+coJKl   •    (s+wJK2~Kl (s+co   )Kn~Kn-l  (22) 

x         d  n 

The  constant  A  can  be  determined  from  any  one  experimental  frequency 
in  the  frequency  spectrum  by  calculating  the  absolute  value  of  the 
transfer  function  from  equation  (21)  at  that  frequency  and  equating  it 
to  the  corresponding  ratio  of  the  magnitude  of  the  peak  output  to  that 
of  the  peak  input  in  the  frequency  spectrum. 

Impulse  Input 

The  mathematical  nature  of  an  impulse,  or  Dirac  delta  function 
input  resolves  the  transfer  function  for  this  system  directly  as  the 
corresponding  output.  This  input  is  characterized  mathematically  as 


,  1/e,  0  <  t  <  e 
6(t)  =  ""  F  (t)  =  \  (23) 

e"°  e      »  0,  t  >  e 
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where , 

<5(t)     =  impulse,  or  Dirac  delta  function 

F  (t)  ■  Heaviside's  unit  function 

It  can  also  be  shown  (reference  5M   that  the  Laplace  transform 
of  the  Dirac  delta  function  is  unity,   or 


1 


[6(t)]  =  1  (2k) 


Consequently  from  the  definition  of  a  transfer  function  it  is  readily 
seen  that  for  an  impulse  input 

0(s)  =  G(s)  (25) 

Therefore,  theoretically  all  that  need  be  done  to  determine  a  system's 
transfer  function  when  the  system  is  excited  by  an  impulse  input  is 
to  measure  its  output  response. 

Real  Convolution 

Both  the  frequency  response  and  impulse  input  methods  for 
obtaining  a  transfer  function  are  frequency,  or  Laplace,  domain 
approaches:  inversion  of  the  solution  is  required  in  order  to  get  a 
time  domain  solution.  Real  convolution  is  essentially  a  time  domain 
approach.   Since  it  is  a  time  domain  response  technique,  convolution 
eliminates  the  transform  and  inverse  transform  steps  associated  with 
the  applications  of  transfer  functions.  Basically,  the  analytical 
concepts  provided  by  convolution  are:   (l)  proof  of  multiplicative 
nature  of  transfer  functions,  and  (2)  duality  between  multiplication 
in  frequency  (or  Laplace)  domain  and  convolution  in  the  time  domain. 
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Of  equal  importance  with  the  concepts  identified  above,  convolution 
can  be  used  as  a  numerical  and/or  graphical  tool  for  computation  of 
responses  without  resorting  to  analytical  expressions  for  input. 

Convolution  introduces  the  concept  of  a  "system  response  function" 
A  system's  transfer  function  is  defined  as 

0(s)  =  G(s)  •  T(s)  (26) 

where , 

0(s)  =  system  operational  output 

I(s)  =  system  operational  input 

G(s)  =  system  transfer  function 
Taking  the  inverse  Laplace  transform  of  equation  (25),  it  follows  that 


0(t)  =^_1  [G(s)  '  I(s)]  (27) 


Equation  27  can  then  be  written  in  the  form  (references  5^,  58,  6l, 
and  62)  of  the  real  convolution  integral  as 

0(t)  =  /  G(t)  •  I  (t-T)  •  dx        (28) 
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APPENDIX  B 

REGRESSION  ANALYSIS-NORMAL  DISTRIBUTION  MODEL 
COMPUTER  PROGRAM 


PROGRAM  NMLFIT( INPUT, OUTPUT, TAPE5=INPUT,TAPE6=OUTPUT) 

C 

C PROGRAM  NMLFIT 

C 

C      

C     PURPOSE 

C        TO  COMPUTE  PARAMETER  VALUES  IN  NORMAL 
C      DISTRIBUTION  MODEL  FOR  PEAK  DYNAMIC  DEFLECTION 
C      BASIN  USING  LEAST  SQUARES  REGRESSION  ANALYSIS 
C     DESCRIPTION  OF  PARAMETERS 

C      X   -VECTOR  OF  VALUES  WHERE  DEFLECTIONS  WERE  MEASURED  (INCHES) 
C      Y   -VECTOR  OF  DEFLECTION  VALUES  TO  BE  FITTED  BY  MODEL 
C      Z   -VECTOR  OF  FITTED  DEFLECTION  VALUES 
C      J   -DATA  SET  NUMBER. 
C      K   -NUMBER  OF  DATA  SETS. 
C      M   -SETS  OF  DATUM  IN  X-Y  ARRAYS. 
C      XMAX-MAXIMUM  VALUE  IN  X  ARRAY. 
DIMENSION  Y(200),  Z(200),  W(200) 
COMMON/AXIS /X( 200 ) 
9  READ(5,100)  J,K,M,XMAX 
100  F0RMAT(2I1,I2,F10.5) 
DO  8  1=1, M 
8  READ(5,102)  Y(I),X(I) 
102  FORMAT(32X,2F8.5) 
DO  36  1=1, M 
36  WRITE(t,110)  X(I),Y(I) 
110  FORMAT(2F10.6) 
F=3.0/XMAX 
DO  30  1=1, M 
30  X(l)=X(l)*F 

CALL  NMLFII(Y,Z,M) 
WRITE  (6,101+) 
10l*  FORMAT(lHl,2UX,7*Xy,10X,#:?0 
DO  35  1=1,72 
X(I)=FL0AT(I) 
C     PRINT  OUT  FITTED  DATA. 
35  WRITE(6,105)  X(I),Z(I) 
105  FORMAT (20X,2F10.  5) 
IF(J.NE.K)  GO  TO  9 
STOP 
END 
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SUBROUTINE  NMLFIl(Y,Z,JJJ) 
C 

C  ...  SUBROUTINE  NMLFII 
C 

C  

C  PURPOSE 

C        LEAST  SQUARE  METHOD  BY  USING  NEC  EXP.  FUNCTION  TO  FIT  VECTOR 

C      OF  VALUES  Y(I)  WITH  THE  HELP  OF  NEWTON  ITERATION  METHOD  TO 

C      SOLVE  THE  ROOT  OF  FUNCTION  GB. 

C  DESCRIPTION  OF  PARAMETERS 

C      X  -THE  AXIS  AT  WHICH  Y(l)  IS  TABULATED. 

C      Y  -VECTOR  OF  VALUES  TO  BE  FITTED  BY  NEG.  EXP.  FUNCTIONS. 

C      Z  -THE  CALCULATED  VALUES  OF  NEG.  EXP.  FUNCTION  ON  X(l) 

C      JJJ-THE  DIMENSION  OF  X,Y,Z. 

C  SUBROUTINE  AND  FUNCTION  SUBPROGRAM  NEEDED. 

C      FUNCTION  GB  AND  SUBROUTINE  NEWTON. 

C  METHOD. 

C      1.  GB  IS  THE  FUNCTION  DEDUCED  FROM  THE  LEAST  SQUARE  METHOD  BY 

C         USING  NEG.  EXP.  FUNCTIONS.   IT  ALSO  COMPUTES  SOME  DERIVATIVES 

C         VALUES  OF  GB  WITH  RESPECT  TO  EXP(-B). 

C      2.  CALL  SUBROUTINE  TO  FIND  OUT  THE  ROOT  OF  GB(EXP(-B))  =  0. 

C      3.  THEN  COMPUTES  A  AND  Z(l). 

C 

C  

c 

COMMON /AXIS /X( 200) 

DIMENSION  Y(l),Z(l),XGB(200) 

J=JJJ 

PRINT  5 

WRITE(6,212) 
212  F0RMAT(1H1,1X,*REC0RD  OF  SMOOTH  OPERATIONS  *//) 

DX=.01 

XM=l.E-7 

XGB(101)=GB(X,Y,XM,J) 

XM=l.E-lU 

XGB(102)=GB(X,Y,XM,J) 

XM=1.E-19 

XGB(103)=GB(X,Y,XM,J) 

DO  UO  1=1,100 

XM=FL0AT(I)«DX 

XGB(I)=GB(X,Y,XM,J) 
1*0  CONTINUE 

PRINT  5 

PRINT  30,(XGB(K),K=1,103) 

CALL  NEWTON(X,Y,ROOT,J) 

IF(ROOT.EQ.O.)  GO  TO  150 
310  B=-ALOG(ROOT) 

VX=0. 

VX=VY=0. 

DO  20  K=2,J 

VX=VX+Y(K)*ROOT**(X(K)*X(K)) 
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VY=VY+R00T**(X(K)*X(K)*2. ) 
20  CONTINUE 

A=VX/VY 

PRINT  5 

PRINT  100 , A, B 
100  FORMAT(*0   A=*,F12.6,*   B=*,F12.6) 

GO  TO  200 
150  PRINT  90 

90  FORMAT(*0   ERROR  IN  SUBROUTINE  ROOTFINDER*) 
200  CONTINUE 

PRINT  6 

J=72 

DO  25  1=1, J 

X(l)=FLOAT(l-l)*(3./72.) 

Z(I)=A*(R00T**(X(I)**2)) 
25  CONTINUE 

PRINT  30,(Z(I),I=1,J) 

PRINT  6 

6  F0RMAT( *0  ====================================================*) 

30  FORMAT (lH0,10F12. 6) 
RETURN 
END 


FUNCTION  GB(X,Y,BX,J) 
C 

C  ...  FUNCTION  GB 
C 

C  

C  PURPOSE  AND  METHOD. 

C  GB  IS  THE  FUNCTION  DERIVED  FROM  THE  LEAST  SQUARE  METHOD  FOR  NEC, 

C  EXPONENTIAL  FUNCTION  A*EXP(-B*X*X)  =  A*(BX**(X*X) ) ,  WHERE 

C  BX=EXP(-B).   DERIVE  TWO  EQUATIONS  FOR  TWO  UNKNOWNS  A  AND 

C  BX.   ELIMINATE  A  AND  OBTAIN  A  FUNCTION  ONLY  DEPENDING 

C  ON  BX.   THIS  IS  THE  GB  FUNCTION. 

C  SUBROUTINE  AND  FUNCTION  SUBPROGRAM  DEEDED. 

C  NONE. 

C  DESCRIPTION  OF  PARAMETERS  (SAME  AS  IN  NMLFII  J=JJJ) 

C  BX  -  THE  INDEPENDENT  VARIABLE  OF  GB.   IT  RANGES  FROM  0.  TO  1. 


C 
C 

c 


DIMENSION  X(l),Y(l) 
XYA=XYB=XYC=XYD=0 . 
XYE=XYF=0. 
L=2 

IF(X(l).NE.O.)  L=l 
DO  100  I=L,J 
XX=X(I)*X(I) 
IF(XX.EQ.O.)  GO  TO  200 
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YY=Y(l)*BX**XX 

ZZ=(BX**XX)*(BX**XX) 

GO  TO  150 
200  YY=Y(l) 

ZZ=1. 
150  CONTINUE 

XYA=XYA+YY*XX 

XYB=XYB+XX*ZZ 

XYC=XYC+YY 

XYD=XYD+ZZ 

XYE=XYE+YY*XX*XX 

XYF=XYF+ZZ*XX*XX 
100   CONTINUE 

GB=XYA/XYB-XYC/XYD 

X(195)=XYA 

X(196)=XYB 

X(197)=XYC 

X(198)=XYD 

X(199)=XYE 

X(200)=XYF 

RETURN 

END 


SUBROUTINE  NEWTON (X,Y, ROOT, J ) 
C 

C    . . .    SUBROUTINE  NEWTON 
C 

C  

C  PURPOSE  AND  METHOD 

C      USE  NEWTON  RALPSON  METHOD  TO  SOLVE  FOR  THE  ROOT  OF  GB  FUNCTION 

C      WITH  THE  INITIAL  VALUE  IN  UU. 

C      WITH  THE  INITIAL  VALUE  IN  UU,  DERIVATIVE  OF  GB  IS  CALCULATED  IN 

C      GB  FUNCTION  ROUTINE.   THE  NO.  OF  ITERATIONS  (ICOUNT)  CANNOT 

C      EXCEED  OVER  150  TIMES.   THE  ROOT  WILL  BE  RETURNED  IN  ROOT. 

C  SUBROUTINE  AND  FUNCTION  SUBPROGRAM  DEEDED. 

C      GB  FUNCTION  ROUTINE. 


C 
C 

c 


DIMENSION  X(l),Y(l) 
UU=. 99999999 
UU=.6 

TEST=l.E-7 
ICOUNT=0 
10  GUU=GB(X,Y,UU,J) 

DUU»X(199)/X(196)-2.«X(200)»X(195)/(X(196)»X(196))-X(195)/X(198) 

1-2.*X(196)»X(197)/(X(198)*X(198)) 

IF(DUU.EQ.O.)  GO  TO  150 

IF(ABS(GUU).LT.TEST)  GO  TO  100 

U=UU*(1.-GUU/DUU) 
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IF(U.LT.O. )   U=UU/2. 

IF(ABS(U-UU).LT.TEST)   GO  TO  100 

IF(IC0UNT.GE.150)   GO  TO  160 

IC0UNT=IC0UNT+1 

UU=U 

PRINT  20,IC0UNT,GUU,U 

GO  TO  10 
100  ROOT=U 

RETURN 
150  PRINT  5 

5  FORMAT(^0  ***  ZERO  DERIVATIVE  OF  GB .EXITS  j) 
STOP 

160  PRINT  6 

6  FORMAT(?<0  ***  ITERATIONS  OVER  150  TIMES ,EXITS  jO 
20   FORMAT(*  +++»,lU,«  +++*,2F12.6) 

GO  TO  100 
END 
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APPENDIX  C 

SIGNATURE  AND  TIME-DEPENDENT  TRANSFER  FUNCTIONS 
COMPUTER  PROGRAM 


PROGRAM  CONVLI  ( INPUT .OUTPUT ,TAPE5=INPUT , TAPE 6= OUTPUT ,PLOT) 

DIMENSION  XINP(150),XOUT(150),XG(150),W(10) 

DIMENSION  DEFL(6,500) ,Z(500) ,GT(6,500) ,Y( 500) ,GF(6,500) 

COMMON /AXIS /X( 500) 

DIMENSION  RDD(lO),RDN(lO),RDF(lO,150),M(lO) 

DIMENSION  IPS(9),NAME(10),PEAK(12),C0NT(10),M0NT(10),SPK(10) 

DIMENSION  Vl(200) ,V2(200) ,V3(200) ,VU(200) ,V5(200) ,V6(200) 

INTEGER  CONT, BEGIN, END 

REAL  MOD.INCR 

C     READ  DATA  BEGIN  AND  END,  DATA  SET  NUMBERS,  NUMBER  OF  SEQUENCE 

C  VALUES,  DITIGIZED  TIME  INCREMENT,  PLOT  COMMAND,  NORMAL 

C  MODEL  PARAMETER  VALUES,  TIME  RANGE,  NUMBER  OF  DESIRED 

C  ARRAY  VALUES,  AND  COMMON  ORIGIN  COMMAND 

READ(5,100)  BEGIN,  END,  N,  DT,  IPLOT,  A,  B,  RANGE,  NY,  IMOD 

N1=N 

IF ( IPLOT. NE.O)  CALL  PLOTS 

CALL  FACTOR (0.6) 

INS=2 

INL=6 

NMY=NY-1 

DY=RANGE / FLOAT ( NMY ) 

DO  UlO  K=1,NY 
UlO  Y(K)=FLOAT(K-l)*DY 
C     BEGIN  LOOP  TO  CALCULATE  SIGNATURE  AND  TDT  FUNCTIONS 

DO  1000  J=BEGIN,  END 
C     READ  IDENTIFICATION  AND  DATA 

WRITE(6,208)  J 

READ  (8,106)  (NAME(I), 1=1,9) 

WRITE(6,107)  (NAME(L),L=1,9) 

DO  30  1=1,6 
30  READ(8)      (DEFL(I,K),K=1,N) 

READ  (8)  ICC 

IF(F0F,8)    U05.50 
1+05  CONTINUE 
C 

C     FOR  EACH  RESPONSE  FUNCTION- 

C        (1)  MODIFY  RESPONSE  FUNCTIONS  TO  OBTAIN  COMMON  ORIGIN (IF  NECESSARY) 
C        (2)  DETERMINE  PEAK  VALUE  OF  EACH  RESPONSE  FUNCTION 
C        (3)  STORE  SEQUENCE  NUMBER  OF  EACH  PEAK  VALUE 
C 

DO  396  1=1,6 

PEAK(I)=0. 

IF(IMOD.EQ.O)  GO  TO  391* 
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AVE=0. 

DO  391  K=l,10 
391  AVE=AVE+DEFL(I,K) 

MOD=AVE/lO. 
393  DO  399  K=1,N 
399  DEFL(l,K)=DEFL(l,K)-MOD 
39U  DO  395  K=1,N 

IF(ABS(PEAK(I)).LT.ABS(DEFL(I,K)))  CONT(l)=K 

395  IF(ABS(PEAK(I)).LT.ABS(DEFL(I,K)))  PEAK(l)=DEFL(l,K) 

396  CONTINUE 
C 

C     FOR  PRIMARY  RESPONSE  FUNCTION  (I.E. ,DEFL(2,K) )- 

C        (l)ESTABLISH  PERCENT  OF  PEAK  VALUE  TO  RETAIN  SEQUENCE  VALUES 

C        (2)TEST  SEQUENCE  VALUES  VERSUS  PERCENT  OF  PEAK(RETAIN  .GT. ) 

C        (3)ZERO  ALL  SEQUENCE  VALUES  LESS  THAN  PERCENT  OF  PEAK  VALUE 

C        (It) REDEFINE  N 

C 

GAP=.002 

DO  273  L=1,N 

IF(ABS(DEFL(INS,L)).GT.GAP.AND.L.GT.2)  GO  TO  290 

LI=L 

DEFL( INS ,1 )=DEFL( INS ,LI ) 
273  CONTINUE 
290  DEFL(INS,2)=DEFL(INS,LI+1) 

LK=C0NT(INS)-LI 

N=3*LK+3 

DO  397  K=3,N 

DEFL(INS,K)=DEFL(INS,K+LI-1) 

IF(DEFL(lNS,K).LT.O.)  DEFL(lNS,K)=0. 

397  CONTINUE 
INB=INS+1 

DO  299  I=INB,INL 

DO  299  K=1,N 

DEFL(I,K)=DEFL(I,K+LI-1) 

IF(DEFL(l,K).LT.O.)DEFL(l,K)=0. 
299  CONTINUE 

NK=N+1 

DO   295  I=INS,INL 

DO  295  K=NK,N1 

DEFL(I,K)=0. 
295  CONTINUE 

DO  298  1=1,6 

PEAK(l)=0. 

DO  298  K=1,N 

IF(ABS(PEAK(I)).LT.ABS(DEFL(I,K)))    CONT(l)=K 

IF(ABS(PEAK(I)).LT.ABS(DEFL(I,K)))   PEAK(l)=DEFL(l,K) 
298  CONTINUE 
C      EXTEND  OR   COMPRESS   DEFL  ARRAY  TO   LENGTH  =  NY 

DZoRANGE / FLOAT ( N-l ) 

IF(N.EQ.NY)   GO  TO  U29 

DO   Itll  K«1,N 
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1+11  Z(K)=FLOAT(K-l)*DZ 
DO  1+20  1=2.6 
IPT=1 

GT(I,1)=DEFL(I,1) 

GT(I,NY)=DEFL(I,N) 

DO  1+21  K=2,NMY 
1+25   IF(Y(K)    .LE.Z(IPT+1))    GO  TO   1*2 U 

IPT=IPT+1 

IF(lPT.GT.N+5)    GO  TO   50 

GO  TO  1+25 
1+21+   0MG=(Y(K)-Z(IPT))/DZ 

1+21  GT(I,K)=0MG«DEFL(I,IPT+1)  +  (1.-0MG)*DEFL(I,IPT) 
1+20  CONTINUE 

DO   1+26  1=2,6 

DO   1+27  K=1,NY 

DEFL(I,K)=GT(I,K) 
1+27  GT(I,K)=0. 

DEFL( I ,NY+1 )=DEFL( I ,NY+2 )=DEFL( I ,NY+3 )=DEFL( I ,NY ) 
1+26  CONTINUE 

N=NY 
1+29   CONTINUE 
1+31  CONTINUE 
C 
C  PRINT  NEW  TIME-DEPENDENT  DEFLECTION  RESPONSE   FUNCTIONS(TRUNCATED) 

C 

JL=0 

1098  DO  1099  I=INS,INL 
W(I)=FL0AT(I-INS)*.75 
RDD(l)=A*EXP(-B*W(l)**2) 
RDN(I)=RDD(I)/PEAK(I) 
RDF(I,1)=RDN(I) 

DO  1100  K=1,N 
1100  DEFL(I,K)=DEFL(I,K)*RDN(1) 
JK=100*J 
DO  800  K=1,N 

800  GT(I,K+JK)=DEFL(I,K) 

1099  CONTINUE 
IF(JL.EQ.l)  GO  TO  802 
IF(J.NE.3)  GO  TO  1001 
DO  801  I=INS,INL 

DO  801  K=1,N 
DEFL(l,K)=(GT(l,K+100)+GT(l,K+200)+GT(l,K+300))/3. 

801  CONTINUE 
JL=1 

DO  805  I=INS,INL 
PEAK(I)=0. 
DO  805  K=1,N 

IF(ABS(PEAK(I)).LT.ABS(DEFL(I,K)))  C0NT(l)=K 
IF(ABS(PEAK(I)).LT.ABS(DEFL(I,K)))  PEAK(I )=DEFL(l ,K) 
805  CONTINUE 
GO  TO  1098 
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802  CONTINUE 

WRITE(6,213) 
WRITE(6,2lU) 
WRITE (6, 215) 
213  F0RMAT(1H1,1X/////,12X*  DEFLECTION  RESPONSE  F 

•UNCTIONS') 

C     INPUT  AIRCRAFT  AND  SITE  IDENTIFICATION 

2lU  FORMAT (38X,'C135  AIRCRAFT-SITE  1') 

215  F0RMAT(1X//,'  TIME      RESPONSE  0   RESPONSE  1   R 
*ESPONSE  2   RESPONSE  3   RESPONSE  U ' ) 

DO  1103  K=1,N 

DEFL(1,K)=0. 
1103  PRING  201,Y(K),(DEFL(l,K),I=2,6) 

DO  llOU  K=1,N 
llOU  PUNCH   108,(DEFL(l,K),l=2,6) 

500  DO   501  1=1,6 
DO   501  K=1,N 

501  GF(I,K)=0.0 
DO   kOk  K=1,N 

kOk   XINP(K)=DEFL(INS,K) 

INB-INS+1 

DO    1+01  I=INB,INL 

DO   1+02  K=1,N 
1+02  XOUT(K)=DEFL(l,K) 

CALL  C0NVL1(XINP,X0UT,XG,DY,N,PEAK(2)) 

DO  1+03  K=1,N 
1+03  GF(I,K)=XG(K) 
1+01   CONTINUE 

WRITE (6,216) 

WRITE( 6,217) 

WRITE(6,2l8) 

216  F0RMAT(1H1,1X/////,'  TDT   FUNCTI 
*ONS ' ) 

C     INPUT  AIRCRAFT  AND  SITE  IDENTIFICATION 

217  FORMAT(32X,'C135  AIRCRAFT-SITE  1') 

218  FORMAT (1X//,'  TIME      FUNCTION  1  FUNCTION  2   FU 
*NCTION  3  FUNCTION  1+*) 

DO  502  K=1,N 

502  PRINT  201,Y(K),(GF(I,K),1=3,6) 
DO  503  K=1,N 

503  PUNCH  108,(GF(I,K),I=2,6) 
IF(IPLOT.EQ.O)  GO  TO  1001 

C 

C     PLOT  TDT  FUNCTIONS  AND  SIGNATURE 

C 

DO  530  K=1,N 

V2(K)=0.0 

V3(K)=0.0 

Vl+(K)=0.0 

V5(K)=0.0 

V6(K)=0.0 
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530  CONTINUE 
C 

C     DETERMINE  SCALING  PARAMETERS  FOR  TDT  FUNCTIONS 

C 

DO  506  K=1,N 

V2(K)=GF(2,K) 

V3(K)=GF(3,K) 

VU(K)=GF(1+,K) 

V5(K)=GF(5,K) 

V6(K)=GF(6,K) 
506  CONTINUE 

CALL  SCALES(l0.0,V2,N,l.V3,N,l.VU,N,l,V5,N,liV6,N,l) 

CALL  SCALE(X,10.0,N,l) 
C 

C     PLOT  AXES  AND  ANNOTATION 

C 

CALL  AXIS(0.0,0.0,12HTIME, SECONDS, -12, 10. 0,0.0,X(N+l),X(N+2),0) 

CALL  AXIS(0.0,0.0,1THRESP0NSE  FUNCTION, IT,  10.0, 90. 0, V2(N+l) ,V2(N+2 

CALL  SYMBOL(3.5,10.0,0.2,l8HRESPONSE  FUNCTIONS, 0.0,18) 

CALL  SYMBOL( 8. 0,8.2,0.08, 1,0.0,-1) 

CALL  SYMBOL(8.1,8.lU,0.08,llH-RESPONSE  1,0.0,11) 

CALL  SYMB0L( 8. 0,8. 1,0.08,2,0. 0,-1) 

CALL  SYMBOL(8.1,8.0U,0.08,11H-RESPONSE  2,0.0,11) 

CALL  SYMBOL(8.0,8.0,0.08,3,0.0,-l) 

CALL  SYMBOL(8.1,7-9H,0.08,11H-RESPONSE  3,0.0,11) 

CALL  SYMB0L(8. 0,7. 9, 0.08,U, 0.0,-1) 

CALL  SYMB0L(8.1,7.8U,0.08,11H-RESP0NSE  U, 0.0,11) 

CALL  SYMB0L(8. 0,7. 8,0.08,5,0.0,-1) 

CALL  SYMBOLCS.l^^.O.OS.llH-RESPONSE  5,0.0,11) 

CALL  SYMBOL(8.0,7.7,0.08,6,0.0,-l) 

CALL  SYMBOL(8.1,7.6i*,0.08,llH-RESPONSE  6,0.0,11) 
C 

C     PLOT  TDT  FUNCTIONS 

C 

DO  517  I=INB,INL 

DO  509  K=1,N 
509  Z(K)=GF(I,K) 

Z(N+1)=V2(N+1) 

Z(N+2)=V2(N+2) 

CALL  LINE(X,Z,N, 1,10,1) 
517  CONTINUE 

CALL  PL0T(18. 0,0. 0,-3) 

DO  615  K=1,N 
615  XG(K)=DEFL(INS,K) 

CALL  SCALE(X,10.0,N,l) 

CALL  SCALE(XG,10.0,N,1) 
C     PLOT  AXES  AND  ANNOTATION 

CALL  AXIS(0.0,0.0,12HTIME, SECONDS, -12, 10. 0,0.0,X(N+l),X(N+2),0) 

CALL  AXIS(0.0,0.0,17HDEFLECTION, INCHES, 17, 10.0, 90. 0,XG(N+l),XG(N+2 

*),-!) 
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CALL  SYMBOL(3.5,10.0,0.2,lUHSIGNATURE  PLOT, 0.0, l1*) 

CALL  SYMBOL(8.O,8.2,O.O8,l,0.0,-l) 

CALL  SYMBOL(8.1,8.lU,0.08,llH-RESPONSE  1,0.0,11) 

CALL  SYMBOL(8.O,8.1,0.08,2,O.O,-l) 

CALL  SYMB0L(8.1,8.0U,0.08,11H-RESP0NSE  2,0.0,11) 

CALL  SYMBOL(8. 0,8. 0,0. 08,3,0. 0,-1) 

CALL  SYMBOL (8. 1,7. 9^,0. 08, 11H-RESP0NSE  3,0.0,11) 

CALL  SYMB0L(8. 0,7. 9,0. 08, i*, 0.0,-1) 

CALL  SYMBOL(8.1,7.8i*,0.08,llH-RESPONSE  U, 0.0, 11) 

CALL  SYMBOL( 8. 0,7. 8, 0.08,5,0.0,-1) 

CALL  SYMB0L(8.1,7.71+,0.08,11H-RESP0NSE  5,0.0,11) 

CALL  SYMB0L(8.0,7.7,0.08,6,0.0,-1) 

CALL  SYMBOL  (8.1, 7. 6lt,0.08,llH-R£SPONSE   6,0.0,11) 

C  PLOT  SIGNATURE 

CALL  LINE(X,XG,N, 1,10,1) 
CALL  PLOT(l8.0,0.0,-3) 
1001  N=N1 
1000   CONTINUE 

IF(lPLOT.NE.O)    CALL  PLOT(0, 0,999) 
C  =============================================================== 

C 

100  F0RMAT( 313, F10. 5,^11) 

102  F0RMAT(8F10.5) 

103  FORMAT (12) 

10U   FORMAT(  '0====   FOURIER    COEFFICIENTS   FOR   G   FUNCTION   NO.  ',12,//, 
•'  A(I)    =  ',10F11.M 

105  FORMAT('0  B(I)    =',10E11.U) 

106  F0RMAT(9R8) 

107  FORMAT(1HO,1X////1X,9R8) 

108  FORMAT(5Fl6.12) 

200  FORMAT (1H1, IX,1    RECORD  OF  RAW  DATA    V/1X,*      TIME    RESPON 
*SE  1  RESPONSE  2   RESPONSE  3  RESPONSE  k     RESPONSE  5   RESPONSE  6 ' ) 

201  F0RMAT(12X,7F12.6) 

208  F0RMAT(1H1,1X/////////////////////////1X, ' 

*  DATA  SET  NUMBER' ,?X, 12) 

210  FORMAT(2l(8E10.3/),7E10.3) 

211  FORMAT(lHl, IX, 'RECORD  OF  RESPONSE   FUNS    7/1X, '  TIME  RESPON 
*SE   L     RESPONSE  2     RESPONSE   3     RESPONSE   U     RESPONSE   5      RESPONSE  6 ' ) 

250  FORMATCO        ILLEGAL  EOF  ENCOUNTERED') 
STOP 
50  PRINT  250 
STOP 
END 


SUBROUTINE   C0NVL1 
C  PURPOSE 

C  TO   COMPUTE  THE  VECTOR  OF  VALUES   OF 

XG   FROM   CONVOLUTION   EQUATION    XINP*XG=XOUT. 
FOR   GIVEN    XINP(I)    AND  XOUT(l) 
C  DISCRIPTION   OF  PARAMETERS 
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C  XINP-INPUT  VECTOR  OF  VALUES 

C  XOUT-OUTPUT  VECTOR  OF  VALUES  AFTER  CONVOLUTION 

C  XG-VECTOR  OF  VALUES  OF  CONVOLUTION  FUNCTION 

C  DT-INCREMENT  OF  THE  SUMMATION  POINTS 

C  N-NUMBER  OF  SUBINTERVAL  PLUS  1. 

C  SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  NEEDED 

C  NONE 

C  METHOD 

C  1.   SET  THE  VALUE  OF  'TEST'.  THEN  TEST  THE  VALUES  OF  XINP(l)  TO 

C  SEE  WHETHER  XINP(l)  IS  GREATER  THAN  TEST  OR  NOT 

C  2.   IF  MAGNITUDE  OF  XINP(K)  IS  GREATER  THAN  'TEST'  THEN  MAKE  UP 

C  K  MORE  VALUES  OF  XOUT(l)  AND  K  MORE  SUMMATION  EQUATIONS 

C  3.   CONSIDER  THE  SUMMATION  EQUATIONS  WHICH  IS  TRANSFERED  FROM 

C  INTEGRATION  EQUATION  OF  CONVOLUTION. 

C  k.      BY  FORWARD  SUBSTITUTION  METHOD  TO 

C  SOLVE  FOR  XG(I),I  =  1,..  .  ,N  . 

SUBROUTINE  CON VLI ( XINP, XOUT, XG, DT, N , PEAK) 
C 

C  

c 

DIMENSION  XINP(l),XOUT(l),XG(l) 

DIMENSION  Z(200),Y(200),W(200) 

A=0. 

TEST=. 000001 

DO  150  1=1, N 

150  XG(I)=0.0 

DO  151  K=1,N 

IF(ABS(XINP(K)).LT.TEST)  GO  TO  151 

L=K 

GO  TO  152 

151  CONTINUE 

152  IN=L-1 
IM=N+L-1 

DO  153  K=1,IN 

153  XOUT(N+K)=XOUT(N) 
DO  155  I=L,IM 
11=1-1 

DO  15k   K=L,II 
15U  A=A+XG(K-L+1)*XINP(I-K+L) 

XINP(L)=10.*PEAK 

XG(I-L+1)=(X0UT(I)-A*DT)/(XINP(L)*DT) 
155  A=0.0 

RETURN 

END 
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APPENDIX  D 

PREDICTING  PAVEMENT  PERFORMANCE 
COMPUTER  PROGRAM 


PROGRAM  CON VLE( INPUT, OUTPUT) 
C 

C  ///////////////////////////////////////////////////////////////// 

C 

C  PROGRAM  CONVLE 

C 

C 

C 

C  ...  PURPOSE 

C  READ  IN  DEFLECTION  RESPONSE  AND  TDT  FUNCTIONS  TO  PREDICT 

C  SIGNATURES  AND  TDT  FUNCTIONS.   THEN  USE  THE  PREDICTED 

C  SIGNATURE  AND  TDT  FUNCTIONS  TO  REPRODUCE  THE  DEFLECTION 

C  FUNCTIONS (  BY  C0NVL2) 
C  ...  DESCRIPTIONS  OF  VARIABLES 

C  XINP,XOUT,XG  -  VARIABLES  FOR  SUBROUTINE  CONVL2 

C  Y  -  THE  COORDINATE  ARRAY 

C  FUNS  -  INPUT  TDT  FUNCTIONS 

C  FIRST  INDEX  IS  THE  SITE  AND  AIRCRAFT  NUMBER 

C  SECOND  INDEX  IS  THE  TDT-X  FUNCTION 

C  THIRD  INDEX  IS  THE  ARRAY  NUMBER  IN  THE  TDT-X  FUNCTION 

C  DEFL  -  INPUT  DEFLECTION  FUNCTIONS 

C  INDICES  ARE  SAME  AS  IN  FUNS 

C  PFUN  -  PREDICTED  TDT  FUNCTIONS 

C  INDICES  ARE  SAME  AS  THE  SECOND  AND  THIRD  INDEX  OF  FUNS 

C  PDEFL  -  PREDICTED  DEFLECTION  FUNCTIONS 

C  INDICES  ARE  SAME  AS  IN  PFUN 

C  PTDT  -  THE  PEAK  VALUES  FOR  EVERY  TDT- 3  FUNCTIONS 

C  FCTOR  -  MULTIPLICATION  FACTORS  FOR  SIGNATURE  PREDICTION 

C  NT  -  TOTAL  NO.  OF  SITES  TIMES  TOTAL  NO.  OF  AIRCRAFTS 

C  N  -  THE  LENGTH  OF  EACH  INPUT  DEFL  AND  TDT  FUNCTIONS 

C  1ST  -  STANDARD  SITE  NUMBER 

C  LAIR  -  1=C135,2=C131,3=C130 

C  JA  -  STANDARD  AIRCRAFT  ON  STANDARD  SITE 

C  JJ  -  STANDARD  AIRCRAFT  ON  CURRENT  SITE 

C  JB  -  CURRENT  AIRCRAFT  ON  CURRENT  SITE 

C  JC  -  CURRENT  AIRCRAFT  ON  STANDARD  SITE 
C 

C  . . .  ROUTINES  AND  FUNCTIONS  NEEDED 

C  CONVL2 
C 

C  . . .  METHOD 

C  1.  USE  FACTOR  MULTIPLICATION  FOR  PREDICTING  SIGNATURE  (  BY 
FCTOR (M,N)) 
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C  2.  USE  FACTOR  MULTIPLICATION  FOR  PREDICTING  TDT  FUNCTIONS 

C  OTHER  THAN  TDT-3  (BY  FTR) 

C  3.  FOR  TDT-3  USE  SAME  METHOD  BUT  WITH  PEAK  AJUSTED  TO  THE 

C  CORRESPONDING  PEAK  OF  STANDARD  AIRCRAFT  ON  CURRENT  SITE 

C  k.    MAIN  LOOP  IS  ON  SITE  1  AIRCRAFT  123  (EXCEPT  IAIR) 

C  THEN  TO  SITE  2,3  (EXCEPT  STANDARD  SITE)   FOR  EXAMPLE, 

C  IST=2,IAIR=2,  THE  PREDICTIONS  AND  REPRODUCTIONS  WILL 

C  BE  ON  SITE  1  AIRCRAFT  1,3  AND  SITE  3  AIRCRAFT  1,3 

C 

C  

c 

DIMENSION  XINP(TO),XOUT(70),XG(TO),Y(5l),FUNS(9,6,5l)  ,PFUN(6,5l) 

DIMENSION  DEFL(9,6,5l),PDEFL(6,5l) 

DIMENSION  IAP(3),PTDT(9) 

DIMENSION  FCTOR(3,3) 

IAP(l)=l+RC135 

IAP(2)=URC131 

IAP(3)=1*RC130 

DO  91  1=1,3 

DO  91  J-1,3 
91  FCTOR(l,J)=l. 
C     FOLLOWING  FOUR  CARDS  ARE  REQUIRED  TO  INPUT  FACTORS  FOR 
C      CALCULATING  FOURTH  SIGNATURE 

FCTOR(l,l)*l.l»9 

FCTOR(l,3)=l.l*9 

FCT0R(3,1)=0.TT 

FCT0R(3,3)=0.TT 

NT=9 

N=U7 

DY=.0i+ 

IST=2 

IAIR=2 
C     INPUT  KNOWN  DEFLECTION  RESPONSE  FUNCTIONS  TDT  FUNCTIONS 

DO  10  1=1, NT 

DO  30  K=1,N 
30  READ  101,(DEFL(I,L,K),L=2,6) 

DO  35  K=1,N 
35  READ  101,(FUNS(I,L,K),L=2,6) 
10  CONTINUE 

DO  1*0  1=1, NT 

PTDT(I)=0. 

DO  1*1  K=1,N 

IF(ABS(PTDT(I)).LT.ABS(FUNS(I,3,K)))  PTDT( l)=FUNS(l ,3,K) 
1+1  CONTINUE 
UO  CONTINUE 
C     USE  FACTOR  MULTIPLICATION  TO  MODIFY  THE  PREDICTED  TDT  FUNCTIONS 

DO  50  1=1,3 

JA=(IST-1)»3+IAIR 

JJ=(I-1)*3+IAIR 
IF(I.EQ.IST)  GO  TO  50 
DO  51  J=l,3 
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JB=(I-1)»3+J 

JC=(IST-1)*3+J 

IF(J.EQ.IAIR)   GO  TO   51 

DO   81  M=3,6 

AJ=DJ=0. 

DO   80  K=1,N 

IF(ABS(AJ).LT.ABS(FUNS(JA,M,K)))   AJ=FUNS( JA,M,K) 

IF(ABS(DJ).LT.ABS(FUNS(JJ,M,K)))    DJ=FUNS( JJ,M,K) 

80  CONTINUE 
IF(AJ.EQ.O.)   AJ=1.E-13 
IF(M.NE.3)    GO  TO   16 

C  FOR  TDT-3   OF   CURRENT  AIRCRAFT  ON   CURRENT  SITE,    AJUST  THE  PEAK 

C  VALUE  TO  BE  SAME  AS  THE   CORRESPONDING  STANDARD  AIRCRAFT 

C  ON   CURRENT  SITE,   THEN   COPY  THE  CORRESPONDING   TDT-3  FROM  THE 

C  CURRENT  AIRCRAFT  ON   STANDARD  SITE  WITH  FACTOR  MULTIPLICATION 

C  TO  AJUST  THE  PEAK  VALUE  SPECIFIED  AS  ABOVE. 

AJ=PTDT(JC) 
DJ=PTDT(JJ) 
16  FTR=DJ/AJ 

DO  lb  K=1,N 
Ik  PFUN(M,K)=FUNS(JC,M,K)*FTR 

81  CONTINUE 

C    PREDICT  SIGNATURE  BY  RATIO  MULTIPLICATION 

DO  15  K=1,N 

PDEFL( 2 ,K )=DEFL( JC,2 ,K ) *FCTOR( I ,J ) 
15   CONTINUE 

DO   82  MY=1,5 

PRINT  216 

PRINT  21T,IAP(J),I 

PRINT  218 

DO  70   LA=1,N 

Y(LA)=FLOAT(LA-l)*DY 
TO  PRINT  201,Y(LA),(FUNS(JB,KA,LA),KA=2,5) 

PRINT  226 

PRINT  21T,IAP(J),I 

PRINT  218 

DO   71  LA=1,N 
71  PRINT  201,Y(LA),(PFUN(KA,LA),KA=2,5) 

82  CONTINUE 

DO   U5   11=3,6 
PEA=0. 
DO   1+6  K=1,N 
XINP(K)=PDEFL(2,K) 

IF(ABS(PEA).LT.ABS(XINP(K)))    PEA=XINP(K) 
1+6  XG(K)=PFUN(II,K) 
MK=0 
M=N 
C  REPRODUCTION  BY   CALLING  CONVL2 

CALL  C0NVL2(XINP,X0UT,XG,DY,M,PEA) 

DO   1+7  K=1,N 

IF(XOUT(K+MK).LT.O.)    XOUT(K+MK)=0. 
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1*7  pdefl(ii,k)=xout(k+mk) 
u5  continue 

DO  83  MY=1,5 
PRINT  213 

PRINT  2lU,IAP(j),I 
PRINT  215 
DO  76  NB=1,N 
76  PRINT  201,Y(NB),(DEFL(JB,KA,NB),KA=2,6) 
PRINT  223 

PRINT  2ll»,IAP(j),I 
PRINT  215 
DO  72  NB=1,N 
72  PRINT  201,Y(NB),(PDEFL(KA,NB),KA=2,6) 
83  CONTINUE 
51  CONTINUE 
50  CONTINUE 
STOP 
101  FORMAT(5Fl6.12) 
201  F0RMAT(13X,F11.2,5F11.6) 

213  F0RMAT(1H1,////36X, ^ORIGINAL  DEFL.  RESPONSE  FUNS?*) 
2lU  F0RMAT(1+0X,RU,?*  AIRCRAFT  SITE-?*, 11,//) 

215  FORMAT(20X,?*TIME     RESP  0     RESP  1     RESP  2     RESP  3      I 
*ESP  k#) 

216  FORMAT (1H1,////33X, ^ORIGINAL  TDT  FUNCTIONS?*) 

217  FORMAT(3UX,RU,?«  AIRCRAFT  SITE-?*, 11,//) 

218  FORMAT (20X,?*TIME     FUNC  1     FUNC  2     FUNC  3     FUNC  U?*) 
223  FORMAT  (lHl,////2l*X,  ^PREDICTED  DEFL.  RESPONSE  FUNS  (BY  RATIO  MULTIP 

*LI  CATION)?*) 
226  F0RMAT(1H1,////33X, ^PREDICTED  TDT  FUNCTIONS?*) 
END 


SUBROUTINE  C0NVL2(XINP,X0UT,XG,DT,N,PEAK) 
C 

C  ///////////////////////////////////////////////////////////////// 

C  SUBROUTINE  C0NVL2 

C  

C  PURPOSE 

C  TO  COMPUTE  THE  VECTOR  OF  VALUES   OF  XOUT(l)    FOR  GIVEN 

C  VECTORS... XINP( I)    AND  XG(l)    BY  USING  CONVOLUTION   TECHNIQUE. 

C  (I.E.    CONVOLUTE  XINP(l)   AND  XG(l)    BY  SUMMATION  METHOD) 

C  DISCRIPTION  OF  PARAMETERS 

C  XINP-INPUT  VECTOR  OF  VALUES 

C  XOUT-OUTPUT  VECTOR  OF  VALUES  AFTER   CONVOLUTION 

C  XG- VECTOR  OF  VALUES  OF  CONVOLUTION   FUNCTION 

C  DT-INCREMENT  OF  THE  SUMMATION  POINTS 

C  N-NUMBER  OF  SUB INTERVAL  PLUS   1. 

C 

C  SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  NEEDED 

C        NONE 

C  METHOD 
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C        TRANSFER  THE  INTEGRATION  FORMULA  OF  CONVOLUTION  INTO 

C        FINITE  SUMMATION  EQUATIONS  WITH  SUBINTERVALS  OF  LENGTH  DT  AND 

C        NUMBER  N-l.  COMPUTE  XOUT(l)  BY  GIVEN  VALUES  OF  XINP 

C        AND  XG  DIRECTLY 

C 

C  

c 

DIMENSION  XINP(l),XOUT(l),XG(l) 

XINP(l)=10.*PEAK 

DO  100  1=1 ,N 

X0UT(I)=0. 

DO  150  J=1,I 

X0UT(I)=X0UT(I)+XG(J)*XINP(I-J+1)*DT 

150  CONTINUE 
100  CONTINUE 

RETURN 

END 


129 


APPENDIX  E 
EXPLORATORY  TESTS 

A  series  of  exploratory  tests  were  conducted  on  an  operational 
taxiway  at  Kirtland  AFB,  New  Mexico  (Albuquerque  International 
Airport).   The  sensing  devices  used  to  measure  deflection  response 
functions  consisted  of  linear  variable  differential  transformers 
(LVDT)  and  accelerometers .   The  placement  of  these  sensing  devices 
relative  to  the  path  of  travel  of  a  prime  mover  are  as  shown  in 
Figure  52.   Figure  53  shows  the  instrumentation  setup  and  the  sensing 
devices  in  place  prior  to  a  data  run  with  a  prime  mover;  both  LVDT 
and  accelerometer  devices  are  depicted. 

The  prime  movers  selected  for  these  tests  included  a  load  cart 
with  a  25,000  pounds  single  wheel  load  and  a  Boeing  727  aircraft  with 
a  65,000  pounds  single  main  gear  load.   Figures  5^  and  55  show  the 
load  cart  and  aircraft,  respectively. 

Installation  and  Calibration 

Linear  Variable  Differential  Transformer  (LVDT).   The  differential 
transformer  was  installed  in  a  holder  which  could  be  removed  when 
not  in  use.   First,  a  six  inch  diameter  hole  was  drilled  through  the 
pavement  system  to  a  depth  of  approximately  16  feet  into  the  subgrade; 
a  four  inch  diameter  galvanized  steel  pipe  casing  was  then  inserted  in 
the  hole;  and,  subsequently,  a  3/8  inch  steel  reference  rod  was  placed 
inside  the  galvanized  steel  pipe  to  serve  as  a  reference  for  the  core 
of  the  differential  transformer  (see  Figure  56  for  a  detail 
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FIGURE    52.    LAYOUT   OF  SENSING   DEVICES 


FIGURE    53.     INSTRUMENTATION  SETUP 
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FIGURE   54.    LOAD   CART 


FIGURE  55.     BOEING  727  AIRCRAFT 
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representation  of  the  assembly).   The  holder  was  composed  of  an 
interior  threaded  aluminum  sleeve  permanently  fixed  to  the  pavement 
wearing  course  with  an  epoxy  (see  Figure  57).   The  differential 
transformer  coil  was  held  by  three  set  screws  in  a  brass  sleeve  which 
mated  with  the  aluminum  sleeve  fixed  to  the  pavement  surface.   A  cap 
with  an  extension  rod  internally  threaded  to  receive  a  LVDT  core  was 
fitted  over  the  reference  rod.   The  brass  sleeve  containing  the  LVDT 
coil  was  then  placed  in  position  within  the  aluminum  sleeve;  wherein, 
the  coil  moved  with  the  aluminum  sleeve  fixed  in  the  pavement  wearing 
course,  while  the  core  remained  in  a  fixed  position  on  top  of  the 
reference  rod.   Components  of  the  holder  assembly  are  shown  in 
Figures  58  and  59. 

Once  the  entire  assembly  was  completed,  calibrations  were  achieved 
by  turning  the  brass  sleeve  a  predetermined  number  of  revolutions  from 
a  null  position  within  the  fixed  aluminum  sleeve.   One  revolution  of 
the  brass  sleeve  provided  a  displacement  of  0.005  inch.   After  the 
calibrations  were  obtained  the  differential  transformer  was  placed  in 
a  null  position.   Data  was  placed  on  magnetic  tape  using  an  Ampex 
FR1300  tape  recorder  as  a  prime  mover  proceeded  along  the  taxiway. 

Accelerometer.   The  accelerometers  were  installed  on  a  removable 
pedestal  to  facililate  removal  of  the  sensing  devices  when  not  in  use 
(see  Figure  53).   The  pavement  surface  was  polished  slightly  with 
carborundum  paper  where  the  accelerometer  was  to  be  placed.   Subsequently, 
a  quick-drying  epoxy  was  used  to  fix  a  small  plug  to  the  pavement 
surface.   The  plug  was  internally  threaded  to  accept  the  external 
threads  on  the  accelerometer.   The  accelerometer  was  precalibratod  no 
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FIGURE   56.     LVDT   INSTALLATION    DETAIL 


FIGURE    57.      ALUMINUM   SLEEVE  (with  cop)  INSTALLED 
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FIGURE    58.    HOLOER   ASSEMBLY  COMPONENTS 


FIGURE    59.     HOLOER   ASSEMBLY  -  PARTI  ALLY  ASSEMBLED 
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that  it  could  be  turned  into  the  plug  and  placed  into  immediate 
operation.   Data  was  recorded  in  the  same  manner  as  described  above 
for  the  LVDT. 

Results  of  Exploratory  Tests 

A  comparison  of  peak  prime  rover  deflections  for  each  sensing 
device  is  summarized  in  Table  1  of  the  Field  Investigation  Section. 
It  was  concluded  that  accelerometers  were  generally  too  insensitive 
to  be  used  at  the  slower  speeds  used  in  the  tests.   Deflection  response 
functions  measured  with  LVDT  sensing  devices  were  good  when  the 
sensitivity  range  of  the  LVDT  was  correctly  selected. 

Figures  60,  6l,  and  62  depict  typical  deflection  response 
functions.   Attentuation  of  the  peak  deflections  were  used  as 
criteria  for  selecting  a  sensitivity  range  for  the  sensing  devices, 
and  additionally,  as  the  criteria  for  spacing  the  deflection  gages 
for  the  field  investigation. 
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APPENDIX  F 

DEFLECTION  GAGES  -  DESCRIPTION,  INSTALLATION,  AND  CALIBRATION 

DESCRIPTION 

The  deflection  gages  described  in  this  Appendix  were  specifically 
designed  and  fabricated  for  use  in  the  field  investigation  associated 
with  this  study  to  predict  pavement  performance  using  a  transfer 
function  approach.   Both  type  I  and  II  installations  employed  commercial 
displacement  transducers  (Models  2UDCDT-250,  2UDCDT-100,  and  2UDCDT-050) 
manufactured  by  Hewlett  Packard.   A  DCDT  coil  assembly  consists  of  a 
differential  transformer  coil,  a  DC-excited  solid-state  oscillator, 
a  phase-sensitive  demodulator,  all  in  one  small  package  0.75  in.  in 
diameter  with  variable  lengths  depending  on  the  particular  model.   The 
oscillator  converts  DC  input  power  to  AC,  which  is  used  to  excite  the 
primary  winding.   The  axial  core  position  determines  the  amount  of 
voltage  induced  in  the  secondary  windings.   Each  of  two  secondary 
circuits  contains  a  secondary  winding,  a  full-wave  bridge,  and  an  RC 
filter.   These  secondary  circuits  are  connected  in  series  opposition 
so  that  the  resultant  output  is  a  DC  voltage  linearly  proportional  to 
the  core  displacement  from  the  electrical  center.   The  polarity  of  the 
voltage  is  a  function  of  the  location  of  the  core  with  respect  to  the 
electrical  center,  which  is  located  approximately  half  way  along  the 
coil  length.   Input  and  output  circuits  are  electrically  isolated  from 
each  other  and  the  coil  assembly  case,  making  them  usable  directly  in 
floating  or  ground- return  systems.   The  DC  output  is  sufficient  to 
drive  most  standard  DC  indicators,  recorders,  or  control  systems.   The 


lUO 


design  of  these  transducers  eliminates  the  usual  phase-shift  correction 
and  harmonic-and-quadrature  null  problems  associated  with  externally 
excited  AC  differential  transformers. 

The  LVDT  elements  utilized  in  the  instrumentation  system  had 
linear  ranges  of  +0.250  in.,  +0.100  in.,  and  +0.050  in.  from  the  core 
center  for  the  2kDCDT-250,  2HDCDT-100  and  21+DCDT-050  models, 
respectively.   Each  instrument  provided  for  a  total  core  movement  of 
twice  the  linear  range  when  the  movement  in  both  the  positive  and 
negative  directions  were  considered. 

For  the  type  I  installation,  the  coil  assembly  was  mounted  in  a 
holder  which  was  attached  to  the  top  of  a  reference  rod.   The  core  was 
mated  to  a  vernier  holder  that  was  secured  to  a  mount  in  the  asphalt 
concrete  wearing  surface.   The  components  of  the  type  I  installation 
are  shown  in  Figure  63.   The  top  row  of  elements  consists  of  the 
components  which  remain  stationary  as  the  pavement  deflects  (viz., 
from  right  to  left  -  casing,  reference  rod,  casing  cap  with  teflon 
bushing,  coil  holder,  and  LVDT  coil).   The  bottom  row  of  elements 
consists  of  the  components  that  more  with  the  asphalt  concrete  surface 
(viz.,  from  right  to  left  -  LVDT  core,  vernier  holder,  vernier, 
pavement  mount ,  and  the  mount  cap ) . 

For  the  type  II  installation,  the  coil  assembly  was  mounted  in  a 
holder  which  was  attached  to  the  asphalt  concrete  surface.   The  core 
was  secured  in  a  beam  that  was  mounted  on  two  reference  rods.   The 
components  of  the  type  II  installation  (excluding  the  beam)  are  shown 
in  Figure  6U.   The  uppermost  and  lefthand  elements  consist  of  the 
components  which  remain  stationary  as  the  pavement  deflects  (viz., 
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FIGURE    63.      COMPONENTS   OF  DEFLECTION    GAGE 
TYPE  I 


FIGURE    64.      COMPONENTS  OF  DEFLECTION    GAGE 
TYPE  H 
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from  right  to  left  and  down  casing  and  cap  with  teflon  bushing, 
reference  rod,  connector,  beam  mount  extension,  beam  mount,  connector 
bolts,  LVDT  core,  and  the  beam  vernier).   The  elements  located  at  the 
bottom  and  to  the  right  of  the  figure  (i.e.,  the  LVDT  coil  holder  and 
LVDT  coil)  were  attached  to  the  asphalt  concrete  surface. 

Figures  65  and  66,   respectively,  show  partially  assembled  views 
of  both  type  installations.   Figures  67  and  68,  respectively,  show 
detailed  design  views  of  each  type  installation.   Figure  69  shows  a 
detailed  view  of  the  support  for  the  reference  beam. 

Both  Type  I  and  Type  II  deflection  gage  installations  were  employed 
at  each  test  site.   Normally  five  gages  were  employed;  however,  at  one 
site  six  gages  were  used  so  that  a  direct  measurement  of  prime  mover 
signatures  could  be  monitored.   A  schematic  representation  of  the 
transient  deflection  measuring  system  is  shown  in  Figure  TO.   Three 
regulated  power  supplies  were  used  to  drive  the  six  LVDT's.   The  time 
of  a  specific  data  run  was  recorded  by  use  of  an  IRIG  B  generator 
(Culton  Industries),  and  a  constant  DC  voltage  trigger  provided  event 
marks  on  the  magnetic  tape  when  a  prime  mover  passed  locations  2k   feet 
on  each  side  of  the  deflection  gages. 

INSTALLATION 

The  reference  rod  positions  were  located  on  18  inch  centers.  Holes 
six  inches  in  diameter  and  approximately  18  feet  deep  were  drilled  at 
each  reference  rod  location  (see  Figures  71  and  72).   Reference  rods 
were  then  cut  and  placed  in  each  hole.   At  each  location  a  reference 
rod  was  fitted  with  a  driving  head  and  driven  to  a  depth  of  six  to 
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FIGURE  65.    PARTIALLY  ASSEMBLED   DEFLECTION 
GAGE  TYPE  I 


FIGURE   66. 


PARTIALLY   ASSEMBLED    DEFLECTION 
GAGE   TYPE   II 
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FIGURE    67.    GAGE    INSTALLATION  -   TYPE    I 
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FIGURE  68.GA6E    INSTALLATION  -  TYPE   I 
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FIGURE    69.    REFERENCE    BEAM    SUPPORT 


1U7 


Power    Supply 
»   I 
Regulated 


Power    Supply 

■   2 

Regulated 


Time   Generator 
(IRIG  B) 


Event    Mark 


Power    Supply 

"    3 

Regulated 


5 


5- 


LVDT 

-  DC    « 


LVDT 
DC    «  2 


LVDT 

-  DC    ■  3 


LVDT 
DC  «   4 


LVDT 

DC   ■    5 


LVDT 
-   DC   ■    6 


-e  Channel     4 


-0  Channel     5 


-©Channel     3(H)S6(L) 


-©Channel     9(H)8I2(L) 


-©  Channel     7 


-©  Channel     8 


-©Channel     14 


-©Channel    I0(L)8I3(H) 


Tape    Recorder 
(H) and  (L)  Refer   to    High   and    Low    Sensitivity 

FIGURE  70.  SCHMATIC     REPRESENTATION     OF    THE     TRANSIENT 
DEFLECTION     MEASURING    SYSTEM. 


1U8 


FIGURE    71.     DRILLING  REFERENCE  ROD  HOLES 
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FIGURE    72.    HOLES   READY  FOR   REFERENCE  RODS 
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nine  inches  below  the  bottom  of  the  drilled  hole.   The  hole  was  then 
back  filled  and  tamped  with  approximately  12  inches  of  subgrade 
material.   A  three  inch  diameter  galvanized  pipe  casing  was  cut  to 
appropriate  length  and  placed  over  and  around  each  reference  rod.   A 
casing  cap  with  a  teflon  bushing  was  placed  over  the  protruding 
portion  of  the  reference  rod  and  turned  into  place  on  the  top  of  the 
pipe  casing.   Each  reference  rod  was  fixed  at  the  bottom  of  the  drilled 
hole  by  virtue  of  driving  it  six  to  nine  inches  and  then  placing 
approximately  12  inches  of  tamped  backfill.   The  casing  cap  fitted 
with  a  teflon  bushing  provided  alignment  and  lateral  support  for  the 
reference  rod  near  the  top  of  a  drilled  hole.   Subsequently,  loose 
subgrade  material  was  placed  in  the  space  between  the  pipe  casing  and 
the  wall  of  each  drilled  hole  to  an  elevation  approximately  one  inch 
below  the  casing  cap. 

To  provide  for  lateral  alignment  of  the  pavement  mount  with 
respect  to  the  reference  rod,  a  bushing  (Figure  73)  was  placed  over 
each  reference  rod,  and  inserted  into  the  pavement  mount.   A  styrofoam 
sleeve  (Figure  7M  was  then  inserted  over  and  around  the  bushing 
between  the  casing  cap  and  the  pavement  mount.   This  sleeve  prevented 
epoxy  (when  poured  to  fix  pavement  mount  to  wearing  surface)  from 
getting  near  the  reference  rod  and  casing  cap  components.   Later,  this 
sleeve  served  to  keep  the  paper  filler  from  interferring  with  placement 
of  a  LVDT  coil,  or  a  reference  rod  extension.   Paper  filler  was  then 
placed  around  the  pavement  mount  and  styrofoam  sieeve  to  approximately 
three  inches  below  the  pavement  surface  (Figure  75)-   This  paper  filler 
served  as  a  compressible  medium  between  reference  components  and 
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FIGURE    73.     BUSHING   FOR   LATERAL  ALIGNMENT 


FIGURE    74.     STYROFOAM    SLEEVE    IN  -  PLACE 
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movable  components  of  the  deflection  gages.   Vertical  alignment  between 
the  pavement  surface  and  the  pavement  mount  was  achieved  by  fixing  a 
horizontal  steel  plate  to  the  mount  cap  (Figure  76).   The  plate 
in- turn  rested  on  the  horizontal  pavement  surface,  thereby,  providing 
vertical  alignment  between  the  pavement  surface  and  the  mount  cap. 
Figures  77  and  78  show  the  unassembled  and  assembled  components, 
respectively,  and  how  both  lateral  and  vertical  alignment  of  the  pave- 
ment mount  with  respect  to  the  reference  rod  was  achieved.   Figures  79 
and  80  depict  views  of  the  mounts  installed  subsequent  to  pouring  and 
curing  the  epoxy. 

An  epoxy  (Hysol,  Resin  R8-2038  and  Hardener  H2-3U75)  and  sand 
mixture  was  poured  in  the  volume  between  the  pavement  mount  paper 
filler  and  the  inside  of  the  drilled  hole.   Herman  Nelson,  Type  H-l 
1+00,000  BTU/Hr  heaters,  were  used  to  cure  the  epoxy  for  a  period  of 
two  hours  at  a  temperature  of  approximately  lUo  F.   Figures  81  and 
82  show  the  epoxy  curing  operation. 

Type  I  Installation  -  The  major  objective  of  this  type 
installation  was  to  achieve  measurements  of  the  deflection  response 
functions  while  at  the  same  time  protect  the  sensing  device  from  the 
tire  of  a  prime  mover.   After  a  pavement  mount  was  securely  fixed  to 
the  pavement's  wearing  course  with  the  epoxy  mixture,  a  LVDT  coil 
holder  and  LVDT  coil  was  placed  on  the  top  threaded  end  of  a 
reference  rod.   Then,  a  LVDT  core,  core  vernier,  and  vernier  holder 
assembly  was  placed  inside  the  pavement  mount.   This  assembly  was 
adjusted  to  an  approximate  null  position  inside  the  pavement  mount  and 
securely  fixed  in  place  by  tightening  the  three  set  screws  in  the 


152 


FIGURE   75.       PAPER    FILLER 


FIGURE    76.     PLATE   FOR   VERTICAL   ALIGNMENT 


153 


FIGURE    77     UNASSEMBLED    ALIGNMENT  COMPONENTS 


FIGURE    78.    ASSEMBLED    ALIGNMENT  COMPONENTS 
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FIGURE    79.      MOUNT   INSTALLATION,  CLOSE-UP 


FIGURE    80.      MOUNT   INSTALLATION, DISTANT 
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FIGURE   81.        EPOXY   CURING,  CLOSE -UP 


FIGURE    82.     EPOXY    CURING,  DISTANT 
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vernier  holder.   The  deflection  gage  was  then  nulled,  calibrated, 
and  the  mount  cap  installed. 

Type  II  Installation  -  The  reference  beam  was  first  put  in  place 
by  placing  a  connector  and  beam  mount  extension  on  top  of  each  of  two 
reference  rods.   Lateral  support  for  the  beam  mount  extension  was 
provided  by  a  mount  cap  fitted  with  a  teflon  bushing  which  was  placed 
over  the  extension  rod  and  secured  in  the  top  of  the  pavement  mount. 
The  beam  was  then  secured  to  the  beam  mount  extensions  an '3  adjusted  to 
a  predetermined  height  of  approximately  three  inches  above  the  pavement 
surface.   At  each  gage  position  on  the  beam  a  LVDT  coil  and  holder 
assembly  was  aligned  with  a  LVDT  core  and  vernier  in  the  beam,  and  the 
holder  glued  to  the  pavement  surface  with  Hysol  608  Epoxy-Patch. 

CALIBRATION 

Each  gage  was  nulled  for  zero  voltage  output  prior  to  calibration. 
The  gages  were  nulled  by  turning  the  core  vernier  until  zero  voltage 
output  was  achieved.   The  gage  output  was  monitored  by  an  oscillscope. 
Once  a  null  position  had  been  obtained  for  each  gage,  the  tape  recorder 
was  started  and  a  zero  reading  was  taken  for  each  gage.   Subsequently, 
the  LVDT  core  for  each  deflection  gage  was  displaced  a  known  amount 
by  turning  the  core  vernier  a  predetermined  rotation.   With  20  threads 
per  inch  on  the  core  vernier,  this  rotation  permitted  calculation  of 
the  actual  displacement  of  the  core.   Figure  83  shows  the  calibration 
devices  used  to  achieve  a  desired  rotation  of  the  core  vernier.   Figure 
8U  represents  actual  manner  in  which  calibrations  were  Derformed.   After 
each  core  had  been  displaced,  the  tape  recorder  was  started  and  a  step 
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FIGURE    83.      CALIBRATION  OEVICES 


FIGURE    84.      METHOD   OF  CALIBRATION 
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calibration  reading  taken  for  each  gage.   This  calibration  procedure 
was  repeated  a  minimum  of  two  times  daily.   Occasionally  a  third 
calibration  step  was  taken  to  check  linearity  of  the  instruments. 
All  linearity  checks  proved  the  gages  be  within  two  percent  error. 
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APPENDIX  G 

DEFLECTION  RESPONSE  AND  TDT  FUNCTIONS 

Data  for  the  deflection  response  functions  and  TDT  functions  are 
contained  in  this  Appendix.   The  deflection  response  function  data 
are  listed  first  for  all  the  aircraft  at  each  test  site.   These  are 
followed  by  the  calculated  TDT  functions.   All  data  was  printed  out 
to  six  decimal  places  from  computer  code  CONVLI. 
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,000041 

.000000 

1.72 

.0^9785 

.001716 

0.000000 

•000040 

.000000 

1.7b 

.009738 

.001712 

0.000000 

•000040 

.000000 

1.8o 

.000747 

.00164R 

0.000000 

,000039 

,000000 

1.8* 

.009540 

.001653 

0.000000 

.000039 

.000000 

DEFLECTION  RESPONSE  FUNCTIONS 
C131  AIRCRAFT-SITE  2 
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TIME 

RESP  0 

RESP  1 

RESP  2 

RESP  3 

HESP  4 

O.Ou 

.003280 

.004297 

,000994 

.000030 

.000000 

.04 

.00^989 

.004784 

.001059 

.000037 

.000000 

.08 

.004639 

.005016 

.001067 

.000041 

.000000 

.12 

.005497 

.005697 

.001126 

.000041 

.000000 

.16 

.006602 

.006461 

.001204 

.000031 

.000000 

.20 

.007653 

.006869 

.001204 

.000036 

.000000 

.24 

.009083 

.007991 

.001276 

.000030 

.000000 

.20 

.010561 

.008448 

.001278 

.000030 

.000000 

.32 

.01?642 

.009350 

.001354 

.000033 

,000000 

.3b 

.0U642 

.010086 

.001391 

.000030 

,000000 

.40 

.016610 

.010419 

,001399 

.000038 

.000000 

.44 

.018532 

.010800 

.001418 

.000045 

.000000 

.48 

.0205l3 

.011532 

,001465 

.000039 

.000000 

.52 

.021914 

.011607 

.0014P7 

.000042 

.000000 

.5b 

.02J727 

.011625 

,001*75 

.000048 

.000000 

.60 

.023000 

.011496 

.001501 

.000036 

.000000 

.64 

.022424 

,010850 

.001437 

.000039 

,000000 

.68 

.021247 

.010640 

.0014(17 

.000042 

,000000 

.72 

.020207 

.009805 

.001365 

.000041 

.000000 

.76 

.018406 

.008780 

.001308 

.000037 

.000000 

.8y 

.016699 

.008077 

.001306 

.000044 

.000000 

.84 

.015102 

.007459 

.001277 

.000040 

.000000 

.88 

.013362 

.006542 

.001182 

.000042 

.000000 

.92 

.011647 

.006052 

,001115 

.000046 

,000000 

.96 

.016663 

,005672 

,001128 

.000035 

.000000 

1.00 

.009174 

.004701 

,001052 

,000046 

.000000 

1.04 

.008065 

.004404 

.001003 

.000046 

.000000 

1.08 

.007048 

,00379f, 

.000972 

.000040 

.000000 

1.12 

.006118 

.003123 

.000891 

.000040 

.000000 

1.1b 

.005295 

.002812 

.000879 

.000039 

.000000 

1.20 

.004483 

.002483 

.000805 

.000039 

.000000 

1,24 

.004033 

,002364 

.000783 

.000040 

.000000 

1.28 

.003728 

.002190 

.000759 

.000043 

.000000 

1.32 

.003298 

.002159 

,000739 

.000044 

.000000 

1,3b 

.002857 

.002056 

,000709 

.000045 

.000000 

1.40 

.002718 

.002046 

,000656 

.000046 

.000000 

1.44 

,002478 

,001632 

,000622 

.000050 

.000000 

1.48 

.002081 

.001589 

,000568 

.000048 

.000000 

1.52 

.002148 

,001389 

.000619 

.000041 

.000000 

1.56 

.002069 

,001417 

.000567 

,000037 

.000000 

1.60 

.001814 

.001271 

.000517 

.000043 

.000000 

1.64 

.001763 

.001315 

.000510 

.000035 

.000000 

1.68 

.001593 

,001187 

.000459 

.000035 

.000000 

1.72 

.001377 

,000922 

.000404 

.000042 

.000000 

1,76 

.001253 

,000981 

.000349 

.000040 

.000000 

1.8U 

,001232 

.000725 

.000305 

.000040 

.000000 

1.84 

,000855 

,000704 

.000226 

.000040 

.000000 

DEFLECTION  RESPONSE  FUNCTIONS 
C130  AIRCRAFT-SITE  2 
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TIME 

0.00 

•  04 

.Ob 

.12 

.16 
.20 
.2* 
.28 
.32 
.36 
.40 
.44 
.48 
.52 
.56 
.60 
.64 
.68 
.72 
.76 
.80 
.84 
.88 
.92 
.96 
1.00 
1.04 
1.08 
1.12 
1.16 

1.20 
1.24 
1.28 
1.32 
1.36 
1.40 
1.44 
1.48 
1.52 
1.56 
1.60 
1.64 
1.68 
1.72 
1.76 

1.80 
1.84 


RESP  0 
.002115 
.002931 
.003922 

.005522 
.00P379 
.01?903 
.019336 
.025234 
.02oifl0 
.02P168 
.03o380 
.033319 
.037199 
.041952 
.0*7268 
.051000 
.056380 
.046558 

.04J770 
.037249 
.033524 

•  03()3o6 
.02PU2 
.02*767 
.02*958 
.02*928 
,02*878 
.02R515 
.030506 
.03?5l5 
.034685 
.03*779 

.036222 
.035237 
•033384 
.031056 
.028383 
.02*558 
•02?970 
•02o6lO 
,018520 
.016797 
•015309 

•  0U064 
.013054 
.01?267 
.011449 


RESP  1 
.002658 
.003185 
.004165 
.005624 
.007690 
.010317 
.013530 
.015946 
.017315 
.018223 
.018765 
.019939 
.021372 
,023120 
.024888 
.025778 
,025294 
,02388c 
.021884 
.019761 
.018059 
.016477 
.015319 
.014715 
.014218 
,013997 
.014237 
.014660 
.015034 
.015547 
.015854 
.015921 
.015869 
.015328 
.014415 
.013330 
.012074 
.010565 
.009147 
.007764 
.006553 
.005267 
.004205 
.003540 
.002816 
.002518 
.002096 


RESP  2 
.000829 
.000578 
,001129 
.001331 
.001609 
.001900 
.002169 
,002415 
.002607 
.002708 
.002761 
.002871 
,002995 
.003116 
.003246 
.003329 
,003307 
.003232 
.003065 
.002864 
.002704 
.0025*9 
.002390 
.002300 
.002224 
.002138 
.002110 
.002115 
.002088 
.002077 
.002063 
.002010 
.001992 
.001924 
.001858 
.001767 
.001636 
.001477 
.001320 
.001171 
.001031 
,000858 
.000767 
.000709 
•000600 
.000488 
.000360 


RESP  3 
.000004 
.000004 
.000010 
.000017 
.000024 
.000032 
.000041 
.000044 
.000044 
.000044 
.000046 
.000048 
.000070 
.000083 
.000079 
.000083 
.000087 
.000094 
.000105 
•000108 
.000109 

.000109 
.000110 
•00010* 
.000107 
.000105 
.000103 
.000099 
.000097 
.00009* 

.000095 
.000093 
.000090 
.000090 
.000089 
•00008P 
.000087 
•000086 
•000084 
•000084 
.000083 
.000083 
•00008] 
•000081 
.000080 
.000079 
.000079 


RESP  4 
•000000 
.000000 
.000000 
.000001 
.000001 
.000001 
.000001 
.000001 
.000001 
,000001 
.000001 
.000001 
.000001 
.000001 
.000001 
.000001 
.000001 
.000001 
.000001 
.000001 
.000001 
.000001 
.000001 

•oooooi 

.000001 
.000001 
.000001 
.000001 
.000001 
.000001 
.000001 
.000001 
.000001 
.000001 
.000001 
.000001 
.000001 
.000001 
.000001 
.000001 
.000001 

.oooooi 

.000000 
.000000 
.(J00000 
.000000 
.000000 
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TIME 

o.oo 

.04 
.08 

•  12 
•16 

•20 
.24 
.28 
.32 
.36 
.40 

•  44 
.48 
.52 
.56 
.60 
.64 
.68 
.72 
.76 
.80 

•  64 
.88 
.92 
.96 

1.00 
1.04 
1.08 
1*12 
1.16 
1.20 
1.24 
1.28 
1.32 
1.36 
1.40 
1.44 
1  .48 
1.52 
1.56 
1*60 
1*64 
1.66 
1.72 
1.76 
1.80 
1.84 


RESP  0 
.00*164 
.003392 
.005238 
•008189 
.012830 
.019777 
.030280 
.042851 
.048505 

.044063 
.038199 
.034642 
.034125 
.037696 
.045787 
.052000 
.047899 
.038441 
.029318 
.021388 
.015337 
.010951 
.007993 
.006115 
.004934 
.004214 
.003875 
.003453 
.003204 
.003008 
.001702 
.001474 
.001469 

.001437 
.001399 
.001399 
•001380 
•001385 
•001362 
.001327 
.001324 
.001324 
.001324 
.001307 
.001298 
.001319 
•001287 


DEFLECTION 
C135  AI 
RESP  1 
.^02195 
.004583 
.008671 
•012884 
•019370 
•026787 
.033371 
.036015 
•  03520<? 
.033683 
.03050? 
.030546 
.03241* 
.035401 
.035994 
.034324 
.029419 
.021910 
.014971 
.009963 
.006363 
.003989 
.002089 
.001543 
•000816 
•000709 
.000760 
.000251 
.000307 
.000089 
.000301 
.000467 
.000290 
•000008 
•000071 
•000202 
•000168 
•000143 
.000167 
•000034 
.000007 

.000299 
•000331 
•000070 
•000373 
•000434 
.000007 


PF^PONSE  HJNCT 

RCRAFT-SHE  3 
HfSP  2 
.000677 
.001166 
.002196 
•003536 
.005182 
.007129 
.009263 
.010894 
.011637 
.011692 
.01163' 
.011517 
.011540 
.011625 
.01188b 
.01196b 
.01 1630 
.010765 
.008831 
.006589 
.004647 
.003062 
.00193a 
.001226 
.000775 
.000489 
.000413 
•000295 

•000221 
•000166 
.000147 
•00014' 
.000147 
.000102 
,000139 
.000132 
•000125 
.000055 

.000009 
.000009 
.000009 
.000055 
.000009 
.000009 
.000009 
.000009 
.000009 


IONS 

HESP  3 
000217 
000279 
000365 
000517 
000724 
00097P 
001299 
001589 
001739 
00179? 
001841 

001856 
001868 
001895 
001907 
001903 
001854 
001801 
001666 
001525 
00133? 
001121 
00090» 
00074? 
00061? 
000507 
000446 
000392 
000357 
000320 
000235 
000201 
000181 
000171 
000163 
000158 
000142 
000132 
000125 
000113 
000112 
000110 
000101 
000096 
000095 
00008P 
OOOOBR 


PESP  4 
.0  00  052 
.000063 
.000075 
.000087 
.000009 
.000112 
.000124 
.000134 
.000139 
.000142 
.000144 

.000145 

e000U6 

.000146 
.000144 
.000139 
.000131 
.0001?? 
.000109 
♦0000^5 
,0000*? 
.000069 
.000058 
.00004^ 
.00004? 
•00003S 
.000031 
.0000?7 
•0000?4 

.0000?1 
.000014 
.000013 
.000011 
,000010 
.000009 
,000008 
.000007 
.000007 
.000006 
.000005 
.000005 
.000005 
•000004 
.000004 
.000003 
.000003 
.000003 
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DEFLECTION  RESPONSE  FUNCTIONS 

C131  AIRCRAFT-SHE  3 

TIME 

RESP  0 

RESP  1 

«ESP  2 

RESP  3 

RESP  4 

0.00 

•001195 

^01558 

000446 

000162 

.000025 

.04 

•001461 

.001550 

000566 

000180 

.000028 

.08 

.001792 

002407 

000725 

000200 

.000030 

.12 

.002238 

003087 

00094b     , 

000233 

•000033 

•  16 

•002823     « 

00394? 

001201 

000260 

•000036 

.20 

.003491 

004301 

001*59 

000290 

•000038 

•  24 

.004309 

005571 

001736 

000324 

•000040 

.28 

•005291 

.006185     « 

002017     , 

.000373 

•000042 

.32 

.006451 

007801 

00230b 

000396 

.000043 

•  36 

.007829     , 

.009196 

002630     < 

.000443 

•000045 

.40 

.009553 

010373 

002982 

000499 

.000046 

•  44 

.011503 

011718 

003325 

.000537 

•000048 

•  48 

.013670 

.012158 

.003624 

.000581 

•000049 

•  52 

.015776 

.012467 

003878     « 

.000625 

•000050 

•  56 

.017313     « 

.012385 

.004048 

.000642 

•000050 

•  60 

.018000 

.012215 

.004142 

.000660 

.000050 

•  64 

.017654 

.011883 

.004134 

.000654 

•000050 

•  68 

•016476     < 

.011069 

.004087     < 

.000639 

•000049 

•  72 

.015069     < 

.010040 

.004007 

.000637 

•000048 

.76 

.013709 

.009384 

,00387/ 

.000625 

•000047 

.80 

.012344 

.008097 

,003679     , 

.000613 

.000046 

.84 

.011069 

.007526 

.003*38 

.000613 

•000045 

•  88 

.009885 

.006574 

.003138 

.000589 

•000043 

•  92 

•008768 

.005543 

.002830 

.000568 

•000041 

•  96 

•007727 

.004986 

,002527     , 

.000540 

•000039 

1*00 

•006809     , 

.004283 

,002222     i 

.000514 

•000037 

1.04 

•006014 

.003939     < 

,001927 

.000479 

.000035 

1.08 

•005285 

.003437 

,001661 

.000455 

.000033 

1*12 

•004664 

.003013 

,001*18 

.000429 

.000030 

1.16 

•004149     ■ 

.002563 

.001197     , 

.000395 

.000028 

1.20 

.003694 

.002313 

.001031 

.000363 

.000025 

1.24 

•003326     i 

.002248 

,000874 

.000334 

•000023 

1.28 

.003103 

.002185 

►000771 

.000317 

.000021 

1*32 

•002856 

.001809 

.000686     , 

.000297 

.000020 

1.36 

•002663 

.001877 

.000612 

.000271 

•000019 

1.40 

.002558     , 

•001686 

.000558 

•000269 

.000018 

1.44 

•002440     , 

.001964 

.000494 

.0002*6 

•000016 

]  .46 

.002337 

•001473     . 

»ooo446   ; 

.000232 

.000015 

1.52 

.002253 

.001544 

•000404     , 

•000220 

•000014 

1*56 

•002158 

•001521 

•000373     » 

•000204 

.000013 

1.60 

•002074 

•001431 

•000345 

.000190 

.000012 

1*64 

•002039 

.001491 

•000316 

•000172 

•000011 

1*66 

•002001 

•001397 

•000301 

•000170 

•000010 

1.72 

.001954 

•001360 

.000275 

•000153 

.000009 

1.76 

•001905 

•001279 

.000263 

•000144 

•000009 

1*80 

.001914 

.001383 

.000259 

•000149 

•000009 

1.84 

.001879 

.001073 

.000240 

•000133 

•000008 
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DEFLECTION  RESPONSE  FUNCTIONS 
C130  AlHCRAFT-SITE  3 


TIME 

RESP  0 

RESP  1 

RESP  2 

RESP  3 

WESP  4 

0.00 

.00?228 

.003533 

.000815 

.00027^ 

,00^049 

.04 

.002836 

.003076 

.000972 

.000296 

.000053 

.0« 

.001548 

.005718 

.001376 

.000336 

.00C057 

.12 

.004614 

.005735 

.001852 

•000399 

.000062 

.16 

.00«5878 

.008751 

.002389 

.000454 

.000067 

.20 

.007304 

.010322 

.002905 

.000511 

.000072 

.24 

.009256 

.012749 

.003500 

.000567 

.000077 

.28 

.011813 

.014582 

.004292 

.000677 

.000082 

.32 

.014704 

.017741 

.005034 

.00078? 

.0  00086 

.36 

.01A237 

.019563 

.005719 

.00088? 

.000091 

.40 

.022295 

.022471 

.006508 

.00101* 

.000094 

.44 

.026881 

.024479 

.007135 

.001111 

.00009b 

.48 

.031484 

.025993 

.007819 

.001205 

.000101 

.52 

.035346 

.026324 

.008258 

.001273 

.000104 

.56 

.037938 

.027011 

.008651 

.001348 

.000)06 

.60 

.039000 

.026126 

.008853 

.001390 

.000108 

.64 

.038570 

.026100 

.008974 

.001421 

.000109 

.68 

.03f,8]9 

.024567 

.008915 

.001430 

.000  109 

.72 

.034487 

.024153 

.008766 

.001414 

.000109 

.76 

.032112 

.022546 

,008649 

.001428 

.000109 

.80 

.029962 

.023248 

.008574 

.001422 

.000108 

.84 

.02*207 

.021018 

.008358 

•001400 

.OOOlOM 

.88 

.02^844 

.021171 

.008222 

.001399 

.000106 

.92 

.026115 

.020098 

.008102 

.001401 

.000106 

.96 

.025926 

.021402 

.007924 

.00139S 

.000104 

1.00 

.026282 

.020089 

.007771 

.001364 

.000102 

1.04 

.027480 

.021715 

.007714 

.001361 

.000100 

1.08 

.029202 

.020858 

.007673 

.001347 

.000098 

1.12 

.036870 

.021358 

.007639 

.001335 

.000096 

1.16 

.031681 

.020317 

.007468 

.001323 

.000093 

1.20 

.031470 

.019285 

.007258 

.001289 

.000090 

1.24 

.03o222 

.017145 

.006965 

.001271 

.000088 

1.28 

.02R258 

.015557 

.006610 

.001222 

.000084 

1.32 

.026040 

.013821 

.006270 

.001187 

.000081 

1.36 

.023557 

.012865 

.005736 

.001149 

.000077 

1.40 

.02o926 

.010308 

.005231 

.001089 

.000073 

1.44 

.018598 

.009020 

.004718 

.001045 

.000069 

1.48 

.016316 

.006737 

.004127 

.000978 

.000065 

1.52 

.0U295 

.005755 

.003529 

.000908 

.000060 

1.56 

.012388 

.004153 

.002971 

.000821 

.000056 

1.60 

.010708 

.003911 

.002454 

.000745 

.000052 

1.64 

.009315 

.002072 

.001916 

.000684 

.000049 

1.68 

.008090 

.002313 

.001588 

.000623 

.000045 

1.72 

.007019 

.000878 

.001209 

.000553 

.000042 

1.76 

.006107 

.001376 

.000929 

.000513 

.000039 

1.80 

.005334 

.000529 

.000631 

.000466 

.000036 

1.84 

.004734 

.000437 

.000448 

.000424 

.000033 
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TOT  FUNCTIONS 

Cl 

35  AIRCRAFT- 

■SITE  1 

TIME 

FUwC  1 

FUNC  ? 

FUNC  3 

FUNC  4 

.00 

•OoqOOO 

.0023*8 

.000010 

.000000 

'  .0* 

.1U1932 

.003173 

.000011 

.000000 

•  OH 

.1^*460 

.0040-37 

.000013 

.000000 

.12 

.16?ig2 

.004907 

.0000  15 

•ooouoo 

•  16 

.2«i7275 

. 0063^4 

.000016 

.000000 

.20 

.30^122 

.007H77 

.000017 

. ooouoo 

.24 

.4JPH37 

.009264 

.000019 

.000000 

.28 

.5J6481 

.0108*7 

.000019 

.000000 

.32 

.61^318 

.0124*5 

.000020 

.ooouoo 

.3b 

.6*s742 

.013259 

.000020 

.ooouoo 

•  4U 

.5o«594 

.0130^1 

.000019 

.000000 

.44 

.49j7]4 

.012357 

.000018 

.000000 

.48 

.4*1401 

.011630 

.000017 

.oouooo 

.52 

.463034 

. 0110*5 

.000015 

.000000 

•  5b 

.518528 

.010616 

.000013 

•ooouoo 

.60 

#48c;399 

.0096*6 

.000011 

.000000 

.64 

.33Q006 

.0079/2 

.000009 

.000000 

.68 

.159824 

.005617 

.000006 

.000000 

.72 

-.011402 

.003002 

.000003 

.000000 

.76 

-.Ib39q5 

.000267 

•  0000c  1 

.000000 

.80 

-.2o-i972 

-.002008 

-.000001 

.000000 

.84 

-.3J7872 

-.003763 

-.000003 

-.000000 

.88 

-.3**405 

-.005104 

-.000004 

-.000000 

•  H<L 

-.4in527 

-.006024 

-.000005 

-.000000 

.96 

-.410740 

-.006417 

-.000006 

-.000000 

1.00 

-.391693 

-.006576 

-.000006 

-.ooouoo 

1.04 

-.336395 

-.006340 

-.O000C6 

-.000000 

1.08 

-.27*573 

-.005654 

-.000007 

-.000000 

1.12 

-.220256 

-.005065 

-.00000b 

-.000000 

1.16 

-.1^0225 

-.00445? 

-.000005 

-.ooouoo 

1.20 

-.U4951 

-,0035b4 

-.300004 

-.000000 

1.24 

-.0^8527 

-.002616 

-.000004 

-.000000 

1.28 

-.021219 

-.001907 

-.000002 

-.ooouoo 

1.32 

.0*1126 

-.000939 

-.000001 

-.ooouoo 

1.36 

.0^3351 

.00011* 

-.000001 

-.ooouoo 

1.40 

.1H371 

.000713 

.000001 

-.000000 

1.44 

.140301 

,0012*H 

.000001 

.000000 

1.48 

.147629 

.001708 

.00C002 

.000000 

1.52 

.1*9*50 

.002102 

.000002 

.000000 

1.56 

.141349 

.00194? 

.000002 

.000000 

1.60 

.123972 

.002333 

.000002 

.ooouoo 

l.b4 

.0"*Q50l 

.001844 

.000001 

.000000 

1.68 

.0^0172 

.001454 

.000002 

.000000 

1.72 

.04*531 

.0011*5 

.000001 

.000000 

1.76 

.021628 

.000917 

-.000000 

.000000 

1.80 

-.009373 

.000216 

.000000 

-.000000 

1.84 

-.029268 

-.C00207 

-.000001 

-.OOOUOO 
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ri 
F  UNC    1 

.'ill  AQfl 
. 0  1 os it 
.012147 
.01610? 
.024081 
.046159 
.087519 
.1532S* 

.33380^ 
,42  7Hc;? 

.51641A 

.639195 

.64900A 
.621*20 
.ts;«740n 
.^,43529 
,48256? 

. TQMr Ifl 

.?O9901 
. 1*S457 
. 039767 
-.  571733 
■.  1  *S7U 
.,246060 
-.31  B^OQ 
.,  ^g?ci 

■.43396P 

■. 434*07 
«.4?957l 

• .  4  0  7  o  0  0 

■.  17?JH 

,,3?460^ 

••  2*7174 

■.202977 

■.  1  3420* 

..062720 

,0049ft4 

.  168*79 

.  12733* 

. 1 7*789 

.?1S?54 

.24318? 

.262325 


TOT    FUNCTI 
31     AIRCRAFT 
FlJNC    2 
.0042*17 

.0  04l<n5 
.005416 

.005^92 

,no6^?4 

.007T9P 
.007*1* 
.00*653 
,009Hrt9 
.010499 
.01 14S7 
.01172* 
.011-01 

.oiimss 

.011/471 
.010*06 
.0O9HS3 

.00RH75 
.007*19 
.00*'+«6 
.00^*03 
,004^fl7 
,003443 
,0024?fl 
.001 39] 
.000^94 
-.000^4 

-.001201 

-.001***. 

-.002127 
-.002397 

-.002*30 
-,00?r>44 
-,00]73* 
-.001911 
-.001^5 
-.000^25 

-.000221 

.000315 
.000^7 
.0014*9 
.001*?^ 
.00?740 
.003109 
,003?97 

.003536 
.00^198 


0NS 
-SITE  1 

FUNC  3 

o.ooooou 

0.000000 

OvOOtrpou 
o.ooooou 

.000006 
.000014 

.oooou 

.000009 
-.000000 
-.000000 
-  .  0  0  0  0  0  0 
-.000001 
-.000001 
-.000001 
^".ODOOOl 
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.OOOOOb 
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.oooou 
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TOUOOOS 
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-.OOOOOJ 
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TT0WI7 

.000001 
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.000011= 

.000001 
-.UOOOOJ 

.U00021 


FUNC  4 

.00  00  '0 
.00O0"O 
.000000 

.000000 
.000000 
.0000  0 o 
•OOOOuO 
.000000 
.000000 

.oooooo 

.000000 

.oooooo 

.000000 
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.oooooo 

.OOOOOO 
.000000 
.000000 

.oooooo 
.oooooo 
.oooooo 
.oooooo 
•oooooo 
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.oooooo 
.oooooo 

•OOOOOO 
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.000000 
.OOOOOO 
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.oooooo 

.000000 
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.OOOOOO 
•OOOOOO 
.000000 

.oooooo 

.000000 
.000000 
•OOOOOO 

•oooooo 
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.2' 

.342216 

00«793 

.  o  o  o  o  l  v 

,00  00'  (! 

.  *  '' 

• 3«3751 

0  0973ft 

,  00  0  01*? 

.000000 

.  ■*'• 

,  44">3  7ft 

,  f!  J  |V->Q7 

.  o oooi  y 

, OUOOuU 

,  >i 

,49333ft 

,  0  1  0 '-•  ft  4 

.  1)0001  * 

.00 00 oo 

.  <*' 

.^^i^gQ 

0  1  1  *  4  5 

,  0  0  0  0 1  a 

.00  00 'iO 

,  u  < 

.S9f.lrt"» 

,on  -v7u- 

,000  0  3  a 

.  OOOOi'O 

.  VP 

.61743] 

nl2  L?l 

.  o  0  0  0  1  1 

.0000  /(/ 

.  s»- 

,ft03?57 

011574 

.OOnoie 

.oono  '0 

.  »■» 

.<^74'<iS?               , 

0  1  0  9  o  3 

.i»00014 

•  00  OOti  tj 

.  *  '■ 

.4974  03 

,  n  1  n    ^  p 

.0000  1  J 

.0000  >0 

.  <1 

,394223 

,  0  0  9  2  ft  ft 

.000011 

,000000 

.  / " 

.287    in 

,  0  0  « ■"«  0  H 

.  0  0  0  0  0  v 

,0000     (i 

.  /- 

.  1  «  1     1  1 

, 0 0  7 «   7  1 

.  o  o  o  o  o  a 

,00001  () 

. K " 

.  0  7  6  1 6  3 

,r  06^ci5 

.000007 

.000000 

.  H" 

-  •  0  0  5  ^  1  ft 

,0  05534 

.  o  0  0  0  0  fc 

.  flu 0  0 »»o 

.  H> 

-.057220 

,  0  0  4  -  7  (i 

.000005 

.030000 

•  M  ' 

-  ,  I  *  *<->  4  g             , 

,  n  o  i  r «  7 

.oooou* 

,0000  10 

m<il 

-.162275 

,002  y  1  ft 

.00000 J 

.00  00  10 

1  .  n 

-.192273 

,00?'l  7? 

.O000OJ 

.0  000  ■;(') 

I      '^  *« 

-,20ft527 

0  0 1 '.  23 

.00000^ 

.oooooo 

1  •  " ' 

-,??^2?7 

OOO^ftl 

.  o  0  0  0  0  d 

,00 00 dO 

1  .  1  '' 

-.232725 

00Q334 

. u  0 0  0  01 

.01  00'  0 

1.1 

-.23573R 

00012ft 

,0  00000 

.  OoOO   10 

I  .2 

-,25352ft 

non-7  / 

-,000000 

.  0  >"/  0  0  f  o 

'  .<!  • 

-,?5?420 

0  011  ft  w 

-.  0O0001 

,0000  .il 

1  .  <" 

-.267616 

.ouhq^ 

-.000001 

-  .  n  o oo  i  o 

1  .  3  ' 

-,?Q7444 

0  02oi<j 

-.000001 

-  .  00  OO'H) 

1  .  i' 

-.33423^ 

,0u?7(S7 

-.  0  000 0<! 

-  .  0  100   iO 

1  .<♦ 

-.35n^l« 

,0  03297 

-.00000* 

-.00  0  0  10 

1  .4* 

-.355185 

,  n  o iHHn 

-.00000*! 

-.n  ioo  mi 

1  .v- 

-.327207 

,  0  0  4  204 

-.00000J 

-.QuPO  10 

1.  .  5  •' 

-,300ft37 

.niH'it;? 

-  .  0  0  0  0  0  J 

-•00 0000 

1  .  V 

-,?7^767 

,0  04*50 

-  .  0  i)  0  0  0  J 

-.OUOO'.MJ 

1  ,f?  ' 

-.253936 

,00534  3 

-  .  0  0  0  0  0  J 

-.0 0  00O0 

1  ,fv' 

-•230N93 

,  0  0  5 '-  <■>  ? 

-  .  0  0  0  0  0  J 

-.0,100' Ml 

1   .ft' 

-•20603ft 

,  005t»02 

-.ooooo j 

-,0000  i  0 

I  ..    1? 

-,177Mcji 

,  00519ft 

-. OO000J 

-.0000 00 

1    .    '' 

-  .  1 4  5  3  M 

,  0  0  ^  n  4  7 

-•00000<£ 

-,0  0  00  if) 

1    .  M 

-.10  77?3 

,O0.(  /<)Q 

-.00000* 

-  .0   100   Ml 

1    .  M'. 

-.065164 

,  0  0  4  "5  1  4 

-.uonoui 

-.0   IOO     0 

172 


TOT  FUNCTIONS 

Cl35  AIKCRAF-T-SITE  2 

TIME 

FUNC  1 

FUNC  2 

FUNC  3 

FUNC  4 

0.00 

•105370 

.036808 

.001013 

.000019 

.0* 

.I3?lfi5 

.043918 

.001340 

.000022 

.Ob 

.181399 

.053357 

.001859 

.000025 

.12 

.24fl802 

.064031 

.002412 

.00002H 

•  lb 

.35?503 

.077733 

.003005 

.000032 

•  20 

.514487 

.095359 

.003596 

.000034 

•  24 

. 7*8338 

.114171 

.004124 

.000037 

•  2fl 

1.024636 

.133632 

.004583 

.000038 

.32 

1.206726 

.147222 

.004833 

.000038 

.36 

1.204223 

.146529 

.004782 

.000037 

.40 

1.0/7595 

.135344 

.004567 

.000035 

.44 

.9JP273 

.122089 

.004291 

.00003? 

.48 

.8b4399 

.110341 

.003981 

.000030 

.b2 

.839253 

.099298 

.003649 

•000027 

.56 

.84q752 

.089923 

♦003290 

•000024 

.60 

.767997 

.075790 

.002887 

.000020 

.64 

.5/3377 

.052412 

.002501 

.000016 

.68 

.2S»78l3 

.024006 

.002149 

.00001? 

.72 

.0J5972 

-.002166 

.001829 

.000009 

.7b 

-.189261 

-.026286 

.001557 

•000006 

.80 

-.369197 

-.046932 

.001303 

•000003 

.Bh 

-.51^327 

-.064227 

.001050 

.000001 

.88 

-.62?941 

-.075612 

.000711 

-.000002 

,92 

-.68?475 

-.081635 

.000344 

-•000004 

.96 

-.6*5831 

-.082672 

.000040 

-•000005 

1.00 

-.657133 

-.078878 

-.000250 

-.000007 

1.C4 

-.583231 

-.071714 

-.000435 

-.000007 

1.06 

-.497624 

-,062468 

-.000565 

-.000008 

1.12 

-.3*9708 

-.052087 

-.000624 

-.000008 

1.16 

-.303856 

-.040533 

-.000663 

-.000008 

1.2H 

-.202529 

-.029179 

-.000598 

-.000007 

1.24 

-.101271 

-.017461 

-.000491 

-.000006 

1.28 

-.005141 

-.006497 

-,000355 

-.000004 

1.32 

.09*557 

.003616 

-.000217 

-.000003 

1.36 

.1^6735 

.011881 

-.000065 

-.000002 

1.40 

.229653 

.018932 

.000075 

-.000001 

1.44 

.267963 

.024297 

.000219 

-.000000 

1.48 

.295408 

.027361 

.000331 

.000000 

1.52 

.2*7090 

.028177 

.000462 

.000001 

1.56 

.2«1614 

.026897 

.000550 

.000001 

1.60 

.253336 

,023801 

.000596 

.000001 

1.64 

.203872 

.019282 

.000628 

.000001 

1.68 

.153768 

.013779 

.000635 

.000001 

1.72 

.0*7461 

,007710 

.000622 

.000001 

1.76 

.045959 

.001476 

.000585 

.000000 

1.80 

-.015170 

-.004559 

.000540 

-.000000 

1.84 

-.O63836 

-.010072 

.000470 

-.000001 
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TUT    FUNCTIONS 
ClJl    AIRCRAFT-SITE    2 


TIME 

FUmC  1 

FUNC  ? 

FUNC  3 

FUNC  4 

0.00 

.4OQ510 

•108575 

.003264 

.000039 

.0* 

.514572 

.113856 

.004033 

.000040 

,Q6 

.520570 

•112398 

.004352 

.000040 

.12 

.59,558 

.116177 

.004230 

.00003* 

.16 

.659043 

•121376 

.003023 

.000037 

.20 

.681867 

,11751ft 

.003506 

.000037 

•  24 

.78?5fl8 

.120786 

.002650 

.000035 

.26 

•80?897 

.115559 

.002496 

.000034 

.32 

.86^9! 3 

.117513 

.002654 

.000032 

.36 

.90^690 

.114193 

•002162 

.000031 

•  40 

.890664 

.106706 

.002773 

.000028 

•  44 

.876625 

.099694 

.003213 

•000026 

•  4<3 

.894077 

.094687 

.002251 

•000022 

.52 

.833639 

.086427 

.002376 

•OOOOIQ 

.56 

.76^037 

.074226 

.002684 

.000018 

.60 

.67^032 

.066426 

.001106 

.000013 

.64 

.52*624 

.049279 

.001160 

•000011 

•  68 

.434085 

.036967 

.001225 

•000009 

.72 

.274367 

.024391 

.000926 

•000007 

.  76 

.101128 

.011828 

.0003(36 

.000006 

•  80 

-.02flU3 

.006771 

.001009 

.000005 

.84 

-.13fl347 

.000450 

.000478 

.000003 

.88 

-.269735 

-.011526 

.000766 

.000003 

.92 

-.3409o8 

-.018420 

.001176 

.000004 

.96 

-•3«7383 

-.015274 

-.000026 

.000003 

1.00 

-.4«4042 

-.0202*4 

.001236 

.000005 

1.04 

-.4^1586 

-.020788 

.001364 

.000006 

1.08 

-.518817 

-.018117 

.000788 

•000006 

1.12 

-.53H6?7 

-.019791 

.000965 

•000009 

1.16 

-.505236 

-.012868 

.000949 

•OOOOOQ 

1.20 

-.4b?206 

-.012092 

.001104 

.000010 

1.24 

-.38*824 

-.005152 

.001392 

.000012 

1.28 

-.31]357 

.001557 

.001828 

.000014 

1.32 

-.216065 

.008575 

.002054 

.000015 

1.36 

-.127155 

.014280 

.002246 

.000016 

1.40 

-.030161 

.016896 

.002475 

.000017 

1.44 

.018021 

.020909 

.002910 

.000017 

1.48 

.101124 

.024168 

.002796 

.000017 

1.52 

.157359 

.033539 

.002003 

.00001* 

1.56 

.227742 

.032511 

.001537 

.000014 

1.60 

.267284 

.030958 

.002197 

.000014 

1.64 

.314182 

.032972 

.001281 

•000012 

1.68 

.328654 

.029062 

.001 188 

•000011 

1.72 

.3H970 

.024040 

.001866 

.000011 

1.76 

.323873 

.018143 

.001562 

.000009 

1.80 

•28s937 

.012746 

.001475 

.000007 

1.84 

.26-^898 

.003124 

.001387 

.000006 

TDT    FUNCTIONS 
CliO    AlRCHAFT-SITE 
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TIME 

FUnC  1 

0.00 

.13H97 

.0^ 

.15*475 

.OB 

.201638 

.12 

.273flo3 

.16 

.37?5?2 

.20 

.49^840 

.2* 

.64f,737 

.28 

.753210 

.32 

.803862 

•  36 

.827386 

•  40 

.827562 

.44 

.85^771 

.48 

.883976 

.52 

.924042 

.56 

.9b«862 

•  60 

.944969 

.64 

.85*970 

.60 

.723282 

.72 

.555481 

.76 

.378969 

.80 

.22?696 

.84 

.0^5234 

.88 

-.044476 

.92 

-.126661 

.96 

-.191576 

1.00 

-.230856 

1.04 

-.236363 

1.0« 

-.223165 

1.12 

-.203483 

1.1b 

-.168301 

1.20 

-.134953 

1.24 

-.105510 

1.28 

-.075423 

1.32 

-.065051 

1.30 

-.071253 

1.40 

-.086657 

1.44 

-.113285 

1.48 

-.156589 

1.52 

-.200034 

1.56 

-.246300 

1.60 

-,287fii7 

1.64 

-.335619 

1.68 

-.373257 

1.72 

-.390379 

1.76 

-.407322 

1.80 

-.398448 

1.84 

-.388855 

FUNC  2 
.040900 
.048017 
.055143 
.064553 
.077414 
.090389 
.101515 
.110728 
.116558 

.117337 

.115047 

.114808 

.114361 

.112914 

.110968 

.105994 

.0954^0 

.082312 

.064752 

.045757 

.029356 

.013973 

-.000532 

-.010323 

-.017990 

-.024698 

-.027279 

-.027156 

-.027694 

-.026601 

-, 024952 
-.024577 
-.021964 
-.021314 
-.020264 
-.020339 
-.022443 

-.025946 
-.029617 
..033029 
-.036176 
-.041029 
-.041771 

-, 040712 
-.041822 
-.042654 
-.043719 


FUNC  3 
.000187 
.000195 
.000514 
.000856 
. 001191 
.001558 
.001975 
.002077 
.002024 

.0019P5 
.001935 
.002010 
,002982 
.003532 
.003195 

.003215 

.003248 

.003405 

.003718 

.003618 

.003394 

.003138 

.002947 

.002627 

.002346 

.002014 

.001666 

.001273 

.000975 

.000724 

.000547 

.000322 

.000051 

-.000061 

-.000172 

-.000230 

-.000284 

-.000298 

-.000333 

-.000287 

-.000227 

-.000155 

-.000121 

-.000035 

.000037 

.000080 
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FUNC  4 
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.000021 
.000025 
.000026 
.000027 
.00002H 
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.00002* 
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.000022 
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-.000007 
-.000006 
-.000006 
-.000006 
-.000006 
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TIME 


.00 

.04 

.06 

•  12 

•  16 

•  20 

1. 

•  24 

1, 

•  28 

1. 

.32 

1. 

•  36 

1. 

•  40 

It 

•  44 

1. 

•  48 

1. 

•  52 

1. 

•  56 

•  60 

•  64 

•  68 

•  72 

-• 

.76 

m     g 

.80 

-. 

•  84 

-1. 

.88 

-1. 

•  92 

-1. 

.96 

-# 

•  00 

-« 

.04 

•  9 

.08 

-• 

•  12 

-. 

•  16 

-, 

•  20 

-. 

•  24 

•  28 

•  32 

.36 

•  40 

•  44 

•  48 

•  52 

.56 

.60 

•  64 

•  68 

•  72 

•  76 

-• 

.80 

«  m 

•  84 

-, 

TD1 

FUNCTIONS 

C135    « 

»IRCRAFT- 

•SITE    3 

FUNC    1 

FUNC    2 

FUNC    3 

FUNC    4 

106026 

0*2703 

.010495 

002532 

220644 

056214 

.013416 

003014 

416254 

105345 

.018386 

00356P 

615563 

16098? 

.024544 

004104 

921077 

246380 

.034131 

004623 

265266 

336632 

.045737 

005097 

559334 

433292 

.060131 

005511 
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501743 

.072370 

005675 

554869 

522560 

.077216 

005561 

408451 
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.076808 

005331 

164607            . 

476384 

.075574 

005060 

061529 
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•07200J 

004730 

035923 

405143 

.067552 

•004360 

061348 

370943 

.063226 
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978078 

344606 

.057730            . 

-003427 

798934 

311249 

.051508            < 

.002750 

477765 

261086 
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.035993 
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361224 
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.025936 

000330 
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.002334 
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.331659 
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625011         -< 
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473515 
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-.031643         -< 
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-.026107         -, 
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.067613 

-.016349 

.000551 
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.01638? 
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260471 

.030863 

-.004118            « 

►000249 

355870 
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.001988 
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444015            . 
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506487 
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►163812 

.015903 

►001221 

517889 

►171894 

.018753 

►001286 

477292 

►170059 
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.001244 

402108 

►159750 

.020832 

.001160 

312138 

.140054 

•0?034V 

.001026 

226184 

•115233 

.018590 

•000814 

124669 
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•015557 
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•047499 
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•000338 

061044 
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134085 
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21535? 

.042973 

.001226 

.000309 
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TIME 

0.00 

.04 

.08 

•  12 

•  16 
.20 
.24 
.28 
.32 
.36 
.40 

•  44 

•  48 
•52 

•  56 
.60 

•  64 
.68 
.72 
.76 
.80 
.84 
.88 
.92 
.96 

1*00 
1.04 
1.08 
1.12 
1*16 
1*20 
1«24 
1.28 
1.32 
1.36 
1.40 
1.44 
1  .48 
1.52 
1.56 
1.60 
1*64 
1.68 
1«72 
1.76 
1.80 
1*84 


CI 
FUNC  1 
.216880 
.213904 
.331101 
.422087 
.535797 
.578351 
.745524 

.817773 
.025755 
.197577 
.332550 
•483479 
.499369 
.488029 
.412907 
•316314 
.189179 

.987923 
.751380 
.562961 
.284441 
•105792 
•123207 
.356459 
•512436 
.673632 
.765130 
.855902 
.911084 
•943805 
•923022 
.852006 
.758864 
•690479 
.545252 
.425989 
.232740 
.151657 
.00690? 
.144015 
.259411 
.379972 
.460549 
.528500 
.568290 
.611541 
.574833 


TDT 
31  A 


FUNCTI 
IRCRAFT 
FUNC  ? 
062023 
07B310 
099671 
129089 
163186 
196770 
232171 
267147 
301B00 
340228 
380399 
417103 
445119 
463974 
467954 
458484 
432206 
398219 
357674 
308637 

249202 
183583 

110063 
036794 
033436 
10067? 
162069 
21408ft 
256627 
289005 
306380 
314865 
308338 
292388 
268360 
235829 
200228 
158963 
11506Q 
068999 
023565 
019252 
060293 
094961 
125600 
150737 
166791 


ONiS 

-SITE  3 

FUNC  3 
.0*2520 
.0?4861 
.02742* 
.03162V 

.034937 
,03859b 

.042606 

•04R401 

.050443 

.05557* 

.0M427 

.06442V 

.067754 

.070491 

•069121 

.067375 

.062080 

.055347 

.05035V 

.043901 

.037530 

.032913 

.025226 

.018136 

.01044V 

.003613 

-.003941 

-.009116 

-.013695 

-.018624 

-.022295 

-.024534 

-.024431 

-.023906 

-.u?355* 

-.019279 

-•017493 

-•014071 

-•010147 

-.00661b 

-.002966 

.000061 

.004855 

.007243 

•01013V 

.014436 

.015071 


FUNC  4 

003491 

003809 

004100 

004443 

004799 

005058 

005226 

005335 

005332 

005366 

005359 

005311 

005067 

004836 

004432 

003991 

003434 

00283? 

002308 

001746 

001258 

000773 

000324 

000080 

000483 

000801 

001105 

001358 

001546 

001674 

001736 

001705 

00153* 

001328 

001088 

000812 

000598 

000364 

000132 

000117 

000345 
000559 
000764 
000855 
000968 
001049 
001058 
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TDT  FUNCTIONS 

130 

1  AIRCRAFT-SITE  3 

TIME 

FUnC  1 

KUNC  2 

FUNC  3 

FUNC  4 

^.00 

.22*576 

.0522*3 

,017648 

.003150 

.04 

.19^646 

.061991 

.018826 

.003405 

.OH 

.363269 

.087354 

.021223 

.003634 

.12 

.3t>0690 

.116993 

.025062 
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.16 

.54Q608 

.150029 

.028268 

.004140 

.20 

.643210 

.18)216 

.031476 
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.24 

•78Q190 

.216678 

.034479 
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.28 
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.040776 
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.306433 

.046705 
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1.H3H9 

.343*67 

.051655 

.005063 

.40 

1.32P521 

.386025 

.058547 

.005067 

•  44 

1.419298 

.415506 

.062709 

.005032 

.48 

1.469353 

.446029 

.066264 

.004897 

.52 

1.433967 

.457960 

.067667 

.004713 

.56 

1.41]776 

.463907 

.069163 

•00443S 

.60 

1.278962 

.454567 

.068137 

.00410S 

.64 

1.19?543 

.437104 

.066017 

.003757 

.68 

1.001230 

.405362 

.062235 

.003336 

.72 

.87036 

.3653*6 

.056589 

.002909 

.76 

.667616 

.325487 

.052665 

.002478 

.80 

.60^314 

.286727 

.047304 

.002074 

.84 

.353646 

.23790P 

.040804 

.001664 

.88 

.257626 

.193915 

.035643 

.00130? 

.92 

.087755 

.151495 

.030655 

.000965 

.96 

.0^8745 

.106594 

.025316 

.000601 

1.00 

-.0H7934 

.065425 

.018569 

.000250 

1.04 

-.0t>2655 

.033244 

.013902 

-.000055 

1.08 

-.I63UI 

,005417 

.008937 

-.000322 

1.12 

-.1  M169 

-.018189 

.004585 

-.000601 

1.16 

-.263300 

-.046432 

.000675 

-.000880 

1.20 

-.3*0024 

-.072823 

-.004004 

-.001133 

1.24 

-.474405 

-.100120 

-.007015 

-.001322 

1.28 

-.56H34 

-.127140 

-.011413 

-.001525 

1.32 

-.647215 

-.149139 

-.014295 

-.001684 

1.36 

-.674663 

-.179862 

-.016748 

-.001841 

1.40 

-.79fl254 

-.205366 

-.020147 

-.001966 

1.44 

-.834278 

-.228251 

-.021997 

-.002040 

1.48 

-.92R803 

-.253293 

-.024867 

-.002073 

1.52 

-.934120 

-.275995 

-.027498 

-.002109 

1.56 

-.97?642 

-.293385 

-.030778 

-.002065 

1*60 

-.91^709 

-,305316 

-.032965 

-.001970 

1.64 

-.951827 

-.315688 

-.033640 

-.001824 

1.68 

-.84?999 

-.3092*1 

-.033936 

-.001666 

1.72 

-.830091 

-.302829 

-.034349 

-.001467 

1.76 

-.680888 

-.286540 

-.032367 

-.001220 

1.80 

-.6U7399 

-.267914 

-.030419 

-.000964 

1.84 

-.4/6253 

-.238721 

-.027727 

-.000682 
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APPENDIX  II 

PREDICTED  DEFLECTION  RESPONSE  FUNCTIONS 

The  predicted  deflection  response  functions  are  contained  in  this 

Appendix.   The  functions  were  calculated  by  using  computer  code  CONVLE. 

All  nine  combinations  of  aircraft  and  test  sites  were  employed,  and 

similar  results  as  shown  in  the  typical  plots  (Figures  h~J   and  U8)  were 

achieved.   Computer  printouts  of  the  predicted  dynamic  deflections  for 

four  of  the  nine  combinations  are  provided.   These  results,  labeled 

cases  1  through  U,    are  representative  of  all  the  combinations.   They 

include: 

Case  1.   C131  Aircraft  and  Site  2  -  Standard 
Predictions  provided  for  C135  and 
C130  aircraft  at  both  sites  1  and  3. 

Case  2.   C135  Aircraft  and  Site  2  -  Standard 
Predictions  provided  for  C131  and 
C130  aircraft  at  both  sites  1  and  3. 

Case  3.   C130  Aircraft  and  Site  2  -  Standard 
Predictions  provided  for  C135  and 
C131  aircraft  at  both  sites  1  and  3. 

Case  h.      C131  Aircraft  and  Site  1  -  Standard 
Predictions  provided  for  C135  and 
C130  aircraft  at  both  sites  2  and  3. 
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Case  1 

C131  Aircraft  and  Site  2  -  Standard 
Predictions  for: 


C135  Aircraft  at  Site  1 
C130  Aircraft  at  Site  1 
C135  Aircraft  at  Site  3 
C130  Aircraft  at  Site  3 
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PREDICTED     DEFL 
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Case  2 

C135  Aircraft  and  Site  2  -   Standard 
Predictions  for: 


C131  Aircraft  at  Site  1 
C130  Aircraft  at  Site  1 
C131  Aircraft  at  Site  3 
C130  Aircraft  at  Site  3 
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TYPICAL   CALCULATED  AND   MEASURED 
DEFLECTION  RESPONSE  FUNCTIONS  FOR  CI3I 
USING  PREDICTED  TDT  FUNCTIONS, WITH  CI35 
AS  STANDARD  PRIME  MOVER  . 
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FIGURE  86.  TYPICAL  CALCULATED  AND  MEASURED  DEFLECTION 
RESPONSE  FUNCTIONS  FOR  CI30  USING  PREDICTED 
TDT  FUNCTIONS  WITH  CI35  AS  STANDARD 
PRIME  MOVER. 
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PPEDICTPD  DFF| 
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HATIO  MULTIPlTCATlOtvi 

C13A  AlRfRAFT  SITE 
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TIME 

RESP  0 

PESP  1 
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.004599 
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.00691  1 
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Cl3l  AIRCRAFT  SITE-3 


TIME 

PESP  0 

RESP  1 

RFSP  2 

RESP  3 

RFSP  4 

'"•OO 

.002525 
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PREDICTED  DEFL.  RESPONSE  FUNS  (BY  HATIO  MULT  IP|_  TCATIOM) 

Cl30  AIRCRAFT  SITF-3 


TIME 

RESP  0 

RESP  1 

fcrSP  2 

RESP  3 

RFSP  4 

n.00 

.001628 

.003*44 

.002280 

000047 

.000044 

.04 

.0022*7 

.004247 

.002690 

00004O 

.000048 

.08 

.003020 

.00*552 

•0O3107 

000130 

•000054 

.1? 

•0042*2 

.007408 

.003662 

000216 

•000O5O 

.16 

•0064*2 

.0102*1 

,004427     , 

000302 

.0000*3 

.20 

.00993? 

•013754 

.00*227 

00039P 

.000066 

•  ?4 

•014889 

•018036 

.00*966 

000508 

.000070 

.?a 

.019430 

.021254 

.00*643     « 

000541 

.000073 

.32 

.021699 

.023076 

•0 071 69 

00053O 

.00007* 

.3^ 

.022459 

.0242*3 

.007447 

000544 

.000077 

.40 

.023393 

.02*000 

.007590 

000563 

•000OB0 

•  44 

.025656 

.02**60 

.007891     . 

.000593 

.0000*4 

•  48 

•028643 

•028463 

.008229 

000863 

.000001 

•  52 

•032303 

.030785 

.008560 

.001032 

.000006 

•  *6 

•03639ft 

•033132 

.008913 

.000983 

.000097 

•  AO 

•039270 

•034308 

.009138 

.001027 

•OOOO09 

•  64 

•03«793 

•033652 

.0O9Q75 

>00107fi 

•000101 

•  68 

•03585c 

•031754 

.008864     « 

.001166 

•000103 

•  7? 

.032163 

.029079 

.008402 

.001301 

.000106 

.76 

•02*68? 

.026237 

.007844     . 

.001338 

.000107 

•  «0 

•0258]4 

.023953 

•0O7402 

.00134* 

.00010* 

.84 

.023336 

.021832 

.006973     < 

.001344 

.00010* 

,P8 

.021670 

.020276 

.006533 

.001356 

.00010* 

.92 

.0206] 1 

.019460 

.006283     < 

•001334 

.000102 

.06 

•019987 

.018790 

.Oo*07l 

•001322 

.000101 

i  .00 

•01006s 

•018490 

.00*833     « 

•00129* 

.000009 

1  .04 

•02069ft 

•0188Q8 

.00*756 

.001264 

•0000°6 

1  .08 

•021957 

•019370 

.00*769 

.001220 

•0000O2 

1  .12 

.023489 

.019868 

.00*695     . 

.001195 

.000090 

1.16 

.025037 

•020554 

.005666 

.001177 

.000088 

i  .20 

.026707 

.020968 

.005629     , 

.001172 

•0000P6 

1.24 

.027550 

.021062 

•0o*485     . 

•00114B 

.000084 

l  .?R 

.027891 

.020998 

.005435 

.001109 

.000081 

i  .32 

.027133 

.020284 

•00*250 

.001102 

.0000*0 

i  .36 

.025706 

.019074 

.005071 

.001080 

.000070 

1  *4Q 

.023913 

.017634 

.004822 

•001082 

.000077 

l  .44 

.021855 

•01*966 

.004461 

.001069 

•00007* 

l  .4? 

.019680 

.013959 

.004028 

.001059 

•00007> 

i  .52 

.017687 

.012075 

.003596 

.001037 

.0000*9 

1  .56 

.015870 

•010235 

.003188 

•001030 

.00006* 

1  .*o 

.014?60 

•0086?3 

•Oo?803 

•001023 

.000062 

1  .64 

.012934 

.006911 

.002328 

.001016 

.000060 

i  .68 

•011788 

.005499 

.002080 

•000997 

•0000*6 

1  .72 

•0108?9 

•004615 

.001921 

.000991 

.000054 

1  .76 

.010051 

.003656 

.001624 

.000982 

.000051 

1  .PO 

.009445 

.003260 

.001318 

.000967 

.000048 

1  .R4 

.008816 

.002703 

.000967 

.00096P 

.000044 
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Case  3 

C130  Aircraft  and  Site  2   -  Standard 
Predictions  for: 


C135  Aircraft  at  Site  1 
C131  Aircraft  at  Site  1 
C135  Aircraft  at  Site  3 
C131  Aircraft  at  Site  3 
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FIGURE    87.     TYPICAL  CALCULATED  ANO  MEASURED  DEFLECTION 
RESPONSE   FUNCTIONS    FOR   CI35  USING 
PREDICTED    TDT   FUNCTIONS   WITH  CI30    AS 
STANDARD    PRIME    MOVER. 
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FIGURE   88.     TYPICAL  CALCULATED  AND  MEASURED 
DEFLECTION  RESPONSE   FUNCTIONS  FOR 
C  131  USING  PREDICTED  TDT  FUNCTIONS 
WITH   CI30  AS   STANDARD   PRIME  MOVER 
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.024516 

. 010152 

.002134 

.000163 

.72 

.030590 

.01  J2  0  0 

. 00J542 

.002173 

.  000163 

.  76 

.u24047 

.0l456o 

. o 0  64  56 

. 002161 

.  000157 

•  is  u 

.01 *640 

.011067 

. 005551 

. 002139 

.000150 

.  j4 

.016210 

. 003430 

. 004313 

. 002095 

.  000142 

.  do 

.013923 

.  006655 

. 00  3437 

. 001993 

.  000132 

.  -j  2 

.012473 

.0-0  5**  9  6 

. 0027o6 

.001852 

.000121 

.  j6 

.  Jll4bl 

.  0o<+5  7  b 

.002264 

.001710 

.000109 

1 .  u  0 

.010320 

. 004023 

. 001857 

. 001549 

.  000097 

1.  04 

.  01032b 

.003660 

. 001510 

. 001411 

.  000086 

1  .  UO 

.009391 

.0o3243 

.001201 

. 001231 

.000076 

1.12 

.  00  J6  0  8 

.003056 

. 000929 

.001169 

.000067 

1.16 

.003371 

.002806 

.000738 

. 001055 

.000053 

1.20 

.00-J157 

.  0u2669 

.000533 

.000975 

. 000052 

1.24 

. 00O977 

.  00249' 

. 000398 

. 0u090o 

. 000047 

1 .  2  o 

.003608 

.  00  2261 

. 000262 

.000851 

.000042 

1.  32 

.  008496 

.002355 

. 000164 

. U00800 

.  000037 

1.3b 

.  0  0  a  3  7  2 

.002108 

.000057 

. 000760 

.000033 

1.4-0 

.OOtJlbb 

. 0u208d 

.000019 

. 000725 

.000030 

1.  44 

. Jool29 

.00192  2 

.000024 

.000700 

.  000027 

1.4  8 

.  008102 

.  001986 

. 000027 

. 000670 

.000024 

1.52 

. J07336 

. 0ol8tto 

. 000023 

. 000659 

.000023 

1.56 

.007797 

.  0dl809 

.000027 

.000641 

.000021 

1.6U 

.007795 

.  0U  1836 

. 000024 

. 00061a 

.000019 

Loft 

.0u7654 

.001647 

.000019 

.000602 

.  000018 

1.  be 

.00752*+ 

.  0 ol669 

.000014 

.000539 

.  000016 

1.  12 

.0U7S35 

.  0!il634 

.000008 

. 000579 

.000015 

1.  76 

.007498 

.001618 

.000001 

. 000570 

.000014 

1.8  0 

.007505 

.  00  1541 

J .  0  0  0  0  0  0 

. 000566 

.000014 

l.o4 

.007346 

.00153  0 

J.  000  000 

.000554 

.000013 
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PRElIoTED    defl 

.    RESPONSE 

FUNS     (BY 

RA1  10    MULTIPLICATION) 

C131     AlK 

C<AFT    SITE 

-3 

TiUE 

slESP    0 

RE.SP    1 

^ESP    2 

RESP    3 
.  000430 

RESP    4 
.  000050 

0.  GO 

.002525 

.005408 

.  003041 

.34 

.003072 

.006021 

.303242 

. 000539 

.000053 

•  Jo 

.003572 

.006312 

.U0  3265 

.000591 

.  0U035<+ 

.12 

.  004233 

. 007168 

.  003445 

. 000583 

.  000053 

.16 

.005083 

.00312* 

.  0036^:3 

.000445 

.  000053 

.Zu 

.u05ay3 

.U03642 

.  003684 

.000523 

.000055 

.24 

.006994 

.010052 

.003902 

.000430 

.C00054 

.20 

.  0  0  o  1  3  2 

.010625 

.  003907 

. 000430 

.0000^5 

.  32 

.00*734 

.011759 

. 004138 

.000474 

.000055 

.36 

.U11274 

.012682 

.004251 

. 000437 

.000056 

.<+0 

.012769 

.013098 

.004274 

.000548 

.000056 

.44 

.J14269 

.013575 

. 004332 

.000642 

. 000057 

.4b 

.015795 

.  014492 

. 004473 

.000554 

. 000056 

.52 

.  J16874 

.014583 

.  004539 

.000610 

.000057 

.56 

.  U17500 

.014601 

.004500 

.0  00691 

.00  30  53 

•  UJ 

.  J17710 

.014434 

. 004560 

. 000523 

.000357 

.64 

.U17267 

.013616 

.0043*2 

.000563 

.  000057 

.  b6 

.016360 

.01334o 

.004269 

.000601 

.000053 

.  72 

.U1556C 

.012293 

.u04153 

.000586 

.  000053 

.76 

.  014173 

.01099  • 

.00  3984 

. 000534 

.001059 

.80 

.012358 

.010110 

. 003975 

.000623 

.  300059 

.  Ik 

.011628 

.00  9331 

.  0  0  3  8  8  -J 

.000567 

.000057 

.  00 

.010289 

.008171+ 

. 003596 

.000607 

.000057 

.92 

.00o968 

.0u7556 

.0  03393 

. 000660 

.000057 

.■36 

.  0  0  o  2 1  0 

.00707o 

. 003431 

.00049O 

.000056 

1.  uu 

.  J07Q64 

.005856 

.  U03199 

. 000654 

.000057 

1.  U4 

.006210 

.  0u548  4 

.003050 

.000661 

.000056 

1.  uo 

.005427 

.004721 

.002955 

.000571 

.  000054 

1.12 

.304711 

.0j3877 

.  002708 

.000577 

.000055 

1.16 

.004077 

.003489 

.002674 

.000557 

.000053 

1.20 

.  003452 

.0  33081 

.002449 

. 000557 

.000052 

1.24 

.003106 

.0029  3  7 

.  U02384 

. 000575 

.000052 

1.28 

.002371 

.  0u2723 

.002312 

. 000615 

.000052 

1.32 

.U02539 

.002690 

.002251 

. 000629 

.000051 

1.36 

.002200 

.002566 

. 002160 

.000641 

.000050 

1  .  <+U 

.01L2093 

.0i,256t 

.002000 

.000661 

.000050 

1.  4<+ 

.  001908 

.002044 

.  001893 

.000711 

.000049 

i.  46 

. 001602 

.001996 

. 001794 

. 000693 

.000049 

1.52 

.  0U1654 

.  001748 

.001391 

. 000589 

.000046 

1.56 

.001593 

.001788 

.  001729 

. 00052O 

. 000044 

1.60 

.0ul397 

.001607 

. 001579 

.000613 

.000044 

L.64 

.001358 

.00166  5 

. 001558 

. 000504 

.  000042 

1.6o 

. 001226 

. 001506 

.  0014U1 

.000493 

.  000041 

1.72 

.001060 

.001173 

.  001234 

. 000596 

.  0  00041 

1.  76 

.  J00965 

.  0ul24o 

. 001066 

. 000568 

.0000  39 

l.oO 

.00  j949 

.0Uu926 

.000330 

.000567 

.000037 

l.o4 

.  0U0658 

.000897 

.  0  0  0  6  6  7 

.  0  00567 

.  0000  36 
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PREDICTED  DEFL 

.  RESPONSE 

"UNS  (BY 

RATIO  MULTIPL 

ICATION) 

C135  AIPC 

<AFT  SITF 

-2 

r  i  m  e 

*ESP  0 

RESP  1 

RESP  2 

RESP  3 

RESP  *4 

0  .  00 

.  001697 

.002731 

. 000673 

. 000061 

. 000001 

.  0*4 

.002255 

. 004043 

.  000912 

. 000067 

. 000001 

.  0  ti 

.00  2944 

.005029 

.  00116*4 

.  000076 

.000001 

.12 

.004263 

.  0  0  646  6 

.001419 

.  000086 

.000001 

.16 

.006980 

. 00908  7 

.  001836 

. 000097 

.000001 

.20 

.012261 

.01240  7 

. 002299 

.000104 

. 000001 

.24 

.021146 

.016259 

. 002726 

.000115 

.000001 

.25 

.03397*4 

. 02060  3 

.003248 

.000121 

.000001 

.  32 

.  049279 

.025211 

.003772 

. 000127 

. 000001 

.36 

.060886 

.027038 

.00*4117 

. 000137 

.000001 

.40 

.055969 

.0  252  9  8 

.004198 

. 000140 

.000001 

.  44 

.042979 

.022422 

.004137 

.000139 

. 000001 

.  48 

.  037554 

.02147  7 

.  00409*4 

. 000137 

.000001 

.52 

.042083 

.023455 

.004121 

.000137 

.000001 

.56 

.  057122 

. 027125 

.004230 

. 000130 

.000001 

.  bU 

.  07070b 

.027506 

.00*4231 

. 000127 

.0  00001 

.6*4 

.062311 

.02343  8 

. 004003 

.  00012*4 

.000001 

.66 

.  045106 

.017679 

. 003588 

.  00011*4 

,  000001 

.72 

.03  3421 

.012456 

.  0  [)  3  0  6  7 

.  000107 

.000001 

.  76 

.  024582 

.00768*4 

.  002*4f  0 

.000098 

.0  00001 

.80 

.019405 

.  00*46  0  2 

.  001975 

.  000P88 

.0  00001 

.64 

.016243 

.002505 

.  001626 

. 000077 

.000001 

.aa 

.014236 

.001070 

.001346 

.  000070 

. 000001 

.92 

.  013036 

.  00070  0 

. 001137 

.  00006*4 

.000001 

.96 

.012369 

. 000731 

.  00  1004 

.000062 

.000001 

l.  oo 

.011734 

.000266 

.  0008*4  i 

.  000048 

.000001 

1.U4 

.011320 

.  000548 

. 000726 

.  0000*41 

.000000 

1.0b 

.011185 

.  000772 

.  00067*4 

. 0U0032 

.  000000 

1.12 

.010768 

.000577 

.  000562 

,000028 

.  o  o  a  n  o  o 

1.  16 

.  01039*4 

.000545 

.000442 

. 000026 

.000000 

1.20 

.010170 

.  000281 

.  0  0  0  3  r<  1 

.  000024 

.000000 

1.24 

.  010016 

.000278 

.  0003*46 

. 000017 

.000000 

1  .2o 

.009631 

.000243 

.000233 

. 000020 

.0  00000 

1.32 

.  009461 

.000152 

.000200 

.  0  0  0  0  1 9 

.  000000 

1.  36 

.00^477 

.000379 

. 000217 

. 000013 

.  000000 

1.40 

.  009225 

.  000200 

.  0  0  0 1  <4  3 

. 000017 

.000000 

1.44 

.009100 

.000275 

. 000118 

. 00001 4 

.000000 

1.4ft 

.009399 

.000051 

.000097 

.  00001*4 

.000000 

1.52 

.008925 

.000034 

.000126 

.  00001*4 

.000000 

1.56 

.008858 

.  0uu032 

.000050 

.000013 

.000000 

1.60 

.009008 

.00002^ 

. 00017  J 

.000011 

. 00000  J 

1.64 

.  006925 

.  000035 

.000096 

. 000OU7 

.000000 

1.68 

.008789 

.000016 

.000077 

. 00001 ? 

. 000000 

1.72 

.  G09Q35 

.000266 

. 000101 

. 000007 

.000000 

1.  76 

.009076 

.  00043*4 

.  0  0  0  1  5  :J 

.  0  Q  0  0  0  <4 

.0  00000 

1.80 

.  j0887*4 

.000280 

.000085 

.000009 

.000000 

1.  64 

. J06793 

.000466 

.000068 

.00000/ 

.000000 
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PREDICTED    DEFL 


RESPONSE    FUNS     (BY    RATIO    MULTIPLICATION) 
C130     AIRCRAFT    SITE-2 


TIME 

RES°  0 

RESP  1 

RESP  2 

RESP  3 

RESP  4 

0.  00 

.  001527 

.  0C3158 

. 00  0942 

. 000060 

. 000001 

.  uu 

.001356 

.00  343  8 

. 001052 

.000063 

.000001 

•  U  0 

.002205 

.003576 

.  001107 

. 000067 

.000001 

.  12 

.002812 

.004457 

.001252 

. 000070 

.  000001 

.16 

.003732 

.  005409 

. 001377 

. 000073 

. 000001 

.  20 

. 005067 

.  006777 

. 001539 

.  00007f> 

.0  00001 

.24 

.006975 

.008218 

.0016H7 

.  000080 

. 0  00001 

.26 

.009702 

.010499 

. 001395 

.  000033 

. 000001 

.32 

.013412 

.011868 

.  002122 

.  000036 

. 000001 

.36 

.  018049 

.013735 

. 0023<*5 

.000039 

. 000001 

.  40 

.  023858 

.  015556 

.002463 

.000091 

.000001 

.1*1* 

.030870 

.017913 

.0026C16 

.000093 

.000001 

.48 

.038536 

.019367 

.  002662 

. 000095 

.000001 

.62 

.045155 

.020513 

. 003011 

. 000097 

.0  00001 

.56 

.049667 

.020712 

.  003035 

. 000099 

.000001 

.  60 

.051100 

.020589 

.0030*47 

. 000100 

.000001 

.64 

.049669 

.  019U98 

.003033 

. 000101 

. 000001 

.  bO 

.  045943 

.016932 

.00  30  39 

. 000100 

.0  00001 

.72 

.041072 

.  014749 

.  003045 

. 000101 

. 0  00001 

.  7b 

. 035908 

.012654 

.003020 

.  oooni 

. 000001 

.80 

.031088 

.010674 

.002955 

. 000100 

.000001 

.0"+ 

.027458 

.  009420 

.  002923 

. 000099 

. 000001 

.  8(3 

.025159 

. 009101 

. 002387 

. 000096 

.000001 

.92 

.023958 

.003896 

. 002877 

.000098 

.000001 

.96 

.024441 

.008274 

. 002839 

.000097 

. 000001 

1.  GU 

.  026436 

.003375 

.  002816 

. 000097 

.000001 

1.  04 

.029688 

.  003772 

. 002760 

.000097 

. 000001 

1.  Ob 

.G33">35 

. 008841 

.002772 

. 000097 

.000001 

1.12 

.036143 

.00  8917 

.002747 

. 000097 

. 000001 

1.16 

.037363 

.006865 

. 002721 

.000095 

. 000001 

1.20 

.03dl63 

.003079 

. 002677 

. 000094 

.000001 

1.24 

.037205 

.007592 

. 002606 

.000092 

. 000001 

1.26 

.035122 

.006406 

.002530 

. 000091 

. 000001 

1.32 

.032140 

. 004634 

. 002417 

. 000089 

. 000001 

1.36 

.026447 

.00  257  3 

. 002243 

. 000088 

. 000001 

1.40 

.024682 

.00110  7 

. 002100 

.000035 

. 000001 

1.  44 

.021058 

.000031 

. 001934 

. 000082 

.000001 

1  .48 

.017960 

a.  ooooo  o 

. 001814 

. 000078 

. 000001 

1.52 

.  015145 

0.000000 

. 001677 

. 000075 

. 000001 

1.56 

.U12773 

0. 000000 

. 001540 

. 000072 

.000001 

1.60 

. 010638 

0. 00000  0 

. 001351 

. 000068 

.000001 

1.6*1 

. 006955 

0. 00000  0 

. 001200 

. 000065 

.000001 

1  .68 

.007667 

0.00000  0 

. 001104 

. 000061 

. 000001 

1.72 

.006861 

0. 00000  0 

. 0010  26 

.000057 

.000001 

1.  7b 

. 006110 

0.00000  0 

.00  0  925 

. 000054 

. 000001 

1  .81) 

.005444 

0.000000 

.000835 

. 000052 

.000001 

1  .84 

.  004862 

0.000000 

. 000741 

. 000050 

. 000001 
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PREDICTED  DFFL 

.  RESPONSE 

^UNS  (BY 

RATIO  MULTIPLICATION) 

C135  ^IRC^AFT  SITE 

-3 

TIME 

RESP  0 

RE3P  1 

RESP  2 

RESP  3 

RESP  4 

o 

.  UO 

.001236 

. 00330  1 

.  001892 

.000717 

.  000053 

.04 

.001643 

.004887 

.  002563 

.000791 

.  0  00062 

.  ul  tl 

. 002145 

. 006079 

.003269 

. 000899 

.000071 

.12 

.  003106 

.  0  0  78  1  i 

.  003967 

.001043 

.000080 

.  16 

.  005085 

.01098  5 

.005156 

. 001150 

.000088 

.20 

.  u0d933 

.014996 

.00645d 

. 001228 

.  000095 

.2'* 

.015407 

.  019654 

.  007664 

.001355 

.000102 

.2d 

.024679 

.0  249  0  5 

.  009124 

. 001430 

. 000107 

.32 

.035903 

. 030475 

.  010596 

. 001502 

.000113 

.  36 

.044360 

.0  3268  4 

.011564 

. 001619 

.000119 

.40 

. 040778 

. 030581 

.  0117O2 

. 001650 

. 000123 

.44 

.  U31313 

.027104 

.  011620 

.0  016  41 

.  000126 

.4ti 

.  027361 

.025962 

.  011500 

.  001616 

. 000126 

.52 

.  030661 

. 02d353 

.  Q1157P, 

.0  01616 

.000123 

.56 

.042127 

.0327d9 

.011882 

.001534 

.000120 

.611 

.  051510 

.033250 

.011865 

. 0  01496 

.000118 

.64 

.045396 

.026333 

.  011257 

.001464 

.000114 

.  68 

.033373 

.021612 

. 010079 

.001350 

. 000110 

.  7? 

. u24349 

. J15057 

.008615 

. 001263 

.003105 

.  76 

.01812ft 

.009530 

.0069f'5 

. 00115H 

. 000097 

.60 

. 014138 

.005563 

. 005649 

. 001039 

.000089 

.64 

.  011834 

.003028 

. 004568 

. 000906 

.  0  0  0  0  °  0 

.  80 

.  010372 

.001293 

.  0  0  7  7  d  0 

. 000830 

.  0000  71 

.92 

.  009498 

.  000846 

.  003192 

. 000757 

.  000063 

.96 

.009012 

.  00088  4 

.  002820 

.000613 

.000056 

l 

.  uO 

. 008549 

.  000322 

.002366 

. 000562 

.000049 

l 

.  U4 

.  00  8?48 

.  000662 

.  002039 

. 000488 

. 0000  43 

l 

.  0  8 

.  0  0  0 1 4  9 

.  000933 

. 001893 

. 000775 

.  000036 

l 

.12 

.007860 

.000698 

.001578 

.000336 

.  000032 

l 

.16 

.007573 

.000658 

. 001241 

.000298 

.000027 

i 

.20 

. 007410 

.  000340 

. 001098 

. 00027b 

.  000023 

l 

.24 

. 007297 

.000337 

. 000978 

. 000205 

.000021 

l 

.28 

. 007Q17 

.  000293 

.000654 

.000235 

. 000019 

l 

.32 

.006908 

.000184 

. 000562 

. 000225 

. 000G16 

l 

.36 

.006905 

.000456 

.0  00611 

.000149 

.  000015 

l 

.40 

. 006721 

.000242 

.000401 

. 000197 

.000014 

l 

.  44 

.006630 

.  000332 

. 000331 

. 000161 

.000013 

1 

.46 

.  006846 

.000061 

. 000273 

. 000166 

.000012 

l 

.52 

.006502 

.000042 

. 000355 

. 000164 

. 000011 

l 

.56 

.  006454 

.000039 

.000139 

.000151 

. 000010 

l 

.60 

.006563 

.000034 

. 000503 

. 000129 

.  000009 

l 

.64 

.  00650  3 

.  000042 

.000274 

.000083 

. 000008 

l 

.  t>d 

.006403 

.000020 

. 000217 

. 000152 

.000008 

l 

.  /z 

.  006563 

.000324 

.  0  0  0  2  «  4 

. 000087 

.000007 

l 

.  76 

. 006613 

.  000524 

.  000446 

. 000043 

. 000006 

l 

.  do 

.006466 

.  000338 

.  000238 

.000106 

. 000006 

i 

.84 

.  006407 

.  000563 

.  000247 

. 000084 

.000006 

20U 


PREOICTEO  DE 

:fl.  response 

CUNS  (BY 

RATIO  MULTIPLICATION) 

C130  AlKC 

MET  SITE 

-3 

TIME 

RESD  0 

RESP  1 

RESP  2 

RESP  3 

DESP  4 

0  .00 

,001113 

.  00361ft 

.  002645 

. 0U0703 

.000054 

.  04 

.  001353 

.004156 

.002954 

.  000741 

.000057 

.0(3 

.001606 

.00*4323 

.  003109 

.000716 

. 0  00060 

.12 

.002049 

.  005364 

.  003517 

. 000821 

.000063 

.16 

.  002719 

.00  653  9 

.  003868 

. 000859 

.000065 

.20 

.0U3706 

•0U8192 

. 004322 

. 000901 

.  000068 

.24 

.  005082 

.009934 

.  00**739 

. 000940 

.  000070 

.28 

.  007066 

.012692 

.00532*4 

. 000978 

.000072 

.32 

.009771 

.014346 

.005962 

. 001011 

.  000074 

.36 

.013150 

.016604 

.  0065R6 

.001046 

.  000076 

.40 

.017362 

.  018804 

. 006975 

.001075 

.000079 

.44 

.022491 

.021654 

.007574 

.001097 

.  000079 

.  46 

.028076 

.023411 

. 0080*41 

. 001125 

.000080 

.52 

.032899 

.  024797 

. 008459 

. 001141 

.000081 

.56 

.036186 

.025037 

.008525 

.001168 

.000083 

.  60 

.037230 

.02488  8 

.008559 

.0  0118  2 

.000085 

.64 

.U36167 

.023086 

.008520 

. 001187 

. 000066 

.68 

.033473 

.02046  7 

.008536 

.001182 

. 000087 

.72 

.029924 

.017828 

.  0  0  6  5  5  k 

. 001167 

.  0  0  00  89 

.76 

.026162 

.015296 

.00^46  3 

.001186 

. 000090 

.80 

.022650 

.012903 

.008302 

.001175 

.000090 

.64 

.020005 

.01138  7 

. 008226 

.001170 

.000091 

.88 

.013257 

.  011002 

. 008110 

.001162 

. 000091 

.92 

.017455 

.01075(4 

. 008080 

.001154 

.000091 

.96 

.017307 

. 01000  1 

.  00797*4 

.001148 

.000090 

1.00 

.019261 

. 010124 

.  007910 

.001150 

. 000090 

1.04 

.  U21630 

.01060  4 

. 007809 

.001141 

. 0  0  00  89 

1.08 

.024214 

.  010688 

. 0077R7 

. 001142 

. 000088 

1.12 

.026333 

.010779 

. 007715 

.001139 

. 000087 

1.16 

.027566 

.  010716 

.0076*42 

. 001124 

.000086 

1.20 

.  027804 

.009766 

.  007520 

. 001108 

.000085 

1.24 

.027106 

.009178 

.007320 

.001090 

. 000083 

1.26 

.1125589 

.  007744 

. 007106 

.001069 

. 000082 

1.32 

.023416 

.005602 

.006790 

. 001048 

.000081 

1.  36 

.020725 

.003111 

.006300 

.  001034 

.  000080 

1.40 

.017963 

.001338 

.005900 

. 001002 

. 000078 

1  .44 

.015342 

.0000  3  7 

. 005^33 

.000971 

.000076 

1.  48 

.013085 

0. 000000 

.  00509*4 

.000922 

. 000073 

1.52 

.011034 

0.00000  0 

.  004710 

. 000886 

.000071 

1.56 

.009306 

0.000000 

.  004326 

. 000648 

.  000069 

1.60 

.007750 

0. 000000 

. 003795 

. 000808 

. 000066 

1  .64 

.  006524 

0. 000000 

.  003369 

.000764 

.000063 

1.60 

. 005600 

0.  00G00  0 

. 003101 

. 000718 

. 000059 

1.72 

.004998 

0.00000  0 

.  002916 

. 000676 

.000056 

1.  76 

.004452 

0.000000 

. 002598 

. 000642 

.  00005*4 

1.80 

. 003966 

0.  000000 

.  0023*4*4 

. 000619 

.  000051 

1.84 

.003542 

0.  00000  0 

.  002081 

.000589 

.0  000  49 
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^1 


